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HAUPTAUFSÄTZE 


Zur Theorie der Windkanalturbulenz. 


1. Mitteilung. 
Von W. Tollmien und Manfred Schäfer in Dresden. 


Nach einem Überblick über frühere Arbeiten wird durch Linearisierung 
der hydrodynamischen Differentialgleichungen ein mathematisches Modell 
der Windkanalturbulenz aufgestellt, ohne die übliche Annahme der Isotropie 
zu machen. Für das Randwertproblem der linearisierten hydrodynamischen 
Differentialgleichungen wird eine Integrationstheorie entwickelt. Dabei 
ergibt sich eine Trennung der Strömung in eine Potential- und Diffusions- 
strömung. Mittelwerte dieser Strömungen werden diskutiert. 


1. Geschichtliche Einleitung. 


Bei den Windkanälen für Versuche an Flugzeugmodellen tritt als störende Erscheinung 
die Turbulenz des Luftstromes auf. Da diese die Übertragbarkeit der Meßergebnisse auf die 
Großausführung gefährden kann, hat man schon seit langem die Windkanalturbule nz quanti- 
tativ zu erfassen gesucht. Während man zunächst eine Vergleichsmöglichkeit für die Turbulenz 
in verschiedenen Windkanälen aus der Beeinflussung des Widerstandsbeiwertes einer Kugel 
ableitete, maß H. L. Dryden [1]') als erster direkt die Intensität der turbulenten Ge- 
schwindigkeitsschwankungen mit einer Hitzdrahtschaltung. Inzwischen sind sehr zahlreiche 
und eingehende experimentelle Untersuchungen der Windkanalturbulenz gefolgt. Dabei hat 
man oft, um die Erscheinung zu verstärken und bequem abwandeln zu können, vor die 
Versuchsstrecke Turbulenzgitter geschaltet, die im einfachsten Fall aus einem Rost mit 
quadratischen Maschen bestehen. 


Diese außerordentlich lebhafte Versuchstätigkeit hat ihren Grund darin, daß theoretische 
Arbeiten von G. I. Taylor [2] seit 1935 ein reges Interesse an der Windkanalturbulenz 
selbst — ganz abgesehen von ihrer Wichtigkeit für Modellversuche — erweckt hatten. Die 
Windkanalturbulenz stellt nämlich den einfachsten Fall einer turbulenten Strömung überhaupt 
dar. Die mittlere Strömungsgeschwindigkeit ist hier ja räumlich konstant, so daß keine 
Energieübertragung von der Haupt- auf die Nebenbewegung erfolgt. Ferner ist die Schwankungs- 
bewegung im Mittel gleichförmig über jeden Querschnitt senkrecht zur Hauptströmung. Taylor 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttumsverzeichnis am Schluß dieser Abhandlung. 
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hat in seinen Untersuchungen noch die Annahme der mittleren Isotropie der Nebenbewegung 
hinzugefügt, die allerdings für die Windkanalturbulenz keineswegs selbstverständlich ist. 
Danach sollen sämtliche Mittelwerte, die mit Hilfe der Geschwindigkeitsschwankungen oder 
ihrer räumlichen Ableitungen in einem bestimmten Raumpunkt gebildet werden können, von 
einer Drehung oder Spiegelung des Koordinatensystems, das seinen Ursprung in dem betreffen- 
den Raumpunkt hat, unabhängig sein. Am Ausbau dieser Theorie der isotropen Turbulenz 
war neben ihrem Begründer G. I. Taylor besonders Th. v. Kärmän [3] beteiligt. 

Wir wollen nun einige typische Probleme 
aus der Theorie der isotropen Turbulenz 
herausgreifen. Zunächst fragte man nach der 
Abnahme der Turbulenzintensität mit wachsen- 
dem Abstand von der Turbulenzquelle, sei diese 
nun ein Gleichrichter oder ein eigens angebrach- 
tes Turbulenzgitter. Dabei hat man als Turbu- 
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200 lenzintensität die Quadratwurzel aus dem mitt- 
um leren Schwankungsquadrat einer Geschwindig- 
keitskomponente definiert. Oft nahm man dazu 

a. die Geschwindigkeitskomponente in der Haupt- 

strömungsrichtung, dasich dann die Messungdes 
ai Schwankungsquadrates besonders einfach 
durchführen ließ. In der Frühzeit der Erfor- 
schung der isotropen Turbulenz glaubte man 
a — allgemein ein Abklingen der Turbulenzintensität 
mit dem reziproken Abstand von der Turbulenz- 
quelle theoretisch erschließen zu können, was 
7 


auch mit dem seinerzeit vorliegenden Versuchs- 
material übereinzustimmen schien. Mit fort- 
| schreitender Kenntnishat sich das Erscheinungs- 

A 80 2 7 bild wesentlich verwickelt. Neuerdings behaup- 

REG m KR tet man nach theoretischen Erwägungen von 
Bild I. Abnahme der Turbulenzintensität nach Messungen Ban >. - . 2 
im California Institute of Technology. Pasadena [4] hinter Th. v. Kärmän [3] nur noch, daß das Ab- 
4 Turbulenzgittern mit quadratischen Maschen vonder Weite klingen nach einem Potenzgesetz erfolgt, ohne 
a no daß eine allgemeine Aussage über den Ex 
einer Potenz x" des Abstandes x vom Turbulenzgitter ab. ponenten möglich ist. Auch die Experimente 
Der Exponent 5 nor heuer wiedergegebenen Messungen [4] unterstützen diese Aussage (siehe Bild 1); 
erte zwischen 0,6 und 0,8. ni i : 
ein allgemeines Abklingungsgesetz existiert 
jedenfalls sicher nicht. Dryden [5] stellte sogar ein logarithmisches Abklingungsgesetz auf, 
indem er von einem anderen Gedankenmodell ausging, und konnte dies Gesetz an den 
Experimenten seines Laboratoriums gut bestätigen. 

Sodann interessierte man sich für die räumliche Ausdehnung der Turbulenzelemente 
(„scale of turbulence*), die man auch wohl als Flüssigkeitsballen bezeichnete. Deshalb 
untersuchte man nach Taylors Vorschlag [6] die Korrelation von zwei Geschwindigkeits- 
schwankungen, die zur gleichen Zeit, aber an zwei verschiedenen Orten beobachtet werden. 
Diese Korrelationskoeffizienten hängen also in erster Linie von dem Abstand der beiden 
Meßpunkte ab. Solche Korrelationskurven sind inzwischen oft gemessen worden. Sie erwiesen 
sich als ein tiefgreifendes Mittel zur statistischen Beschreibung der isotropen Turbulenz. So 
konnte G. I. Taylor [2] einen Zusammenhang der durch Reibung dissipierten Turbulenz- 
energie mit der Krümmung dieser Korrelationskurven aufstellen. Th. v. Kärmän [3] führte 
für die isotrope Turbulenz sämtliche möglichen Korrelationen dieser Art auf zwei Haupt- 
korrelationen zurück, die ihrerseits durch eine Differentialgleichung verknüpft sind. Solche 
theoretischen Beziehungen konnten in manchen Fällen experimentell bestätigt werden. Es 
seierwähnt, daß Keller und Friedmann [7] solche Korrelationen benutzten, um das Programm 
einer allgemeinen Turbulenztheorie zu entwerfen, freilich mußten sie dabei gewisse Vernach- 
lässigungen in Kauf nehmen, deren Tragweite umstritten ist. An diese Arbeit anknüpfend, 
hat Millionshtchikov [8] jüngst einen interessanten Beitrag zur Theorie der isotropen 
Turbulenz geliefert. 

Man hat auch das Spektrum [9], [5] der Windkanalturbulenz experimentell untersucht, 
d. h. die Verteilung der Geschwindigkeitsschwankungen auf zeitliche Frequenzen. Beziehungen 
dieses Turbulenzspektrums zu dem mittleren Schwankungsquadrat der betreffenden Ge- 
schwindigkeitskomponente und zu einer Korrelationsfunktion wurden von Taylor [10] unter 
gewissen Annahmen aufgestellt und mit guter Annäherung experimentell bestätigt. Für diese 
Zusammenhänge ist sogar die sonst wesentliche Voraussetzung der Isotıopie entbehrlich. 
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Wir wollen unseren Überblick über die bisherige Erforschung der Windkanalturbulenz 
abbrechen, indem wir für weitere Angaben über die schon sehr ausgedehnte Literatur auf den 
zusammenfassenden Vortrag von G. I. Taylor [11] auf dem Internationalen Mechanikkongreß 
in Cambridge-Mass. 1935 verweisen. In den Verhandlungen dieses Kongresses finden sich 
außerdem wichtige theoretische und experimentelle Einzelbeiträge zur Windkanalturbulenz. 
Im übrigen werden wir selbst in dieser Abhandlung auf manche Ergebnisse der bisherigen 
Forschung noch eingehender am passenden Orte zu sprechen kommen. 


2. Die Problemstellung. 

In allen bisherigen Untersuchungen hatte man stets nur Mittelwerte und (ihre Beziehungen 
untereinander betrachtet. Es liegt nun der Versuch nahe, in den Mechanismus dieser ein- 
fachsten Turbulenz auf dem Wege tiefer einzudringen, daß ein angemessen vereinfachtes 
Modell dieser Erscheinung aufgestellt und in seinem kausal bedingten Verhalten betrachtet wird. 
In Durchführung dieses Programms linearisierten wir die Navier-Stokesschen Differential- 
gleichungen, so daß alle Glieder, welche die Geschwindigkeitsschwankungen quadratisch ent- 
halten, gestrichen werden. Nun ist freilich die Linearisierung der hydrodynamischen Diffe- 
rentialgleiehungen im allgemeinen ein schwerer Eingriff. Dieses Vorgehen scheint aber im 
vorliegenden Fall weniger bedenklich; denn einmal sind die Geschwindigkeitsschwankungen 
durchschnittlich */,, der Hauptgeschwindigkeit, sodann sind keine Wände vorhanden, die wegen 
der Haftbedingung hohe Gradienten der Geschwindigkeitsschwankungen erzwingen würden °?). 
Für die linearisierten hydrodynamischen Differentialgleichungen wird in dieser Abhandlung 
eine neue Integrationsmethode entworfen. 

Als zweites Unterscheidungsmerkmal unserer Untersuchung möchten wir hervorheben, 
daß wir die Annahme der mittleren Isotropie fallen ließen, die für alle bisherigen theoretischen 
Versuche charakteristisch war. Nach den Versuchsbedingungen für die Windkanalturbulenz 
ist von vornherein nieht einzusehen, weshalb die Turbulenz unbedingt isotrop ausfallen muß, 
da die Richtung der Hauptströmung eindeutig ausgezeichnet ist. Der Übergang von der 
Anisotropie zur Isotropie kann allenfalls in größerer Entfernung vom Turbulenzgitter sich 
vollziehen. 


3. Integrationsmethode für das Randwertproblem der nicht-stationären linearisierten 
hydrodynamischen Differentialgleichungen. 

Zur formelmäßigen Darlegung unseres Problems wählen wir folgende Bezeichnungen. 
Die konstante Hauptgeschwindigkeit bezeichnen wir mit «„. In ihre Richtung legen wir 
die »-Achse eines cartesischen Koordinatensystems. Die Schwankungsgeschwindigkeit habe die 
Komponenten « in der «-Richtung, » in der y-Richtung, »w in der z-Richtung: ? bedeutet die 
Zeit. Die Druckschwankung werde mit p bezeichnet. Mit o als Dichte, v als kinematischer 
Zähigkeit liefern die Navier-Stokesschen Differentialgleichungen bei Unterdrückung aller 
Glieder, die quadratisch in den Geschwindigkeitsschwankungen sind: 
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Hinzu tritt die Kontinuitätsgleichung 


) )v dm 
ou dv On 
+3; #; a a a ee Er 
0x dy dz 
Beiläufig führen wir die mittleren Bewegungsgleichungen nach OÖ. Reynolds an. Die 
Mittelwerte werden über genügend große Zeiten genommen und durch Überstreichen gekenn- 
zeichnet. Da die mittlere Geschwindigkeit konstant gleich «,, ist und die Mittelwerte nur 
von x abhängen sollen, erhält man aus den Revnoldsschen Gleichungen 
dp dur dur d um 


.) 
0 s o 5 o a A a re 
0x "dx - 92 - dx 


2) Taylor [12] hat das Auftreten hoher Gradienten der Geschwindigkeit behauptet. Da das Hauptargument 
von Taylor die 4. Ableitung einer empirischen Funktion benutzt, dürfte diese Taylorsche Behauptung einigen 
Zweifeln ausgesetzt sein. 
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Aus den beiden letzten Gleichungen folgt, wenn man das Verschwinden der Geschwindig- 
keitsschwankungen in sehr großer Entfernung von der Turbulenzquelle, also für große x, 
hinzunimmt, 

ur—d, EEE a ee ie ee ee 


Es sei bemerkt, daß wegen der mittleren Gleichförmigkeit der Schwankungsbewegung über 
jeden Querschnitt senkrecht zur »-Richtung gelegentlich nicht nur über die Zeit, sondern 
außerdem räumlich über diese Querschnitte, also über y und z, gemittelt werden wird. 

Es läge nun der Gedanke nahe, das Differentialgleichungssystem (1) durch Einführung 
der unabhängigen Variablen =x—u,„t an Stelle von x zu vereinfachen. Man würde 
dann für die’ Differentialgleichungen des Systems (1) bekommen: 

0u_ dp, Su du u) 
Pi Zr FE Tr IT 


und 2 weitere entsprechende Differentialgleichungen. In der Tat ist man bisher in den 
theoretischen Untersuchungen von einem derartigen Differentialgleichungssystem ausgegangen, 
wobei man manchmal auch die in den Schwankungsgeschwindigkeiten quadratischen Glieder 
mit anschrieb. Zur weiteren Bestimmung des Problems nahm man die Kenntnis des Strömungs- 
zustandes im ganzen Bereich zur Zeit t=0 an und suchte Aussagen über die weitere zeit- 
liche Entwicklung zu gewinnen. Bei dieser Formulierung des Problems als Anfangswertproblem 
ist die Annahme der Isotropie fast eine Selbstverständlichkeit. 

Demgegenüber schien uns nach den wirklichen Versuchsbedingungen eine andere 
Formulierung des Problems angezeigt. Die Windkanalturbulenz wird im wesentlichen durch 
das Turbulenzgitter erregt. Diese Störung der homogenen Strömung betrachten wir für alle 
Beobachtungszeiten als eingeprägt; also werden in der Ebene des Turbulenzgitters <= x, 
Randbedingungen der Strömung für den ganzen Zeitbereich vorgeschrieben. Gefragt wird 
nach der Verteilung der Geschwindigkeitsschwankungen stromabwärts vom Turbulenzgitter. 
Bei diesem Randwertproblem drängt sich die Frage nach dem Übergang von der anisotropen 
zur isotropen Turbulenz unmittelbar auf, da die «-Richtung deutlich gegenüber der y- und 
z-Richtung ausgezeichnet ist. 

Allerdings läßt sich dies Randwertproblem noch nicht mit unseren linearisierten Diffe- 
rentialgleichungen behandeln, da in der Nähe des Turbulenzgitters die Geschwindigkeits- 
schwankungen nicht mehr klein gegen die Hauptgeschwindigkeit «,„ sind. Im übrigen sind 
in der Nähe des Turbulenzgitters weder die Hauptströmung ®) noch die Mittelwerte der 
Schwankungsbewegung gleichförmig über die Querschnitte verteilt. Wir wählen als Ersatz 
ein möglichst verwandtes Randwertproblem. Die Randbedingungen sollen der Strömung nicht 
mehr an der Ebene #—=x, des Turbulenzgitters aufgeprägt werden, sondern an einer Ebene 
x —=const., die schon genügend weit stromabwärts liegt, um mit linearisierten Differential- 
gleichungen arbeiten zu können. An dieser Ebene, die durch «=0 gekennzeichnet sei, steht 
die Strömung in unserem Gedankenmodell sozusagen ständig unter Beobachtung, um die 
Randbedingungen bei @=0 zu erhalten. Außer an der Ebene z=0 sind noch im Unendlichen 
triviale Bedingungen des Verschwindens von Geschwindigkeits- und Druckschwankungen vor- 
zuschreiben. Gefragt ist nach der Verteilung von Geschwindigkeits- und Druckschwankungen 
stromabwärts von der Beobachtungsebene «=0, also im Halbraum x >0. Es ist klar, daß 
bei diesem Randwertproblem die Einführung der Variablen «’ keine Vereinfachung bedeuten würde. 

Um die Integrationstheorie unseres Problems zu entwickeln, denken wir uns der Reihe 
nach die Differentialgleichungen des Systems (1) nach &, y und z differenziert und addiert, 
wodurch man unter Beachtung von (2) 


te ae Ye En ni 
erhält. 
Demnach ist.p eine harmonische Funktion H (x, y,2,t), die im Unendlichen, inbesondere 
bei 2= verschwinden soll. Für diese harmonische Funktion H machen wir mit einer 
Funktion x den Ansatz: 


ee 
p=H=-o tm 


ENT 
Es ist sofort klar, daß man für 9 nur eine harmonische Funktion zu wählen braucht, die im 
Unendlichen verschwindet, um ein H mit den gewünschten Eigenschaften zu erhalten. Die 
Umkehrung aber, daß man allgemein bei harmonischem H auch für 9 eine harmonische 
Funktion nehmen kann, die im Unendlichen verschwindet, bedarf eines Beweises. Erst 
dadurch wird die Allgemeingültigkeit unseres Ansatzes (5) gewährleistet. 


3) Den Übergang von der ungleichförmig verteilten Hauptströmung in die bei uns vorausgesetzte gleichförmige 
Hauptgeschwindigkeit “u hat L. Prandt] [13] theoretisch untersucht. 
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Von der harmonischen Funktion H, der Druckschwankung, setzen wir voraus, daß sie 
bei =» mitsamt ihren ersten und zweiten Ableitungen so stark verschwindet, daß die im 
folgenden vorgenommenen Integral- und Differentialoperationen durchführbar sind. 
Die Differentialgleichung (5) für 9 schreiben wir jetzt in der Form: 
0q og H(z,y.2,®) 


+u = 
d t m Ri) x o 


FE Pe  " 


wobei die Reihenfolge der unabhängigen Variablen in H im folgenden streng beizubehalten 
ist. H wird als Funktionssymbol behandelt, nicht etwa wie p und @ als eine physikalische 
Größe, für die wir stets dieselbe Bezeichnung verwenden, gleichgültig, wie die unabhängigen 
Variablen gewählt werden. H ist, was beachtet sein will, nur für positive Werte seines 
Argumentes definiert. Für t führen wir die neue Variable 


27 


u m 


ein. Dann wird aus (5a) für p (x,y.2,f): 


dpla,y,2,f) 


u 
m dx Pr 


Unter Berücksichtigung der Randbedingung bei = x erhält man: 


Ki & 
E,y,.,t ES 
? Ol \ ’ | 3 ® Um ’ 
oder in den alten Variablen: 
& 

E:: £ Ee— x 
y= HIE,uy,2,t- Me a, | 
e OUm \ | s it Um 


Um zu beweisen, daß x harmonisch ist, führen wir in (6) die Variable 


el x 


u m 


ein, wodurch wir 

u 
57T E 
\Hia+ ut, BEER: er 


Es 


qJ = 


erhalten. Diese Lösungsform läßt sich auch unmittelbar durch Einsetzen in (5) mitsamt den 
Randbedingungen leicht bestätigen. Aus (7) folgt sofort, daß Ay mit 14H verschwindet, 
vorausgesetzt, daß die zweiten Ableitungen von H beim Unendlichwerden des ersten Argumentes 
genügend stark verschwinden, wie wir schon eingangs betonten, so daß die Differentiationen 
unter dem Integralzeichen ausführbar sind. 

Wir tun der Allgemeinheit der Lösung also keinen Abbruch, wenn wir g als harmo- 
nische Funktion ansetzen, die bei =» verschwindet. 

Die Geschwindigkeitskomponenten «, v, w setzen wir aus je 2 Teilen up, up bzw. vp. vn 
und wp,wp zusammen. Für die ersten Teile können wir 

04 04 04 


"p= P: vwp= (8) 
ray 


Eye 02 

nehmen. Diese Partikularlösung genügt (1) und (2) und der Randbedingung des Verschwin- 
dens im Unendlichen. Der erste Strömungsteil ist demnach eine Potentialströmung. Diese 
Potentialströmung ist bestimmt, wenn wir neben der schon erwähnten Bedingung im Unend- 
lichen @ bei z=0 als Funktion von £, y,z vorschreiben. Unser Ansatz (5) erfährt nun nach- 
träglich seine physikalische Deutung als linearisierte Bernoullische Gleichung für nicht- 


stationäre Potentialströmungen. 


4, Auffallend an der Formel (6) ist, daß % in einem bestimmten Zeitpunkt durch das Verhalten von H für spätere 
Zeiten bestimmt wird. Allerdings ist das Verhalten von H für sehr große Zeiten wieder gleichgültig, da H für große 
Werte des ersten Argumentes stark verschwindet. Immerhin sieht man, daß die Formel (6) wenig zur Berechnung der 
Größe 9 aus H geeignet ist. Wir brauchen sie auch nur zur Herleitung eines allgemeinen Satzes. Die erwähnte 
merkwürdige Tatsache unterstreicht nochmals, daß wir es nieht mit einem Anfangswertproblem, sondern mit einem 
Randwertproblem zu tun haben. 
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Die restlichen Anteile «wp,”p.Wwp müssen dann den folgenden Differentialeleichunzen 
D, UD, Wı 


genügen! 
0 u Ö u | g 
u =? u nr BREI Br .. WIR, 
h) } m d x D 
Or 327) 91 
Un v " a A ee ren a 
a ”" dx D 
Day? mp | g 
ru =r Aw ee u Se 
hi N m Rh) 2 D 


up , den, dwp & (10) 
E > E a ea 2 2 Are 

Wir können dies homogene Differentialgleichungssystem lösen, indem wir aus den beiden 
nichtgekoppelten Differentialgleichungen (9b) und (9e) v, und ww, bestimmen. Zur eindeutigen 
Festlegung von vr) bzw. mw) haben wir dabei neben der Bedingung des Verschwindens im 
Unendlichen v) und wpn bei «0 als Funktion von f, y und z vorzuschreiben. Durch den 
Ansatz 

E37 


on Bi 
Otp cewp\ 
u + a 
D (a, 2, 


s 
werden die Kontinuitätsgleichung (10) und die Randbedingungen im Unendlieben befriedigt. 
Gleichzeitig wird dadurch auch die letzte noch ausstehende Differentialgleichung (9a) erfüllt: 
denn nach (11) wird 


x 
177) (2777) \ rn, Owl dep Om "op , ’arp 
ru —- du r de u " + H 
N mr D i loto Y otodzf m N y m 2 \drd y or0z 
I 
x , i 
| (0° "p 0° mp 0° "p 0° | 1 
- y — — — „+ : dr 
FU NER y oy?0z dydz 02° | 
L 
v % 
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Dabei ist die Existenz der Integrale (11), sowie mehrfach ihre Differenzierbarkeit unter dem 
Integralzeichen vorausgesetzt. Von dem Zutreffen dieser Voraussetzungen (genügend starkes 
Verschwinden von vn und wp bei e=x) kann man sich nach der tatsächlichen Herstellung 
der Lösungen überzeugen. 

Über eine allgemeine Integrationsmethode hinaus erhielten wir eine physikalisch inter- 
essante Scheidung der Strömung in einen Potentialanteil und einen Teil, den wir Diffusions- 
anteil nennen wollen, da die Differentialgleichungen (9) für einen Diffusionsvorgang, aller- 
dings mit Konvektion, charakteristisch sind ®). Die Bezeichnung ‚Reibungsströmung“ wäre 
mißverständlich Für diese Trennung der Strömung in Potential- und Diffusionsanteil wurden 
dabei nur Bedingungen starken Verschwindens im Unendlichen für Druck und Geschwindig- 
keitskomponenten benutzt. Trotzdem das Differentialgleichungssystem (1). (2) häufig diskutiert 
ist, insbesondere von Oseen [15], so ist uns die vorstehende Integrationsmethode aus der 
Literatur nicht bekannt. Die Übertragung auf Fälle mit anderen Rand- oder Anfangsbedin- 
gungen scheint leicht möglich ®). 

Die zur vollständigen Bestimmung unseres Problems bei #0 vorgeschriebenen Rand- 
verteilungen lassen sich auch leicht physikalisch deuten. Einmal wird 9=g, gefordert. 
Diese Potentialverteilung bei #«—0 hängt aber unmittelbar mit der dortigen Verteilung der 
Druckschwankungen p, zusammen, vgl. etwa die später aufgestellte Formel (42). Nach dieser 
Festlegung des Potentialanteiles der Strömung hat man auch #> = rp, und ırp = p, bei z — 0. 
Der Rest der Querschwankungen bei z—=0: @,— rp, bzw. w, — mp, muß gleich den Rand- 
werten der Diffusionsströmung "7, und ıp, sein. Die Vorgabe der drei Randverteilungen 
Po, ® Do: Wp, kommt also darauf hinaus, daß die Druck- und Quergeschwindigkeitsschwankungen 
bei 20 vorgeschrieben sind. 

5) Die Trenunng der linearisierten Strömung in Potential- und Diffusionsanteil wurde zuerst 1938 in einer 
Diskussionsbemerkung [14] auf dem Mechanikkongreß in Cambridge-Mass. von dem ersien Verfasser dieser Abhand 
lung mitgeteilt. 

‘) L. Prandtl [16] hat neuerdings aus Messungen von H. Reiehardt im geschlossenen Kanal vermntet, 


daß die Schwankungensich aus einem zähigkeitsunabhängigen und aus einem zähigkeitsabhängigen Bestandteil aufbauen. 
Diese Vermutung steht in Analogie zu unserem obigen Satz von der Trennung der Schw ankungsströmung in Potential 


und Diffusionsanteil, 
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In einer Abhandlung von F. Reißner [17] wird die Differentialgleichung 


RT] jr u. u, du | 
N) t Q x: Ö y hi) 2? 


betrachtet, der die Geschwindigkeitsschwankung « gehorchen soll. Zur Zeit *=0 wird die 
Anfangsverteilung von « im ganzen «-y-z-Raum vorgeschrieben. Wegen der angenommenen 
Isotropie braucht nur eine Komponente der Geschwindigkeitsschwankung betrachtet werden. 
Es wird sodann unter anderem die zeitliche Abhängigkeit des räumlichen Mittelwertes von 
u? berechnet. Die Reißnersche Arbeit beschäftigt sich also mit jenem Teil der Strömung, den 
wir als Diffusionsanteil bezeiehneten. Die Herkunft dieser Diffusionsströmung aus den line- 
arisierten Navier-Stokesschen Differentialgleichungen ist bei E. Reißner nicht aufgezeigt. 
Da weder Druckgefälle noch Kontinuitätsgleichung berücksichtigt sind, handelt es sich bei 
diesem in seiner Einfachheit bemerkenswerten Gedankenbild der isotropen Turbulenz eigent- 
lieh nieht um ein hydrodynamisches Modell. 


4 Mittelwertberechnung nach dem Parsevalschen Satz für Fourierintegrale. 
Sätze von G. I. Taylor und J. Kampe de Feriet über das Turbulenzspektrum 
und Meßergebnisse. 

Diebeiden klassischen Rand wertprobleme, aufdie unser Integrationsverfahren führte, werden 
wir teilweise nach der Greenschen Methode, teilweise nach der Fourierintegralmethode be- 
handeln. Für die Abhängigkeit von der Zeit werden wir stets ein Fourierintegral ansetzen; 
einmal, weil die Verteilung der Geschwindigkeitsschwankungen auf zeitliche Frequenzen, das 
„[urbulenzspektrum*, auch experimentell untersucht ist und wir einen Zusammenhang zwischen 
Turbulenzspektrum und Abklingen der Turbulenzintensität feststellen werden, zum anderen, 
weil die Berechnung der zeitlichen Mittelwerte durch den Parsevalschen Satz für Fourier- 
integrale stark vereinfacht wird. 

Da verschiedene Bezeichnungen in der Theorie der Fourierintegrale nebeneinander 
üblich sind, setzen wir zunächst eine bestimmte Bezeichnung fest. FA) ist die Fourier- 
transformierte einer Funktion f(f), wenn 


=SFWwe"d... u RE: 
gilt”). Dabei ist 
Fi) s.\fihe ie RR? (12a). 


Für eine reellwertige Funktion f(f) ist der Realteil von F (4) eine gerade Funktion von 4, der 
Imaginärteil eine ungerade. Sodann führen wir die Parsevalschen Formeln an. Es seien 
F(i) und @ (4) Fouriertransformierte von f(f) bzw. von g(f). Dann ist 

L x + % 

ae Erd: :.: . 2er ©. 

- L 
wobei durch eine darübergesetzte Tilde () der Übergang zur konjugiert komplexen Funktion 
angezeigt werden soll’). Für g(f)  f(f) ergibt sich die Parsevalsche Formel im engeren 
Sinne: 


& +. 
re EUR . >... SEE . (14). 
Es "o 
Ferner sei der Faltungssatz erwähnt. Zu der durch „Faltung“ entstandenen Funktion von E: 
+ 


\figit-wdu 


gehört die Fouriertransformierte 27 F(4) G (4). Für die Gültigkeit dieser Beziehungen ist 
hinreichend, daß fund y für absolut sehr große ft genügend stark verschwinden, was wir bei 
den Funktionen unserer Anwendungen stets voraussetzen können. 
Die entsprechenden Formeln für Funktionen von 3 Veränderlichen f, y,z werden zum 
Teil auch benötigt. Fi, u,,4,) ist die Fouriertransformieite von fit,y, 2), falls 
it 


fit,y.2) = J | I Fin.) di, - -» -» (EB) 


z X Rn 


7) Allgemein werden Fouriertransformierte mit großen Buchstaben, die ursprünglichen Funktionen mit den 
entspreehenden kleinen Buchstaben bezeichnet werden. 


», Man beachte, daß 6 (#) nieht die Fouriertransformierte von gt), sondern von g{ t) ist! 
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eilt. Dann ist 
| | % ' . . . 
2a) \ \ \ fit,y,2)e-WHt+timytimddtdydz . . . . -(16). 


Fi, u,, u,) 


E. Es u. 


Die allgemeine Parsevalsche Formel lautet jetzt: 


2 + + 
\ \ \fit,y,2)g(t,y,z2)dtdydz | 
SL e. L- = 
(17): 
+o +0 + r . 
(rn) \ \ \ Fir, n,.u,) @lA,u,,u,)didu, du, | 
“ 2- x j a 
die spezielle Parsevalsche Formel ist 
1877. too 4% += Lo + 
\ \ IIft,yw, 2)’ dtdydz=(2n) \ \ \FÜ, nu.) didu, du, . (18). 
x I 2 — —. @ 2 


Die Parsevalschen Formeln werden wir zur Berechnung von Mittelwerten verwenden. 
Gedanklich würde es am meisten befriedigen, wenn der Mittelwert von, sagen wir, «? definiert 
wäre durch: 


Hr 

ni ; LI 

“= Um 2 DEREN ea nt 7 U 
PR 


Nun würde eine solche Definition höhere Hilfsmittel aus der neueren Theorie der Fourier- 
integrale erfordern Vor allem entspricht diese Definition durchaus nicht dem Vorgehen beim 
Experimentieren. Deshalb nehmen wir bei unserer Mittelwertdefinition nicht den Grenz- 
übergang zu unendlichem 7 vor, sondern bleiben bei einer endlichen, wenn auch sehr großen 
Beobachtungszeit T. Da wir aber für die Methode der Fourierintegrale den zeitlichen Ver- 
lauf für alle Zeiten brauchen, setzen wir die betreffenden Funktionen, oben also «?, außerhalb 
der Beobachtungszeit 0; in Formeln also 





Ya 
"=z\wdt ET BR EEE SA 
wobei 
. T T 
“=0 für > und I<—z. 


Als erste Anwendung dieser Überlegungen bringen wir einen Beweis zweier Sätze von 
G. I. Taylor [10] über den Zusammenhang des sogenannten Turbulenzspektrums mit der 
Fourierintegraldarstellung der Geschwindigkeitsschwankung « und einem Korrelations- 
koeffizienten, welcher die räumliche Ausdehnung der Wirbelgebilde charakterisiert. Da diese 
beiden Sätze zu den wenigen gehören, welche über die Windkanalturbulenz ohne die An- 
nahme der Isotropie vorhanden sind, dürfte eine vereinfachte Ableitung mit leicht zugäng- 
lichen Beweismitteln willkommen sein. 

Betrachtet werde zunächst die Geschwindigkeitsschwankung « an einem festen Raum- 
punkt, so daß « eine Funktion der Zeit # allein ist: 


a ee a an we 2 A 
Die Fouriertransformierte von f(f) sei F(A). Bei der Messung des Turbulenzspektrums 
wird bestimmt, welcher Beitrag zu «u? von Schwingungen der verschiedenen zeitlichen Frequenzen 
n geliefert wird. Der Beitrag von den Frequenzen zwischen n und n+dn wird mit 
uw q(n)dn angesetzt, so daß die zu messende Funktion q(n) der Bedingung 


ao 


5 ein a ar 


0 
genügen muß. Zwischen der Kreisfrequenz i und der Frequenz n besteht der Zusammenhang 
ae sata wre KIISE 
Nach dem Parsevalschen Satz (14) ist 


= +  “ 
we T=Sf}dt=2a)|FW’di=4m’||Fean)’dn .... . (MM). 


R ao — © 
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Da f(t) eine reellwertige Funktion ist, wird im Anschluß an eine frühere Bemerkung |F(A)? 
eine gerade Funktion von /, so daß wir das Integrationsintervall von 4 bzw. n nur von 0 
bis > gehen lassen können: 
_. 
Br. m = 
= lIFWamPÄan. . . ».... + (D). 
u 


Sn’|Fi2nn)” 


Demnach ist dn der Beitrag, den die harmonischen Komponenten zwischen n 


und n+dn zu «? liefern, also gleich «® q(n) dn. Der gesuchte Zusammenhang zwischen 
Turbulenzspektrum und Fourierintegralzerlegung ist daher 


Sn?|Fean)’ _ - 

T RE 54 
Da |F(2an)|’ eine gerade Funktion von n ist, muß die Spektralfunktion q(n), die zunächst 
nur für positive » definiert ist, für negative n als gerade Funktion fortgesetzt werden, damit 
die Beziehung (26) allgemein gültig ist. 

Jetzt soll in die laufende Auseinandersetzung noch der Korrelationskoeffizient einbezogen 
werden, der zwischen zwei gleichzeitigen Geschwindigkeitskomponenten u gilt, die an zwei 
in der x-Koordinate verschiedenen Plätzen beobachtet werden. Streng genommen, müßte man 
dafür etwa, wenn man von @=0 ausgeht, 





u(0,y,2,8) u(e,y,2,®) 





VYw(O,9,2,1) Vulz,y,z2,t) 


nehmen. Statt dessen nimmt man wegen der kleinen «, für die man diese Größe überhaupt 
von O verschieden findet, mit guter Annäherung 





ud,y,2,t) ule,y,z,P) 


ux (27). 


wi0,y,2,®) 

Wir setzen wie vorhin #(0,y,2,4)=fit). Taylor macht nun die einfache Annahme, daß 
das zeitliche Schwankungsbild unverändert mit der Hauptströmung u, stromabwärts trans- 
portiert wird, jedenfalls über die kleinen Entfernungen, die bei der Bildung der Korrelationen 
in Frage kommen. Folglich soll 


u(2,y,2,®) rlt — De a er Er I 
m; 
sein. Damit wird 
BrR=\f@rlt-—\dt ... 2.2.2222... 
. m’ 


Ist nun wiederum F(A) die Fouriertransformierte von f(f), d..h. 
+% 
fib=|\Fieittda, 


e. 


so folgt, indem man ? durch { “_ ersetzt, 


Im 
+» 
» KL [ PT —i4 2: ; . ‘ 
fit )=\ Fine me u > A HE 54 | ° 
Um 
L 


ve ru k . } | B r 
m.a. W., F(i)e”'"wm ist die Fouriertransformierte von rlt ). Nach dem allgemeinen 





Um 
Parsevalschen Satz (13) ist daher 
+. =“ 
“ a a 0 Pe \ R >> 24 . 
uw TRx: 27 FiÜ)F(Ü)e mdi=2n \ F (A)? e "m nt eh a 


Da | F(i)|’ eine gerade Funktion von 4 ist, kann man auch schreiben: 


n 
FEN r ara AxX . 
u? 1 Rk=4a\ \F (A)? cos di, 
J Um 
0 
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oder mit » statt 4: 

I 
E ») . 
— ,n r R. e e ETN a 
“„TR,=8r \ Fi2an)|? cos dn 
. m 
u 


Mit der Beziehung (26) zwischen 


os 
i { an 
R.=\4 (n) cos An 
. m 
0 
Umgekehrt gilt 
Es 
D) 
anna 
yn)= \ Rx cos dx. 
Um . Um 
) 


R, und q sind demnach wechselseitig Fouriertransformierte. 
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Bild 2. Turbulenzspektrum und Korrelation Ry nach Messimgen von L.F.G. Simmons 
[10] im Abstand von !/3 Maschenweiten hinter einem Turbulenzgitter aus Flachstäben 
mit einer Weite der quadratischen Maschen von 3 Zoll. 





F\? und der Spektralfunktion y ergibt sich schließlich 


133). 


(34). 


Diese Beziehung ist 
experimentell recht 
gut bestätigt an Mes- 
sungen von Sim- 
mons und Salter 
[9] (siehe Bild 2), fer- 
ner an Messungen 
vonDryden[5]und 
seinen Mitarbeitern. 
Offenbar ist dieser 
Satz nicht einmal auf 
die Windkanalturbu- 
lenz beschränkt, da 
ja nur eine gewisse 
Erhaltungstendenz 
des turbulenten 
Schwankungsbildes 
vorausgesetzt wird. 
So komnten Rei- 
chardt und Motz- 
feld [18] seine Gül- 
tigkeit auch im recht- 
eckigen Rohr nach- 
weisen. 
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Korrelationen 
National Bureau of Standards, Washington [5] im Ab- 
stand von 251/; Maschenweiten hinter einem quadratischen 
Turbulenzgitter der Maschenweite 1 Zoll bei einer Wind- 


25 B/ 7 


33mm %0 


Rx und Ry im 


Pf T | Bild 3 
| ı (iımks). 
| link 
| | | Iz | | Turbulenzspek 
| I T | trum nach Messun- 
| | | gen von L. F.G. 
z | | | Simmons|[I0] im 
7 — — to 
q TI | | 1 | höheren Frequenz- 
| | | | | | bereich. 
_b— .- 
I | 
N | | 
RE r——tH- 
3 | | 
| | 
| | 
— 1 
| ) 70 15 
| Am 
| 
Me are — 
| Bild 4. Messungen der 
| | 
27) 60 
4213 
geschwindigkeit von 40 feet/sec. 
Sehr bedeutsam scheint für die Windkanalturbulenz die in diesem Zusammenhang 


getroffene Feststellung G. 1. 
der Geschwindigkeit u,,. 
eine Funktion von allein. 
m 


(Geschwindigkeiten «,, im großen und ganzen ähnlich bleiben. 


Taylors über die Abhängigkeit der Spektralfunktion q von 
Nach den Experimenten ist «„qg fast im ganzen Bereich von n 


Dies bedeutet, daß die Strömungen bei allen untersuchten 


Praktisch (etwa zu 99"/,) ist 


7 


DE ee ee 





XUM 


ar re 


7 





sw ıaAA 


Z. angew. Math. Mech. 2 - { ae a . . “ 
Bd. 21 Nr.1 Feb. 1341 lTollmien u. Schäfer, Zur Theorie der Windkanalturbulenz 11 





die ganze Turbulenzenergie auf harmonische Komponenten zurückzuführen, deren mittlere 
Verteilung ähnlich bleibt (siehe Bild 2). Der nicht ähnliche Rest des Spektrums (siehe Bild 3) liegt 
beiden hohen Frequenzen; er trägt, wie gesagt, zur Turbulenzintensität außerordentlich wenig bei, 
dagegen sehr viel zur Dissipation (in manchen Fällen 73°%,). Schließlich hat Dryden [19] 
das Turbulenzspektrum in 2 verschiedenen Abständen vom Gitter gemessen. Danach scheint 
die Energie aller Frequenzen ungefähr gleich stark abzunehmen. 

Aus dem raschen Abklingen von R, mit = konnte J. Kampe de Feriet [20] den 
wichtigen Schluß ziehen, daß das Turbulenzspektrum ein linienfreies kontinuierliches Spektrum 
darstellt. Ferner machte J. Kampe de Feriet über die Spektralfunktion q(n) die Aussage, 
daß sie bei »n—=0 nicht verschwindet. Während in der Originalfassung der Beweis dieser 
Sätze mittels der Summenfunktion von q(n) und Fourier-Stieltjesscher Integrale geführt 
wird, wollen wir uns hier mit unseren einfachen Hilfsmitteln Einsicht in den Sachverhalt 
verschaffen. Wären Linien in dem Turbulenzspektrum vorhanden, so müßte die für kontinuier- 
liches Spektrum definierte Spektralfunktion q(n) an einzelnen Stellen unendlich werden. 
Nach der Formel (34) tritt dies wegen des experimentell festgestellten raschen Abklingens 
von R, mit. aber nirgends ein, also ist kein Linienspektrum zu erwarten. Dieselbe Formel 
liefert 

I 
40) = , EEE 
Um. 
u 
Da AR, nach allen bisherigen Messungen für die Windkanalturbulenz sogar stets positiv ist), 
muß entsprechend der Behauptung von Kamp« de Feriet 4(0) von O verschieden sein. 


5. Die Potentialströmung und einige ihrer Mittelwerte. 

Nach diesem Zwischenspiel, das uns mit wichtigem Erfahrungsstoff bekannt gemacht 
hat, wenden wir uns zur Lösung der klassischen Randwertprobleme, auf die wir unser 
hydrodynamisches Problem zurückgeführt haben. Zunächst stellen wir den Potentialanteil 
der Strömung her, d.h. eine Lösung g von 


lq N a a N EEE A 
im Halbraum & >0 mit den Randbedingungen 
y=Yg (Yy,2, 0) bei Pin 0 
und den Bedingungen des Verschwindens im Unendlichen. 
Die zeitliche Fouriertransformierte von 9, (y,2,f) sei P,(4. 9.2), d.h 
r x 
ON ih * EEE | ° 
F. 
Die zeitliche Fouriertransformierte P(i,x,y.z) von p(x#,9,2,#) bestimmt sich dann nach 
dem Poissonschen Integral zu 


+0 + 


P=5,\ \06,0 


L L 


dndc 
- weg > 


+ — nf +le — Sy 


Daraus lassen sich die Geschwindigkeitsbeiträge der Potentialströmung und insbesondere der 
Druck berechnen. Abschätzungen der mittleren Quadrate dieser Geschwindigkeitskomponenten 
wurden mit der Sehwarzschen Ungleichung vorgenommen. Da diese Diskussion aber 
bisher nicht zu interessanten Resultaten geführt hat, wollen wir uns gleich der zweiten be- 
nutzten Lösungsmethode zuwenden. 
Der Randverteilung 9,(y,2,f) wird jetzt als zeitlich-räumliche Fouriertransformierte 
D, (A, u,. 4,) zugeordnet, d.h. 
I + X + u 
923, 0)=\ \ \ Di, u,,,) dttimytrimzdidudu, - - . . (39. 
: I 2 — 
Die zeitlich-räumliche Fouriertransformierte P iA, u,, ,) von p(y,z,#) ergibt sich zu 
Ba EP TR 
°) Dagegen kaun Ry, der Korrelationskoeflizient zwischen zwei gleichzeitigen Geschwindigkeitsschwankungen u 


an zwei Punkten, die um y quer zum Strom auseinanderliegen, nach manchen Beobachtungen auch schwach negativ 
werden, vgl. etwa [10], Bild 4 oder [11], Bild 2. 
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oder auszeschrieben, 


+ + u +] 
q \ \ \ D, (4 4 u 
L .. u. 


, t,)e Vu? +ug2a rikt+imyrisezdidu,du, . . . (40a). 


Die zeitlich-räumlichen Fouriertransformierten der Geschwindigkeitskomponenten der 
Potentialströmung sind 


Up Yu?+ u? D(k, u,, u,)e= VYartns?x, | 
„= iu, D,(i, u,, u,) e-Vru®?+tus?x. (41) 
Wr= Nu, P,(i. u,,u,)e Vu? +us?x | 
Die zeitlich-räumliche Fouriertransformierte des Druckes ist 
P=-oliäi-unVur te Dlü,u,u,J)e-Yurtutx . . .... (4). 


Bei der Berechnung der Mittelwerte nehmen wir an, daß die Randverteilung 9, (y,2,!) für 
il> -i y\> £ 2 u identisch O0 wird. Wir setzen deshalb den quadratischen zeit- 
lich-räumlichen Mittelwert von up 
a erg 
up = FrZ\ \ \urdtayaz a ee ae selig, 5 MER 


—_ 2 E- Es 
vgl. dazu unsere früheren Bemerkungen über Mittelbildung bei der Gleichung (20). Gewisse Beden- 
ken gegen diese Mittelwertdefinition sollen für den Augenblick zurückgestellt werden. Nach dem 
Parsevalschen Satz (18) wird aus dieser Mittelwertdefinition 








er 
| i EURER 
pr yzea®\ \ \ ran. ee ee: 
eh - x 2 
Führen wir 
tr = 
eu... ern 
ee) 
ein, so hat man schließlich 
+ x r 2 
"pP= TYy ze \ \ (u? + u) Pr lu,, u,)e”” Ver? + us? x du, du,, 
i u 
Es I 
pP=o yz' a" \ \ We D*lu,,n,e?Vmttns? x du, du,, EN... \ 
"x u 
+2 + 
nPonyy 2 | \ Ru DP*lu,.u,)e ? Ve®rus®x Au, du, 
L 
Aus diesen Gleichungen folgt sofort die merkwürdige Beziehung 
a u © SORTE NEE See er ;' 


Dieser Satz ist zunächst für die zeitlich-räumlichen Mittelwerte abgeleitet, d.h. die Mittelungen 
sind über die Zeit und den ganzen (Juerschnitt (senkrecht zu x oder w,„) genommen. Da 
aber die rein zeitlichen Mittelwerte der ganzen Problemlage nach im wesentlichen konstant 
über den ganzen (Juerschnitt sind, gilt dieser Satz auch für die rein zeitlichen Mittelwerte. 
Weil y vor z nicht irgendwie ausgezeichnet ist, wird außerdem vp? = #p? sein. Der qua- 
dratische Mittelwert der Längsschwankung der Potentialströmung ist also 
das Doppelte des quadratischen Mittelwertes einer Querschwankung. Die 
Potentialströmung kann also für sich allein nicht isotrop sein. 

Wie schon erwähnt, kann man gegen unsere Mittelwertsdefinition einen Einwand 
machen. Bei dieser Definition war vorausgesetzt, daß die zu mittelnden Größen im wesent- 





lichen nur innerhalb der Querschnittabmessungen Y, Z und innerhalb des Zeitabschnittes 7 


von 0 verschieden sein sollen, während sie außerhalb dieser Grenzen praktisch verschwinden 
sollen. Freilich kann man diese Eigenschaften in der Randebene x = 0 erzwingen, aber für 
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größere x breitet sich die zunächst auf - Zey=5, >>2=5 beschränkte Erregung 


seitwärts aus. Eine zeitliche Ausbreitung spielt bei der Potentialströmung keine Rolle, da 
Ip=t0 die Zeit nicht enthält. 

Um einen Begrift zu bekommen, inwieweit diese seitliche Ausbreitung die Mittelwerte 
und ihr besonders interessantes Abklingen mit x verfälscht, studieren wir ein ganz einfaches 
Beispiel. Wir nehmen zur ersten Orientierung sogar eine ebene Potentialströmung in der 
x-y-Ebene und betrachten die Mittelwerte von up’. Da up ebenfalls der Laplaceschen 
Potentialgleichung genügt, können wir leicht die Verteilung von u» anschreiben, die bei x = 0 für 
= LER <= bestimmte Randwerte, sagen wir up,=sinay, außerhalb den Wert 0 ebenso 
wie im Unendlichen annimmt. 

Die gesuchte Lösung stellen wir am besten durch Superposition her. Als Elementar- 
lösung nehmen wir dasjenige Potential, das bei #—=0 für positive y die Werte up,=sinay 
annimmt, dagegen für negative y verschwindet. Wir erhalten für diese Elementarlösung nach 
dem Poissonschen Integral für die Halbebene 


Es 
x sinan 
tp \ r = _— dn, (48). 
a) +iy— 9) 
v0 
oder mit y- ,=n, 
.. 
x | sinaln — My) 
“er=— ä "aaa. ee Bee u nr 
, n \ x +n’ . 


U 


Eine Abschätzung von up läßt sich leicht für „<0 mit Hilfe einer alternierenden Reihe 
finden. Der Zähler des Integranden hat abwechselndes Vorzeichen und ist periodisch, während 
der Nenner stets zunimmt. Zerlegt man das Integrationsintervall in solche Abschnitte, in 
denen der Integrand einheitliches Vorzeichen hat, so erhält man dadurch eine Darstellung 
des Integrals durch eine unendliche Reihe von Integralen mit abwechselndem Vorzeichen. Ab- 
gesehen vom ersten dieser Integrale, dessen Intervall bei y beginnt, nimmt der absolute Wert 
dieser Integrale ab. Eines dieser Integrale ist etwa 
Hi a 
= (sinaln— 9) 
— \ : Be 


nm) a? 


Yi 


a® rn 
wobei sina(y —y)=0 ist. Da y; negativ ist, hat der Nenner als kleinsten Wert «+ yj?. 
Ersetzt man den Nenner durch diesen Wert und integriert den entstandenen Ausdruck, so 
bekommt man einen absolut größeren Wert. Als Absolutschranke dieses Integrals bekommt 
> 
x - R: ee 
man also - >. Da Iy;'=y. so erhält man schließlich als Absolutschranke der alter- 
an(x” + yi?) ; 
nierenden Reihe, durch die wir up darstellen, 
) 


Pe 1 


"pi< (O0). 


anni? — y) 
Für 4 >0 läßt sich die Abschätzung ganz ähnlich durchführen, wenn man das Integral 
für up zunächst zerlegt: 
+» +» 
x (sina(y—n) "sin a(y-- n) 
Nee \ dı 


I 
= ——- dn —— = r ‘ 
x’ +9? / n \ x” +1? . 


— & Y 


Das erste Integrai hat den exakten Wert e”-“*sinay, während das zweite wieder mittels 


) 
. . . R. . . - u 
einer alternierenden Reihe abgeschätzt wird und die Absolutschranke ae top) hat '®). 
anrnıX Tr 4 . 
' . R i a Y Y i 
Diese Elementarlösung wird zweimal mit nach — 5 bzw. + verlegtem Nullpunkt 


angesetzt, sodann wird die bekannte Lösung e=-“*sinay für die Randbedingung u, = sin ay 
10) Eine andere naheliegende Abschätzung dieser Elementarlösung up erhält man nach (48), wenn man die 

“ 1 1 R ; ’ ; r : 
Absolutschranke 1 des sinan, einsetzt: up <<, + „ arctg z . Diese Abschätzung ist bei kleinem x und kleinem 4 


besser als die obenstehende. 
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für -e<y<+® hinzugefügt. Es sei sina. =0. Durch Überlagerung dieser «drei Lösun- 
gen mit passendem Vorzeichen ist es leicht möglich, die gewünschte Lösung herzustellen, die 


als Randwerte bei «= 0 up, = sin ayfür —, =y F5 und u&,=0 für |y > — hat. 


Danach erhält man angenähert 


L 


| 
\urdyzoYe-:ex, 


wobei der Fehler dieser Näherung absolut genommen kleiner als 


l tx Se-—«ex 
T F areig —, arcte 
da T x 


ist. Als Mittelwert erhält man 


L 


| | 
y \updy a See er (1) 


mit einem Fehler, der absolut genommen kleiner als 


t | l er 1x de er 
xz 420°-+Y7° Yaz 


; = aretg Be BER Er We 
Yo x 


ist. Dieser Fehler''‘) wird mit größer werdendem Y beliebig klein. Als Mittelwert nach 
unserer Definition hat sich also bei großem Y gerade der Wert ergeben, den man vernünftiger- 
weise unter Verschwinden des Randeinflusses bei unserem leicht übersehbaren Problem er- 
warten muß. 

Die Abschätzung des Fehlers ist bei kleinem x sehr schlecht und könnte leicht durch 
Heranziehung der in Anmerk.'’) erwähnten Abschätzungsmethode bei kleinem « und |y| verbessert 
werden. In dieser Gegend interessiert der Fehler uns aber sowieso nicht. 

Bei der Diskussion des räumlichen Problems kann man etwa so voreehen, daß man ein 
spezielles up= ei"? f(x,y) wählt und für f(#,») ein ebenes Randwertproblem ansetzt, 
wiederum mit Randwerten bei #«—0, die nur im Abschnitt |y! < - von 0 verschieden sind. 


Man kann aber auch direkt mit der Poissonschen Formel für den Halbraum arbeiten und 
ganz entsprechende Schlußweisen wie in dem eben diskutierten zweidimensionalen Fall an- 
wenden. Die zeitliche Abhängigkeit hätte man durch Zufügung eines zeitabhängigen Faktors 
mitberücksichtigen können, sie macht aber bei der Mittelung einer Potentialströmung nach 
unserer Definition gar keine Schwierigkeiten. 

Im übrigen wird man bei Verwendung höherer Hilfsmittel aus der Theorie der Fourier- 
integrale direkt mit der Mittelwertsdefinition 


1 r Z 
’ +) ” > 
e— im — 4 - \ \ \ uw dtdydz 
ee be i 
v>« et vw % 
1>« u 5) F} 


arbeiten können. Die vorige Betrachtung hat uns aber einen gewissen Einblick in die 
wenig erforschte Ausbreitung von Regularitätsdefekten bei der Potentialgleichung gewährt. 


6. Die Diffusionsströmung. 


Um den zweiten Teil der Strömung, der von uns Diffusionsströmung genannt wurde, ‘ 
zu berechnen, ist die Differentialgleichung ! 
5 
j 
Orp vn R 
+ vie a ee A a a a a u h 
d t 1 m d r D ) 
im Halbraum x >0 zu lösen, wobei vp (#,y,2,t) für <=0 eine vorgeschriebene Verteilung # 
fi . 2 z x “ ä 
"po (Y,2,t) darstellen soll, während es im Unendlichen verschwinden soll. Ein entsprechendes £ 
Randwertproblem ist für wj zu lösen. £ 


11) Bei der Durchführung dieser Fehlerreehnung wurden wir dankenswerterweise durch Herrn Dipl.Ing. 
Kurt Sacher unterstützt. 
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Wir zerlegen zunächst wp nur nach der Zeit spektral; d. h. setzen 


ie 
592,4, ii, . . .-. .» - se 
es 
entsprechend 
r 2 
"m (4:2:9=J) Vm(y,2.Meittdä. 53a) 
ee) 
so daß wir für V„ die Differentialgleichung 
u 
er a) r e 
II Vp + Um Pi yAVp Et EP a a Auge A 
erhalten. Setzt man 
Hm 
IE. ee EI EERE Be 
so kommt man auf die Differentialgleichung 
a Be iz 
ı»+-7 ia a 


v +»? 


Die Randbedingung bei «=0 ist V*—= V,,; außerdem muß die Randbedingung des Ver- 
schwindens im Unendlichen von V* ebenfalls erfüllt sein. 
Nach der Greenschen Methode ergibt sich als Lös ieses Randwertproblens 
Nach der ( | Methode ergibt sich als Lösung di Randwertprobl 
+0 +0 
;. e #5 nf ö Be 
\ / / Von alye). „Ands a N: (7). 
zz —ı 
Dabei ist die Greensche Funktion @ durch Spiegelung an der Ebene #0 aus der Grund- 
lösung 49 (r) hergestellt: 
G=g(r— ger), 





r=yle — 2?’ +y—n®’+le— 7%, 
’ Swen P) -i9 
r =y(x +5) ! nn’ ri2—c6y, en 
J a; y ! (SS), 
yirn= ep" iir Ir): 
I tar P\ , mn Er 
exp (...) bedeutet et +») 
Damit ergibt sich 
oO 
e ”\ N 1+ 0 / | 0 4 Mn 
= | rs r|\eXD e 5 R . (99), 
leo Zar] 2, 7 I 23 o 


wobei 

dr: Hm | 4, 16r # a 

} = u m 

und in dieser Formel 

yet - Fra -—F 
wird. Die schließlich aus diesen Ansätzen sich ohne weiteres ergebenden Formeln für v), 
wp, tp haben ihre Eignung besonders bei der Erörterung des Einflusses eines endlichen 
Störungsbereiches bewiesen. In diesem Falle sind nur in einem sehr kleinen Gebiete der 
Ebene © =0 Randwerte, die von 0 verschieden sind, vorgeschrieben. Es haben sich dabei 
interessante Beziehungen zu den mit anderen Methoden durchgeführten Untersuchungen von 
Millionshtchikov [8] ergeben. 

Bei der zweiten Lösungsmethode verwenden wir zeitlich-räumliche Fouriertransformierte, 

d. h. es wird 


+22 +o +w 
vp(2,y,2,)=J J JS Vpla,d,u,n)etttriuytimsedidudu,. . -» . . (60) 
R 4 -% 


gesetzt. Die Fouriertransformierte der Randverteilung vn, sei Vs: 


+» +00 + 
v.(Y,2..9=J) J J Vnli,n.,a)ettinstinmzdidudn, . . . . . (60a). 
nn 2 nn x 
Man erhält dann 
FL = 5 = = = = r . 
Yn=exp & „Lem Vum +tiri+d ARE EEE: ' 


wobei das negative Vorzeichen vor der Wurzel andeuten soll, daß die Wurzel mit dem nega- 
tiven Realteil entsprechend der Randbedingung bei «=» zu nehmen ist. Oft empfiehlt es 
sich, das Argument der Exponentialfunktion in Real- und Imaginärteil zu zerspalten: 
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r Fe Um . 4 \r * 
—= exp | i—x)\ een 2 607 A: 
\ D € X] 2 y 6 Do ( ) 
Dabei ist o jetzt definiert durch 
a 16»? 4? . 
0 | n | | Be en ee 
mit 
et Te ee 


Entsprechend wird die zeitlich-räumliche Fouriertransformierte von wp: 


| " 2 
x ie) Wo. ee. 


i oO Um 
u nD=eXp Iy 


Nach der Kontinuitätsgleichung läßt sich, wie in Abschnitt 3, Gl. (11), dargelegt wurde, nun 
auch die Fouriertransformierte von u bestimmen: 


2ri0o -+ilo — u,) 0? 


0? (0 Un) + iR [u, Vo+ u, W)] ea re ee 
m . 


Up =. 2» 


Als zeitlich-räumliche Mittelwerte der Geschwindigkeitsquadrate der Diffusionsströmung 
allein erhält man nach dem Parsevalschen Satz: 


tx +0 + 


u Sn’ F - G— Um ü 
vp = 57 \ \ \ :exp| z\didu,du,, 
D T} z\ . } Do } y PaWw pP 
-n - 2 - 

ro +n + 
Sn’ ö i 6 \ 


np = TYy z\ \ \ Wp. exp (- — «) didu,du,, 








A (65). 
—ı —ı 2 
= +0 2 
—— Sn’ \ \ \ 4»? 0° 
“D TYZ) ) oo -un”’ +44 
o — 00 — 
u” Vpo + ug Wo” +2 1, My Vo r Wo» r+ Vpsi Wi) di du, du, 
Dabei ist nach Zerlegung in Real- und Imaginärteil 
Yn= Vpor +i Vosi | (66) 
IO 


Wpo= Wper+i Wosi | 


gesetzt. Eine allgemeine Beziehung zwischen den mittleren Quadraten der Geschwindigkeits- 
komponenten wie bei der Potentialströmung haben wir bei der Diffusionsströmung noch nicht 
aufgestellt. 

Wir wollen noch eine kurze Betrachtung über Zusammenwirken von Potential- und 
Diffusionsströmung anstellen. Bei der Bildung der mittleren Quadrate der gesamten Ge- 
schwindigkeitskomponenten treten auch mittlere Produkte auf, die aus Potential- und Diffu- 
sionsanteilen gemischt sind. Nehmen wir wieder zeitlich-räumliche Mittelwerte, so handelt 
es sich etwa um das mittlere Produkt 


Fo +0 + 


"Peip=r yz\ \ \vr "pdtdydz. 


m — a — 0 


Nach dem allgemeinen Parse valschen Satz (17) bekommt man dafür mit zeitlich-räumlichen 
Fouriertransformierten: 


u; ” e 


—— (2n) 

"pvp TYZ \ \ \ 

. L z 

Allgemeine Sätze über die gemischten Mittelwerte von Potential- und Diffusionsströmung 

sind vorderhand ebensowenig wie für die reine Diffusionsströmung gewonnen. Solche Be- 

ziehungen sind am ehesten zu erwarten, wenn man in der Diffusionsströmung, etwa an 

Hand der Formel (61), den Grenzübergang zu kleinem » vornimmt. Die nötigen Aussagen 

über die praktisch in Frage kommenden Bereiche von 4, u,, 4, vermögen die Messungen des 

Turbulenzspektrums und der Korrelationen R, zu liefern. Auch die Untersuchungen zur 
Dissipation der Turbulenzenergie würden heranzuziehen sein. 


VpVpdi du, du,. 
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r Über den Einfluß der Düse (oder des Auffangtrichters) 
1d auf Widerstandsmessungen im Freistrahl. 
er Von Dietrich Küchemann und Friedrich Vandrey in Göttingen. 
h (Bericht der Aerodynamischen Versuchsanstalt Göttingen e. V., Institut für theoretische Aerodynamik.) 
Es wird die räumliche Potentialströmung berechnet um Halbkörper und Quell- 
senkenkörper in einem Windkanal von kreisförmigem (Querschnitt, der zum 
5 Teil von festen Wänden (Düse, Auffangtrichter) begrenzt wird und zum anderen 
g Teil aus einem freien Strahl besteht. — Zur Erfüllung der Randbedingungen 
; wird die feste Wand mit Quellringen und die Strahlgrenze mit Wirbelringen 
N H belegt. Die Integralgleichung für diese Belegung kann durch Iteration nume- 
H risch gelöst werden. Aus der Belegung werden die zur eventuellen Korrektur der 
gemessenen Widerstandsbeiwerte benötigten Zusatzgeschwindigkeiten, sowie 
der durch den Einfluß der Düse hervorgerufene zusätzliche Widerstand 
bestimmt. Dieser Zusatzwiderstand liegt in allen praktisch wichtigen Fällen 
1g unterhalb der Meßgenauigkeit, während das Geschwindigkeitsfeld in der Um- 
in gebung des Körpers durch die Kanalwände und freien Strahlgrenzen oft 
u. merklich gestört wird. 
en I. Aufgabenstellung und Überblick über die Möglichkeiten zur Lösung der Aufgabe. 
es In der vorliegenden Arbeit soll die Frage untersucht werden, wie sich der Einfluß der Düse 
ur 


auf die Widerstands-Messungen in einem Windkanal mit freier Meßstrecke auswirkt. Der 
Kanalquerschnitt wird der Einfachheit halber kreisförmig vorausgesetzt und als Widerstands- 
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körper werden nur solche genommen, die sich durch eine Quell- und Senkenbelegung auf der 
Kanalachse darstellen lassen. Durch diese Annahmen wird die Strömung drehsymmetrisch. 
Der Widerstandskörper ruft auf der Begrenzung Normal- und Tangentialgeschwindigkeiten 
hervor, und wir haben diejenige Zusatzströmung zu bestimmen, welche auf der festen Wand 
die Normalgeschwindigkeiten und auf der Begrenzung des freien Strahls die Tangentialge- 
schwindigkeiten kompensiert. Die entsprechende Aufgabe für den reinen Freistrahl und für 
den geschlossenen Kreiskanal ist schon früher von H. Lamb u. a.') gelöst worden. Im 
Gegensatz zu diesen Arbeiten liegt bei der vorliegenden Aufgabe ein gemischtes Randwert- 
problem vor, auf das sich die bekannten Lösungsmethoden nicht ohne weiteres übertragen 
lassen. 

Es liegt nahe, die Aufgabe zunächst so zu vereinfachen, daß man einen unendlich langen 
zylindrischen Kanal behandelt, der zum Teil aus fester Wand und zum anderen aus einem 
Freistrahl besteht. Als störende Singularität nimmt man zunächst eine einzelne Quelle in 
der Kanalachse an und löst die Randwertaufgabe für verschiedene Lagen der Quelle zur Düse. 
Durch Überlagerung mehrerer solcher Lösungen lassen sich dann die Strömungen um endliche 
Widerstandskörper bestimmen. 

Es sind zwei grundsätzlich verschiedene Lösungswege zu unterscheiden. Entweder be- 
legt man die Kanalwände und Freistrahlgrenzen mit geeigneten Singularitäten und erhält aus 
den Randbedingungen eine Integralgleichung für diese Belegung, oder man sucht das Ge- 
schwindigkeitspotential oder die Geschwindigkeitskomponenten selbst durch Superposition 
geeigneter Partikularlösungen der Differentialgleichung der drehsymmetrischen Potentialströ- 
mung zu gewinnen. Beide Wege wurden von Herrn W. Flügge und nach dessen Weggang 
aus Göttingen von den Verfassern untersucht. Über diese Versuche soll zunächt berichtet 
werden, da sie die Schwierigkeiten dieser und ähnlicher Aufgaben erkennen lassen, den ein- 
geschlagenen Lösungsweg erklären und die noch offenstehenden Fragen mathematischer Art 
aufzeigen. 

1. Differentialgleichungsmethoden. Es ist das gemischte Randwertproblem der 
Potentialgleichung für drehsymmetrische Strömungen zu lösen: Auf einem Teil des Randes 
ist das Potential vorgeschrieben (Freistrahl), auf dem anderen seine Normalableitung (Kanal- 
wand). — Um eine einheitliche Lösung für das ganze Gebiet zu erhalten, verschafft man sich 
eine Funktion F (x), die auf der festen Wand O0 und auf den Freistrahlgrenzen l ist. Dann 
lautet die Randbedingung für das Potential: 

oo» 

F-D+kil-—-F) —( 

gn 
mit beliebigem, aber von 0 und © verschiedenem k. Stellt man sich nun einen Windstrom 
vor, der in periodischem Wechsel Freistrahlmeßstrecken und Kanalstrecken enthält, so kann 
man F (x) als Fourierreihe darstellen und auch für das gesuchte Potential eine Fourierreihe 
mit unbekannten Koeffizienten ansetzen. Die Lösung scheitert daran, daß eine Umordnung 
der Reihen wegen der ungleichmäßigen Konvergenz von F(x) nicht erlaubt ist. Beim Über- 
gang zu unendlicher Periodenlänge ergeben sich bei der Behandlung der Fourier-Integrale 
entsprechende Schwierigkeiten. 

Weiter wurde für Kanal- und Freistrahlstrecken ein getrennter Ansatz versucht. In 
dem Freistrahl sitzt die störende Quelle @,. Zu dem Quellpotential ist eine aus homogenen 
Lösungen gebildete Reihe zu addieren, die auf dem Strahlrand die axiale Zusatzgeschwindig- 
keit w=0 macht, und für den Kanalbereich ist eine Reihe aus Gliedern anzusetzen, für die 
auf dem Rande die radiale Zusatzgeschwindigkeit e—=0 ist. In der Ebene der Kanalmündung 
müssen die beiden Lösungen zusammengeflickt werden. Dabei erhält man ein unendliches 
Gleichungssystem für die noch freien Koeffizienten. Es zeigt sich, daß die auf diese Weise 
gewonnenen Darstellungen nicht konvergieren. Als physikalische Ursache für die Divergenz 
beider Verfahren kann vermutet werden, daß die Reihen nicht geeignet sind, den Verlauf der 
Strömung an der Kanalmündung darzustellen. 

Ein weiterer Versuch wäre eine andere Bereichsaufteilung, so daß eine schmale Zone 
beiderseits der Kanalmündung als besonderer Bereich behandelt wird, für den eine Lösung 
angesetzt werden kann, die alle Randbedingungen an der Übergangsstelle erfüllt. Der ein- 
facheren Rechnung wegen wurde die störende Quelle auch in die Mündungsebene gelegt. Das 
Zusammenflicken in den beiden Ebenen, die die drei Bereiche (Kanal — Übergang — F reistrahl) 
trennen, wurde hier nicht durch Reihenentwicklung und Koeffizientenvergleich vorgenommen, 


1) H. Lamb: On the effect of the walls of an experimental tank on the resistance of a model. ARC-Rep. and 
Mem. 1010 (1926). — C. N. H. Lock: The interference of a wind tunnel on a symmetrical body. ARC-Rep. and Mem. 
1275 (1980). — C.N. Watson: The use of series of Bessel functions in problems eonneeted with eylindrical wind 
tunnels. Proc. Roy. Soc. Bd. 130 (1930). — C.N.H. Lock und F. C. Johansen: Wind tnnnel interference on 
streamline bodies. ARC-Rep. and Mem. 1451 (1931). — H. Glauert: Interference of a wind tunnel on a symmetrical 
body. ARC-Rep. and Mem. 1544 (1933). 
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sondern durch Erfüllung der Übergangsbedingungen in einigen Punkten. Die Durchrechnung 
ergab, daß zwischen diesen Punkten die Stetigkeitsbedingungen so schlecht erfüllt waren, 
daß die Rechnung untauglich ist. Nach diesen Ergebnissen ist es unwahrscheinlich, daß 
eine besondere Behandlung der Kanalmündung die Lösung des Problems durch Reihenansatz 
erleichtert. 

2. Integralgleichungsmethoden. Als Singularitäten wählt man zweckmäßig Quellen 
auf der festen Wand und Wirbel auf den Grenzen des freien Strahls, da man auf diese Weise 
eine Integralgleichung zweiter Art erhält, die eine Lösung durch Iteration gestattet. Bei 
einer anderen Art der Belegung würde die Integralgleichung zumindest in einem der Bereiche 
von erster Art sein. Wir wollen zunächst bemerken, daß sich ein Ansatz der Belegungs- 
dichte als Fourierintegral unbrauchbar erweist, weil nicht für diesen Ansatz Randwerte zu 
erfüllen sind, sondern für ein Integral darüber. Dieser Integrationsprozeß macht die übliche, 
auf der Örthogonalität der Kreisfunktionen beruhende Koeffizientenbestimmung unmöglich. 

Eine weitere Möglichkeit wäre der Übergang zu unendlich vielen Variablen, indem man 
die stetige Quell- und Wirbelbelegung durch eine unstetige Verteilung, und zwar zweckmäßig 
durch eine stückweise lineare, und das Integral durch eine Summe ersetzt. Man erhält so 
unendlich viele Gleichungen für die Belegungsdichte in ausgewählten Punkten. Mit diesem 
Verfahren ist es nicht möglich, die Feinheiten der Strömung am Kanalrand selbst zu 
erfassen. (Sie könnten vielleicht nachträglich durch eine besondere Rechnung klar- 
gestellt werden.) Dagegen lieferte die Durchrechnung von Beispielen für die Geschwindig- 
keit in.der Kanalachse ein gutes Bild. Trotzdem sich die Ergebnisse bei einem erheblichen 
Mehraufwand an Gleichungen nicht wesentlich änderten, ist die numerische Genauigkeit bei 
diesem Verfahren nicht ausreichend, wie ein Vergleich mit den späteren Ergebnissen zeigte. 
Der Grund hierfür ist in der ungenauen Kenntnis der Belegung selbst zu suchen. Bei vor- 
sichtiger Wahl der Intervallaufteilung war der erforderliche Rechenaufwand 8 bis 10 Glei- 
chungen, die sich mit dem Rechenschieber eliminieren ließen. 

Schließlich haben wir die Möglichkeit, die Integralgleichung auf übliche Weise dureh 
Iteration zu lösen. Die Integralgleichung läßt sich in folgender Form schreiben: 


x = (feste Wand): 
1 1 1 








2 Jva + I Joy 3 q = v. 
(1). 
x =0 (Freistrahl): 
1 1 1 
garage gr 


Hierin bedeutet 4 die Quelldichte, y die Wirbeldichte der gesuchten Belegung, «, und v, die 
durch die Quelle auf dem Rande hervorgerufenen Störungsgeschwindigkeiten in axialer und 
radialer Richtung und z.B. 37 va die durch die Quellbelegung auf dem Rande induzierte 


Radialgeschwindigkeit. Setzt man eine geschätzte Lösung q‘”, y“” in die Integrale ein, so 
kann man aus der Gleichung eine neue Lösung g°”, y‘” berechnen, die unter gewissen Vor- 
aussetzungen besser ist. Die Rechnung stößt hier auf zwei Schwierigkeiten: Wechsel der 
Randbedingung an der Kanalmündung und unendliches Integrationsintervall. Der Wechsel 
der Randbedingung ist von grundsätzlicher Bedeutung. Endet nämlich an der Kanalmündung 
die Quellbelegung mit einer endlichen Dichte, so ist dort die Axialgeschwindigkeit u=x, 
und ebenso gehört zu einer endlichen Wirbeldichte am Freistrahlanfang v=». Es müssen 
also beide Belegungen an der Übergangsstelle Null werden. 

Läßt man die Ausgangsnäherung in der Umgebung dieser Stelle linear nach Null ab- 
fallen, so genügt das Iterationsergebnis nach unseren Rechnungen dieser Forderung nur dann, 
wenn der Abfall sehr steil gewählt wird. Wie sich die Belegungen in diesem Gebiet genau 
verhalten, ließ sich nicht feststellen. In einem entsprechenden ebenen Problem ergab sich 
ein Abfall proportional Y|x.. Für die numerische Rechnung genügte jedoch die erwähnte 
Annalıme eines genügend steilen Abfalls. — Das unendliche Integrationsintervall jedoch machte 
eine Bestimmung der Quellbelegung durch Iteration unmöglich. Für große negative x lautet 
die Integralgleichung (1) näherungsweise: 

1 
3; Jya 3 q U, 
oder 
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Nimmt man als Ausgangsnäherung in diesem Gebiete q"—=2v,, so geht J,, gegen 2zv, und 
damit g“’ gegen 4v,. Entsprechend ergibt sich 

WHO. 2.20.00 ‚g”’>2(n+1l)v,. 
Mit q”=v, als Ausgangsnäherung erhält man q"”>(2n+1)v,. Das Iterationsverfahren 
divergiert also. Bei der Wirbelbelegung ergibt sich diese Schwierigkeit nicht, denn J., geht 
mit «, gebildet gegen zu, und damit 

ou 


an) 


>u,: 

Das unangenehme Verhalten von 4 kann man sich plausibel machen: Betrachtet man 
das Strömungsbild eines Quellringes ?) (Bild 1), so zeigt sich, daß außerhalb des Ringes auf 
der Kanalwand nach außen gerichtete Normalgeschwindigkeiten vorhanden sind. Diese treten 
neben die Normalgeschwindigkeiten der störenden Quelle und müssen ebenfalls kompensiert 
werden. Die Quellringe heben sich also in ihrer Wirkung zumindest teilweise auf. Im 
Gegensatz hierzu unterstützen sich die Wirbelringe in ihrer Wirkung gegenseitig. Es ist zu 
vermuten, daß diese Schwierigkeit grundsätzlich an dem Ansatz von Quellen auf einer zylindri- 
schen Wand liegt. 


„0m 
Dye re... ) 























j a 
! ”n FR 
n U 
| nd a A 
| a Ss 
Düse? 
Bild 1. Stromlinienbild eines Quellringes (nach E. Trefftz, a. a. 0.?),S. 25). Bild 2. Untersuchte Düsenformen. 


Auf Grund dieser Erkenntnisse haben wir die Annahme einer zylindrischen Kanalwand 
aufgegeben und folgende Form des festen Kanals unseren Rechnungen zugrunde gelegt: Nach 
einem kurzen zylindrischen Stutzen erweitert sich der Kanal und geht in eine zur Achse 
senkrechte Wand über (Bild 2). Wir behandeln also den Ausfluß aus einer Druckkammer 
mit Stutzen. Für einen auf der senkrechten Wand liegenden Quellring ist diese Stromlinie, 
so daß sich benachbarte Quellen dort nicht beeinflussen. 

Damit sind die oben angegebenen Schwierigkeiten umgangen. — An Stelle der Parallel- 
strömung im zylindrischen Kanal tritt jetzt allerdings eine kompliziertere Grundströmung. 
Aus den Rechnungen von E. Trefftz°) geht jedoch hervor, daß sich im Stutzen alsbald eine 
Strömung herausbildet, die mit genügender Genauigkeit als Parallelströmung angesehen 
werden kann. 

II. Aufstellung der Integralgleichung. 

1. Berechnung des Kernes der Integralgleichung. Wir belegen die Oberfläche des Freistrahls 
mit Wirbelringen und die feste Wand der Düse mit Quellringen, deren Achse in der Kanalachse 
liegt. Wirbelstärke und Quellergiebigkeit sind auf einem solchen Ring konstant. Wir brauchen 
zunächst das Geschwindigkeitsfeld dieser Quell- und Wirbelringe. Es sei x die Koordinate 
längs der Kanalachse mit dem Nullpunkt in der Mündung der Düse und r die radial ge- 
richtete Koordinate, « und =» die Komponenten der induzierten Geschwindigkeit in diesen 
Richtungen. 

2) vgl. E. Trefftz: Über die Kontraktion kreisförmiger Flüssigkeitsstrahlen. Diss. Straßburg 1913 und Z. f. 
Math. und Phys. Bd. 64 (1916) S. 34 bis 61, 

3))a.2.0.2. 
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Ein Quellring im Punkte (x’; r’) mit der Ergiebigkeit 4 je Längeneinheit ıuft dann im 
Punkte (x; r) folgende Geschwindigkeiten hervor (gesamte Ergiebigkeit E=2anr'.g): 


Ip „’ . 
q 22 —r).ı . 
u— TraR er A N a ee a 
2n[le — x "+(r— r’P]yie — a? +ir+r') E(k) 
.’ [ Ielr .) 1 
2 1 h4 | K (A) MN Kl > { j m E (k) . . . (3). 
2aryle- a? Hr +’) \ (ae — a’? +lr — r’)J 





Dabei sind E und K vollständige elliptische Integrale mit dem Argument: 
2yrr 

Va’ +(r+ r' 

Für einen Wirbelring mit der Wirbelstärke y ergibt sich entsprechend mit demselben k: 


k - 


y 1 IK n I, Zr (r- r) Ein f 
2a yle— a” +(lr+r') | en‘ rg ey’ +(r — r’’J (“) LT 1“, 
y (2 — x’) | 1 2rr' - 
—- I: z = Rees +»), 
Ir r] (r Hr +r')%]| Itz- x ”+(r r ‚f Eik)| (») 


Für unser Problem brauchen wir nun die Kenntnis dieser Geschwindigkeiten nur auf 
der gegebenen Berandung. Diese setzt sich zusammen aus der zur Achse senkrechten Wand, 


b i “ie i i D 
dem Einlauf (der im Längsschnitt aus zwei Viertelkreisen vom Radius o 5 besteht) und dem 


zylindrischen Stutzen und Freistrahl vom Durchmesser D. Es ist zweckmäßig, statt der 
beiden Koordinaten =; r bzw. x’; r’ für den Aufpunkt bzw. den laufenden Punkt die mit D/2 
dimensionslos gemachte Bogenlänge o bzw. 0’ längs der Berandung einzuführen (siehe Bild 2). 
Der Nullpunkt von o wird in den Punkt gelegt, wo die senkrechte Wand in den Einlauf 
übergeht. Dieser habe die Koordinaten 

) 


X 

D N 
2r 

u 1-+o. 


Die Berandung wird dann durch folgende Beziehungen festgelegt: 
Il. Senkrechte Wand: 


(6). 


Il. Einlauf: 


III. Stutzen und Freistrahl: 


n0_ £4 
5 = St: r ee u EV 
E rn | 


D 
wobei o (7/2) 1) +ce=s gesetzt ist. 

In dem Bereich -—<o<s ist die Berandung mit Quellringen belegt, im Bereich 
s<o<-+» mit Wirbelringen. In dem ersten Bereich soll die Berandung Stromlinie werden 
und eine feste Wand darstellen, im zweiten Bereich dagegen Potentiallinie und Begrenzung 
des Freistrahls, auf der verlangt wird, daß keine zusätzliche axial gerichtete Geschwindigkeit 
aufritt. In dem ersten Bereich soll durch die Belegung die Komponente senkrecht zur Be- 
grenzung der durch einen Widerstandskörper hervorgerufenen Geschwindigkeit sowie die 
durch die Belegung selbst erzeugte normal gerichtete Geschwindigkeit kompensiert werden; 
im zweiten Bereich gilt das Entsprechende für die Axialkomponenten der Geschwindigkeiten. 
Die durch den störenden Körper hervorgerufenen Geschwindigkeiten sollen später betrachtet 
werden. Die durch die Belegung selbst induzierten Geschwindigkeiten erhalten wir durch 
Einsetzen der Koordinatenbeziehungen (6), (7) und (8) in die Gl. (2) bis (5). Eine Integration 
über 0’ von 7 =— » bis 0’ —=-+x verschafft uns dann die gesamte, von allen Quellringen 
und Wirbelringen induzierte Geschwindigkeit für jeden Aufpunkt o. Wir wollen die Kerne 
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dieser Integrale nicht alle hinschreiben. Es sind insgesamt 16 verschiedene Kerne, nämlich 
in jedem der vier Gebiete — senkrechte Wand, Einlauf, Stutzen, Freistrahl — der Einfluß 
auf sich selbst und die der übrigen drei Gebiete. Man kann sich die Kerne leicht selbst aus- 
rechnen. Wesentliche Umformungen wurden bis auf wenige Ausnahmen, die wir noch be- 
sprechen wollen, nicht gemacht, da sie die numerische Rechnung nicht vereinfachen würden. 

Besonders betrachten wollen wir zunächst den Einfluß der Quellringe der senkrechten 
Wand auf sich selbst. Wir brauchen dort die axial gerichtete Geschwindigkeit «, die sich 
nach (2) unter Benutzung von (6) bestimmt. Nach (6) ist z=.«"= -3. ce = konst., also 
= —- x’=0, Damit ist nach (2) dieser Anteil an der Normalgeschwindigkeit gleich Null in 
dem ganzen Gebiet -—=o=0. Die Quellringe der senkrechten Wand rufen an der Wand 
selbst keine normal gerichteten Geschwindigkeiten hervor, ein Ergebnis, das durch die An- 
schauung unmittelbar bestätigt wird. 

Als nächstes wollen wir das Gebiet 0o<o<7 0/2, den Einlauf, untersuchen. Hier müssen 
wir die gewünschte Normalkomponente w der Geschwindigkeit aus « und " zusammensetzen: 


oO 4 o 
WERE ee a ee ea Me 
[E) (2) 


Besondere numerische Schwierigkeiten macht die Berechnung des Einflusses der Quellringe 
dieses Gebietes auf sich selbst. Wie man vermuten konnte, lösten sich diese Schwierig: 
keiten durch Umformung des Kernes unter wesentlicher Benutzung der Tatsache, daß r (x) ein 
Kreis ist, Formel (7). Durch Einsetzen von (7) in (2) und (3) unter Anwendung von (9) er- 
halten wir von den Quellen dieses Gebietes her den Geschwindigkeitsanteil: 


To 


’ ) l j ’ 
w’ (0) = 0’): 
3,\4 ) 
u 
’ r] 

[2] [77 u. [7] [77 

2[-c0s— + cos- \(1+oli sin —|) E(k) - cos . 

0 0 | o| O 


o% ia Fi Be 
-)+{2+20 )|sın — —sın- 
[6 


( 
So, ı 
’ 





| [/ 0 2 / o'\? 


o 
COS ——+Cos + | sin ° sin .) 
o 0 0 0 





R o 
0 | COS — ur cos 
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De 


(ref —- sin )sin I > 

















(1+o sin - )] 0? (—cos i + c0s -)+(2+20 sin „tsin # j 
un rel) 


en: E(k)|:do’ 
o 6 Eh et 
o [ cos — + cos —) +| sın — —+sın = ) | 
Y 0 0 \ V 0 j] 
mit 
vor 
2 Ute 1-— sin - sylt+eli sin — 
( 
Bu -——— — = = 
RR | ir | EL RE : u 
| 0: | cos + cos 2 7 u 20—0 jun ——+sın ee 
m si; o Oo 
Der Kern läßt sich in dieser Form 2 schlecht numerisch auswerten, besonders in der Nähe 
der Stelle o—=o'. Man kann ihn jedoch umformen, so daß sich ergibt: 


10 r 








? d+oll-sin | K(k)-sin” 2E(k) 
m ()=5,\ ae) | z - . “2 7 „do (10) 
z I+o|l-sin- | e’|- eos, + cos =) + [2+20 -o sin“ „rin : 
‘ 














(k wie oben). In dieser a ist der Kern numerisch zu gebrauchen. 

Bei den übrigen Kernen sind Umformungen im allgemeinen nicht nötig. Nur in einigen 
Fällen empfiehlt es sich, von der bekannten Beziehung 

K=E+%k’-D 

Gebrauch zu machen, um bei der numerischen Auswertung nicht mit kleinen Differenzen 
nahezu gleicher großer Zahlen rechnen zu müssen. 

Schließlich möchten wir noch die Formeln für die induzierten Geschwindigkeiten im 
Stutzen und Freistrahl hinschreiben, die von den Quellen und Wirbeln in diesen Gebieten 
herrühren, da man sie vielleicht auch anderweitig verwenden kann. 
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Die Quellringe im Stutzen induzieren im Stutzen die Radialgeschwindigkeit 
RU. 2, | 
5,\16 ) Vo glK(k)-Elk]de ...... MM). 
Aue ee 


v(O)r > 


Die Wirbelringe im Freistrahl induzieren im Stutzen die Radialgeschwindigkeit 


. an er ' G a’ . 25 SE | £ | RN 5) 
vo, = \y(o Ki) 14 ng Elh|de . . . (12). 


) 5 | 
) (a —o)" +4| 
D s 


Die Quellringe im Stutzen induzieren im Freistrahl die Axialgeschwindigkeit 
5 


u(0),, Fr) )y oe gE(k)do u ° 


D To 


Die Wirbelringe im Freistrahl induzieren im Freistrahl die Axialgeschwindigkeit 








ra. 
l 
u (0), — y (0') 7 K(k) de EEE ...: A. 
D : 
Dabei ist in allen Fällen: 
») 
k= 
(oo ”+4 
Damit können wir diese Betrachtungen abschließen. Zusammenfassend können wir die 
durch die Belegung selbst auf der Berandung hervorgerufenen Geschwindigkeiten schreiben : 
+ 2 
w(0)=; \ “(o’\ K (0,0') do’ SENT FT TEEN 
sn) 
. ” = 
Dabei ist: 
(0) für -m<o<s<O 
o A 4] n To 
u (0)- cos -—+v(o)sin für O<o<5 
0 0 _ be 
nm (6) = ae . (16). 
je >» T 0 
v(0) für ><soss 
u(6) für s<o<+x 
| 1) für Bug 
u(o) = . 
ly (0) für s<o<+x. 


2. Singularitäten des Kerns. Wir 
müssen uns jetzt noch mit den be- 
sonderen Eig+ aschaften und den 
Singularitäten des Kernes K (o, 0°) 
befassen. Zur Veranschaulichung 
haben wir den Kern für ein Zahlen- 
beispiel (o=1;s=3)in Form eines 
Reliefs aufgetragen (Bild 3). 
Zunächst hat K bei o=s und 


Bild 3. Relief des Kernes der Integral- 
gleichung (Düse 1). 
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o’—s eine Sprungstelle. Außerdem hat der Kern an dieser Stelle eine unangenehme Sin- 
gularität. Betrachten wir zunächst die von den Quellringen des Stutzens in den Punkten 
des Freistrahls induzierte Axialgeschwindigkeit (Formel (13)). Im allgemeinen ist hier = 0’ 
ausgeschlossen, da s<o<+tx und zo/2<o’<s ist. Eine Ausnahme bildet jedoch die 
Stelle o=0'=s, wo k>1, E(k)>1 geht und im Nenner o — 0’ stehenbleibt*). Damit das 
Integral an dieser Stelle existiert, muß 4 (0’) = u (0°) gegen Null gehen. Diese Forderung ist 
in der Einleitung bereits ausführlich diskutiert. Genau entsprechend liegen die Tatsachen 
bei der von den Wirbelringen des Freistrahls im Stutzen induzierten Radialgeschwindigkeit 
bei o=o'=s (Formel (12)). Im Reliefbild sind diese Stellen durch die besonders hohen 
Spitzen gekennzeichnet. 

Wesentlich weniger unangenehm sind die Singularitäten, die bei den Gl. (10), (11) und (14) 
auftreten. Hier ist bei o=o’ wieder k=1 und der Kern wird wegen K logarithmisch un- 
endlich wie 

4 

Rn: IP 
Rechentechnisch läßt sich diese Singularität am einfachsten behandeln, indem man das Quadrat 
des Abstandes von der Stelle o=0’ als neue Veränderliche einführt’). Im Reliefbild werden 
diese Stellen durch die schräg bei o=0’ durch das Bild laufende Erhebung oder Senke ver- 
anschaulicht. Hätten wir noch den unendlich langen zylindrischen festen Kanal, so würde 
diese Erhebung auch bei negativen o und 0’ immer weiter laufen, nun aber bricht sie durch 
die neue Kanalform im Gebiet des Einlaufs ab und für negative o und 0’ ist der Kern über- 
haupt identisch Null. 

Schließlich erwähnen wir noch, daß in (15) über 0’ von — bis +» integriert wird. 
Schwierigkeiten treten dabei nicht auf, da sowohl « wie K genügend stark abnehmen. Auch 
hier wendet man am zweckmäßigsten das eben angegebene Verfahren an. Damit sind die 
durch die Belegung selbst hervorgerufenen Geschwindigkeiten bestimmt. 

3. Die Störungsfunktion. Aufgabe der Belegung ist es nun, die durch einen im Freistrahl 
befindlichen Körper auf der Berandung erzeugten Tangential- bzw. Normalgeschwindigkeiten m, 
zu kompensieren. w, setzt sich aus «, und », entsprechend (16) zusammen. Wir erhalten 
also die Integralgleichung 


17 
3:\ 


L 


% 
u (0) K (0, 0’) do’ — R. u()=—n,(0) .... a 
eine Gleichung, welche für jedes o erfüllt sein muß. Es handelt sich um eine Integralgleichung 
zweiter Art, die durch Iteration gelöst werden kann, da sie gestattet, aus einer Näherung u‘ 
durch Berechnung von w (u””) nach (15) bei bekanntem w, ein neues «“” +" zu bestimmen. 
Wir haben als störenden Körper zunächst einen durch eine Quelle der Ergiebigkeit E, 
in einer Parallelströmung erzeugten Halbkörper genommen. Die Quelle liege auf der Achse 
im Abstand @=a von der Düsenmündung. Quellsenkenkörper erhalten wir dann durch Über- 
lagerung einer Quelle in «=a und einer Senke in e=a-l in einer Parallelströmung. Die 
Störungsgeschwindigkeiten einer Quelle der Ergiebigkeit E, sind gegeben durch: 





gr—a 
E, eE 
un = . 
BP N 4) r\? 
[9 > 
5) +85] 
(18). 
ar 
E, - 
= 2 a = 
aD (2*= + (er) 
| a 


Durch Umreehnung auf die Veränderliche o kann daraus unter Beachtung von (16) w, (0) be- 
stimmt werden. Unsere Aufgabe ist es nun, für gegebenes w,(o) aus (17) die Belegung u (0) 
zu berechnen. Ist «(o) dann bekannt, so können wir aus den Gl. (2) bis (5) das Geschwin- 
digkeitsfeld dieser Belegung bestimmen. Wir kennen dann die bei Anwesenheit eines Wider- 
standskörpers durch die Düse und die freien Strahlgrenzen hervorgerufene Zusatzgeschwindig- 
keit in jedem Punkte des Kanals, insbesondere auf der Kanalachse. Für die letzteren ergeben 
sich übrigens einfache Gleichungen. Die Zusatzgeschwindigkeiten der Quellringe auf der Achse 
r=0 werden: 





4) Es handelt sich nieht um ein Integral mit Hauptwert, da # nur von einer Seite gegen v—s geht. 
5) Vgl. F. Vandrey: Luftf.-Forschg. Bd. 14 (1937), S. 355. 
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Für die Zusatzgeschwindigkeiten der Wirbelringe auf der Achse r—=0 erhält man: 


\r ya 2),+1ı 


nr” 
3 193 d.’ | 
| rer 


III. Ergebnisse. 

1. Die Singularitätsbelegungen des Randes. Wir haben nun einige Zahlenbeispiele durchgerech- 
net und zwar zwei verschiedene Längen des Stutzens, um dessen Einfluß festzustellen, und zwei ver- 
schiedene Lagen der störenden Quelle, um die Möglichkeit zu haben, durch Überlagerung auch einen 
Quellsenkenkörper als Widerstandskörper zu erhalten. In allen Fällen ist o— 1 gesetzt, d.h. 
die Öffnung in der senkrechten Wand hat den doppelten Durchmesser wie der Stutzen unıd 
der Freistrahl. In dem einen Beispiel (in der Rechnung mit Düse 1 bezeichnet) ist die Länge 


) 
des Stutzens (3—- 7/2)-(D/2)=1,12%-5 


und damit ce =2,4292; s=3, in dem anderen 
(in der Rechnung mit Düse 2 bezeichnet) 
\ , D 
ist die Länge des Stutzens (#5 — /2)5 
—3,4292:(D/2) und damit ce=4,4292; s=5 
(siehe Bild 2). Als Abstand der störenden 
Quelle von der Mündung der Düse wurde 
a=D/2 und a=D gewählt. — Die Haupt- 
rechenarbeit bestand in der numerischen 
Auswertung des Kerns der Integral- 
gleichung. Für das Beispiel Düse 1 ist 
der Kern in Bild 3 aufgetragen und im 
vorigen Abschnitt bereits besprochen. 
Da viele dieser Rechnungen möglicher- 
weise auch bei anderen Gelegenheiten 
gebraucht werden können (so z. B. das 
Geschwindigkeitsfeld eines Quellringes 
und das eines Wirbelringes), ist beabsich- 
tigt, die auftretenden Funktionen, die 
sich im wesentlichen aus elliptischen 
Integralen zusammensetzen, in einer 
besonderen Arbeit in Tabellenform zu 
bringen. — Die Integralgleichung (17) 
wurde durch Iteration gelöst, und zwar 
wurde im allgemeinen als Ausgangs- 
näherung für die Belegung auf der festen 
Wand das Doppelte von w, und auf 
dem Freistrahl w, selbst genommen. 
In Bild 4 sind die Belegungen u 
für 4 Zahlenbeispiele über o aufgetragen ®). 
In dem untersten Bild haben wir durch 
Überlagerung der zu einer Quelle im Ab- 
stand a=D/2 von der Düsenmündung 
(Beispiel 1) und der zu einer Senke im 
Abstand a —= D (Beispiel 2 mit entgegen- 
gesetzten Vorzeichen) gehörigen Bele- 
gungen die Belegung eingetragen, die 
von einem Quellsenkenkörper hervorge- 
6) In den Bildern 4 und 5 ist « mit E,/7 D? 


dimensionslos gemacht, jedoch die alte Bezeich- 
nung / beibehalten. 
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Bild 4. Belegungsfunktionen für 4 Zahlenbeispiele. 
Die Ausgangsnäherung ist gestrichelt eingetragen. 
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rufen wird. Man kann erkennen, in 














welcher Weise die endgültige Belegung II BEN 
von der Ausgangsnäherung (gestrichelt f\ Her vr 
eingetragen), abweicht. WR 

Die Konvergenz des Iterationsver- Zen | Een 
. . r . Y. 7 SE > ZI? 7  ; -2T 7 e 7 
fahrens ist gut. Nur in der Nähe der un RE 2 oe 
.. . 7 2; B I / — - Sf 
Übergangsstelle wurde bis zur fünften us : fi” 
Tu 7 , a) . . . WUSP | n , 
Näherung gerechnet. Es zeigte sich je- ri SL / he 

: R a a = Se 

doch keine Divergenz des Verfahrens, di I ie 





wie schon im ersten Abschnitt ange- 
kündigt wurde. Der Wert von u bei Bild 5. Belegungsfunktionen für 3 Zahlenbeispiele als Funktion 
g=e wurde nicht berechnet, sondern 4 Abıtunler ya der Mu, Damm Verbin I a nelseuner 
immer zu Null angenommen. Für große 

positive und negative o fällt die endgültige Belegung praktisch mit der Ausgangsnäherung 
zusammen. Um die Belegungen besser miteinander vergleichen zu können, sind die ersten 
drei Beispiele gemeinsam in Bild 5 über den Abstand von der störenden Quelle aufgetragen 
und zum Vergleich dazu die Belegung, die die Quelle auf einem unendlich langen zylindrischen 
Freistrahl hervorruft. 

Bemerkenswert ist besonders, daß sich sämtliche Belegungen rechts von der Quelle 
nicht merklich von der Freistrahlbelegung unterscheiden. Bis dahin reicht der Einfluß der 
Düse augenscheinlich nicht. Auch eine Änderung der Stutzenlängen macht in der Nähe der 
störenden Quelle nicht viel aus. Auffallend ist jedoch die starke Zunahme der Quellstärke 
in der Nähe der Düsenmündung, wenn die störende Quelle näher an die Düse rückt. 

2. Die Zusatzströmung. Um eine anschauliche Vorstellung von der durch die Belegung 
hervorgerufenen Zusatzströmung zu bekommen, ist im Bild 6 das Stromlinienbild der Zusatz- 
strömung für ein Beispiel (Düse 1;2 a / D=2) gezeichnet. Es ist an diesem Bild bemerkenswert, daß 
die Zusatzströmung in unmittelbarer Nähe der störenden Quelle aufder Achse einen Staupunkt be- 
sitzt. Bild 6 zeigt außerdem das Stromlinienbild für Zusatzströmung und Quellströmung zusammen. 
Man sieht, daß die Bedingung, daß der Freistrahlrand Potentiallinie ist, recht gut erfüllt ist. An 
der Übergangsstelle «=s oder in deren Nähe auf der Begrenzung des freien Strahls ergibt 
sich ein Staupunkt. Damit die Strömung physikalisch sinnvoll wird, ist zu erwarten. daß 
dieser genau bei o=s liegt, denn sonst ergibt sich zwischen dem Staupunkt und der Über- 
gangsstelle eine nach innen gerichtete Strömung, welche die feste Wand an der Übergangs- 
stelle mit unendlich hoher Geschwindigkeit umfließt. Wie die Verhältnisse an dieser Stelle 
wirklich liegen, und welche Lösung die gestellte Potentialaufgabe in diesem Gebiet genau 
hat, konnte durch numerische Rechnung nicht festgestellt werden. Es bedarf dazu einer ge- 
naueren Untersuchung der Übergangsstelle. In unserem Zahlenbeispiel Jag der Staupunkt 
nicht genau bei o=s. Da die Rechnung jedoch in diesem Gebiet nicht so genau durch- 
geführt werden konnte, um derartige Abweichungen zu garantieren, haben wir in Bild 6 das 
Stromlinienbild so abgeändert, daß sich eine physikalische sinnvolle Strömung mit dem Stau- 
punkt an der Übergangsstelle ergibt. Weiter im Stutzen ist der gewünschte Strömungs- 
verlauf wieder gut erfüllt. Es ergibt sich im Stutzen praktisch eine negativ gerichtete 
Parallelströmung von kleinem Betrage. Schließlich haben wir in Bild 7 das Stromlinienbild 
der Gesamtströmung um einen Halbkörper vor uns. Wir haben die Geschwindigkeit U der 
Parallelströmung so gewählt, daß sich für den Durchmesser d des Halbkörpers nach 





a a EEE ee _ 
al 
der Wert d=0.707 D ergibt. Wir wählten einen solchen dieken Körper, damit die an sich 
kleinen Änderungen im Stromlinienbild gut herauskommen. Da wir diese Strömung durch 
einfache Überlagerung von Zusatzströmung, Quellströmung und Parallelströmung erhalten, 
können wir durch Umkehrung der Parallelströmung außer dem Ausfluß aus einer Düse gegen 
einen Halbkörper (oberes Bild) auch die Strömung um einen Halbkörper im Freistrahl mit 
Auffangtrichter konstruieren (unteres Bild). Wie die Strömung weiter hinten im Auffang- 
trichter verläuft, können wir der besonderen Form der festen Wand wegen jedoch nicht 
feststellen. — Beide Stromlinienbilder unterscheiden sich nicht viel von der ungestörten 
Strömung in unendlich ausgedehnter Flüssigkeit (gestrichelt eingetragen). _ Insbesondere 
machen sich die Besonderheiten der Zusatzströmung in der Umgebung der Übergangsstelle 
(s=s) kaum bemerkbar. Der Staupunkt am Halbkörper ist etwas nach hinten verschoben. 
3. Der Einfluß der Zusatzströmung auf Widerstandsmessungen. Damit kommen wir auf 
die Wirkungen, welche eine durch Düse und Freistrahl abgeänderte Strömung auf einen 
Widerstandskörper ausübt und auf die Frage, in welcher Weise Widerstandsmessungen in 
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Bild 6. Stromlinienbilder. Oben: Zusatzströmung. Unten: Zusatzströmung + Quellströmung (Düse 1;22/D—2). Zwischen je 
zwei benachbarten ausgezogenen Stromlinien strömt 1/50 der Quellflüssigkeit hindurch, zwischen je zwei gestriehelten 1/400. 
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Bild 7. Stromlinienbilder von der Umströmung eines Halbkörpers. Oben: Ausfluß aus einer Düse. Unten: Vor dem 
Auffangtrichter. Gestrichelt eingetragen ist das Bild in unendlich ausgedehnter Flüssigkeit (Düse 1; ?a/D —?). 
einem Windkanal mit freier Meßstrecke eventuell korrigiert werden müssen. Bei einer solchen 
Messung wird die Kraft W auf den Körper gewogen, die Anströmungsgeschwindigkeit 7, bei 
der diese Kraft auftritt, gemessen und daraus gebildet 
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wobei o die Dichte der Flüssigkeit und F eine Bezugsfläche, hier die größte Querschnitts- 
fläche des Körpers, ist. 

Zunächst einmal muß bemerkt werden, daß der Körper im Kanal eine andere Druck- 
verteilung besitzt, als in unendlich ausgedehnter Strömung. Man müßte daher diejenige 
Körperform bestimmen, welche in unendlich ausgedehnter Strömung dieselbe Druckverteilung 
hat wie der gemessene Körper im Kanal, und hätte dann den Körper, dessen Widerstand 
man gemessen hat. Diese Anderung der Druckverteilung wird jedoch im wesentlichen von 
den freien Strahlgrenzen herrühren und nur in 














Bi ganz geringem Maße von der Düse. 

FA u./0/° Man sieht dies aus der in Bild S gezeigten 
F \ v Id Geschwindigkeitsverteilung, aus der hervorgeht, 
7 Pi daß gerade im hinteren Teil des Körpers, der 
: We | für den Widerstand maßgebend ist, ein Einfluß 
X) Zu der Düse nicht mehr zu erkennen ist. Im Rahmen 
Je ar unserer Aufgabenstellung kann dieser Einfluß 

Ri / f außer acht gelassen werden. 
Si J / Im Zusammenhang mit der Änderung der 
/ / ‘  Körperform ist noch darauf hinzuweisen, daß 
= 2 n 7 ; auch die Bezugsfläche F durch Düse und freie 
—.32 Strahlgrenzen bei endlichen Körpern (nicht beim 
Diser £ -2 Halbkörper) geändert wird. Diese Änderung 
a Er können wir nach den Ergebnissen unserer 
2 Rechnung abschätzen: Es ergibt sich eine merk- 
_——— 2 v2 2 - ‘  liehe Änderung von F von etwa 1 v. H. und mehr 
EM Freistrahl | bei Quellsenkenkörpern, deren Vorderkante etwa 


einen Kanalradius von der Düse entfernt ist, 
| erst bei Dickenverhältnissen d/ D über 0,5. Auch 
| dieser Einfluß ist fast ausschließlich durch den 
Freistrahl und nicht durch die Düse bedingt. 
Bild 8. Zusatzgeschwindigkeiten in der Achse als Einen Einfluß wird bei allen diesen Über- 
Funktion des Abstandes von der Quelle. legungen auch die Deformation des freien 
Strahlshaben, was hier nicht berücksichtigt wurde. 
Außer der oben bereits erwähnten Änderung des Widerstands infolge anderer Geschwindig- 
keitsverteilung und damit anderer Ablösebedingungen kann auch schon in der Potentialströmung ein 
zusätzlicher Widerstand oder Vortrieb entstehen. Hierzu könnten wir grundsätzlich die be- 
kannte Tatsache °) heranziehen, daß eine Quelle der Ergiebigkeit E, in einem beliebigen dureh 
Singularitäten hervorgerufenen Geschwindigkeitsfeld die Kraft — « E, u* erfährt, wobei o die 
Dichte der Flüssigkeit und «* die Geschwindigkeit am Ort der Quelle bedeutet. In unseren 
Falle hat aus Symmetriegründen «* die Richtung der Kanalachse, und wir erhalten auch eine 
Kraft in Richtung der Kanalachse, den Vortrieb 


A2058 


> -008 — 


en rn: 25 


weil nach der Kontinuitätsbedingung E,= F- U ist. Da nun in unmittelbarer Nähe der Quelle 
ein Staupunkt der Zusatzströmung liegt, wird «* und damit 4 W immer verhältnismäßig 
klein sein und liegt innerhalb der Rechenunsicherheit. Genauer können wir 4 W für Quell- 
senkenkörper auch mit Hilfe einer Impulsbetrachtung bestimmen. Wir benutzen dazu die 
Tatsache, daß ein Teil der Quellflüssigkeit in den festen Kanal hineinströmt (vgl. Bild 6 
unten), wodurch dort eine der Anströmung entgegengesetzte Parallelströmung hervorgerufen 
wird. Diese habe in der Düse die mittlere Geschwindigkeit @ (an anderen Stellen des festen 
Kanals eine dem veränderten Querschnitt entsprechende). Nach der Kontinuitätsbedingung 
ist die in die Düse einströmende Quellflüssigkeit 
] 1 -cosg): I D:?.ü 
wobei p der Winkel zwischen der auf den Staupunkt am Düsenende zulaufenden Stromlinie 
und der Achse am Ort der Quelle ist (siehe Bild 6). Daraus erhalten wir unter Benutzung 
von (21): 

„H 


l I\? 
jryi cos) 7) Be ED 


?), Vgl. A. Betz: Singularitätenverfahren zur Ermittlung der Kräfte und Momente auf Körper ın Potential- 


strömnng. Ing.-Areh. Bd. 3 (1932), S. 454 bis 462. 
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Dabei ist x von der Lage des Widerstandskörpers zur Düse abhängig, also von 2a/D. In 
dem gerechneten Beispiel mit Düse 1:2a/d=2 ist 9= 14", so daß «/U= 0,015 (d/D)” wird. 
Legen wir nun zwei Kontrollflächen senkrecht zur Achse durch den Stutzen (Geschwindigkeit 
u,=U-—ü, Druck p,, Querschnitt F,) und weit hinter dem Quellsenkenkörper durch den 
Freistrahl (u,= U; p,; F,), so liefert der Impulssatz eine zusätzliche Kraft auf den Körper 


AW=F,(u-+p, - Pp)- F,ou}. 


Diese Beziehung könnten wir mit Hilfe der Bernoullischen Gleichung und der Kontinui- 
tätsgleichung umformen in 


or D? 


3 aa)? 


[8 
= ' W.. 
AW=—F,(w,-u)- 
Nun ist u, — u, aber die oben eingeführte Geschwindigkeit @ und nach (23) 


. Bi (1—cosp) (d\’ 
w-gu.r Tel). :--.22022: 0 
In dem Zahlenbeispiel mit 9= 14° wird (1— cos g) ?/4 = 0,0002, so daß 4 W wiederum ver- 
nachlässigbar klein ist. Befindet sich der Körper näher an der Düse, wird g größer, was 
man jedoch in jedem Falle leicht abschätzen kann. 

Schließlich besteht die Möglichkeit, daß die Geschwindigkeit U korrigiert werden müßte, 
In der Praxis wird die Anströmgeschwindigkeit meist in der Weise bestimmt, daß der stati- 
sche Druck p; in der Düsenvorkammer gemessen und daraus durch Umrechnung mittels des 
Düsenfaktors 1/[1 — (Fa/F7)?] die Geschwindigkeit U erhalten wird aus: 


eo JR Pk Po von 
I s A „ee, 2er Se Dr . . R . . . . . A, 
1 (7) 


wobei p, der Druck im Außenraum, Fj die Querschnittfläche der Düsenvorkammer und 
F4=nD?°j4 der Düsenaustrittsquerschnitt ist. Nach den Voraussetzungen, unter denen diese 
Beziehung abgeleitet ist, ist 7/ die Geschwindigkeit am Strahlrand, was sich auch durch die 
Zusatzströmung nicht ändert, da diese an der Düsenmündung auf dem Rande einen Staupunkt 
hat und so bestimmt ist, daß auf dem Strahlrand keine zusätzliche Axialgeschwindigkeit 
auftritt. Ist diese Geschwindigkeit auf dem Strahlrand maßgebend für den gemessenen 
Widerstand, so ist gegen die Messung der Anströmgeschwindigkeit mittels des Düsenfaktors 
nichts einzuwenden. — Es besteht jedoch die Möglichkeit, daß nicht die Verhältnisse auf 
dem Strahlrand, sondern die Mittelwerte von Druck und Geschwindigkeit in der Düsenmündung 
maßgebend sind, eine Geschwindigkeit, die, wie oben gezeigt wurde, um den Betrag ü nach 
(23) vermindert ist. Dies würde jedoch nach dem angegebenen Wert für 9 den Staudruck 
erst bei einem Durchmesserverhältnis d/D = 0.58 um 1 v. H. ändern. Wir schließen daraus, 
daß auch ohne Klärung der Frage, welche Geschwindigkeit maßgebend ist, gegen die übliche 
Methode, den ungestörten Staudruck zu bestimmen, nach unseren Ergebnissen im allgemeinen 
nichts einzuwenden ist. j 

Weiterhin ist es möglich, daß die Anströmgeschwindigkeit mittels Staurohr an einer 
Stelle gemessen wird, an welcher durch die Quell- und Wirbelbelegung noch eine Zusatz- 
geschwindigkeit «’ auftritt. Anstatt die ungestörte Anströmgeschwindigkeit U zu messen, mißt 
man dann U+ uw’ und muß nachträglich korrigieren. 

Den Wert von «’ können wir aus unseren Belegungen berechnen. «’ setzt sich additiv 
zusammen aus den entsprechenden Werten der Zusatzströmung und der Quellströmung. Uns 
interessiert hier zunächst der Anteil, welchen die Zusatzströmung liefert, da man sich den 
der Quellströmung leicht berechnen kann. In Bild 8 ist daher die Zusatzgeschwindigkeit in 
der Kanalachse aufgetragen. Man erkennt wieder, daß die Länge des Stutzens keinen großen 
Einfluß hat, wohl aber die Lage der störenden Quelle. Rechts von der Quelle unterscheiden 
sich jedoch wie früher die Werte kaum von denen, die man in einem unendlich langen Frei- 
strahl erhalten würde. Der Einfluß der Düse macht sich deutlich in einem Absinken der 
Werte bemerkbar. In Bild 9 ist die Zusatzgeschwindigkeit in der Strahlachse noch einmal 
für verschiedene Abstände der störenden Quelle von der Düse aufgetragen. Die Kurven für 
2a/D=0,8; 1,2; 1,4; 1,6; 1,5 sind dabei linear interpoliert, was sich wohl verantworten läßt, 
da die lineare Interpolation rechts von der Quelle exakt richtig ist. Um zu sehen, in welcher 
Weise die Zusatzgeschwindigkeiten im ganzen Windstrom von den Werten auf der Achse 
abweichen, wurde Bild 10 gezeichnet, allerdings nur für das Zahlenbeispiel mit Düse 1, 
2a/D==2. Für andere Lagen der störenden Quelle werden die Abweichungen ähnlich sein. 
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a | & T are T sr Bild 9 (links). Zusatzgeschwindigkeiten in der 

4.0) | Achse als Funktion des Abstandes von der Düse 

5 | für verschiedene Abstände der Quelle von der Düse. 

| Die dünn ausgezogenen Kurven sind interpoliert 
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Bild 10 (unten). Zusatzgeschwindigkeiten als Funk- 
tionen des Abstandes von der Düse für verschiedene 
2r|D (Düse 1; 2a4/D=2) 
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Zusatzgeschwindigkeiten in der 


Bild 12. 


Zusatzgeschwindigkeiten in der Achse als Funktion des 


Bild 11. 

Abstandes von der Düse für Quellsenkenkörper verschiedener Achse. a Zusatzströmung. 5b Quell- 

Sehlankheit. Die Quelle liegt bei 2a/D=| Der Abstand 21/D strömung. ce Zusatzströmung + Quell- 
der Senke von der Quelle ist variiert (Düse 1). strömung (Düse 1: 2a/D —2). 
Schließlich wurden durch Überlagerung von einer Quelle und einer Senke der gleichen 

Ergiebigkeit E, im Abstand I die Zusatzgeschwindigkeiten « in der Kanalachse für Quell- 


senkenkörper verschiedener Länge gewonnen (Bild 11). 


Dabei wurde die Dicke d des Quellsenkenkörpers nach Lamb °) aus 


0.) Formel (2). 
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urn _ + El e 
a de RE EEE N (27) 


bestimmt und damit die Geschwindigkeit E,/z D?, mit welcher dimensionslos gemacht wurde, zu 


bo ® v(5)] 13 d 2 
aD 4 \D +7) 
Es ergeben sich, wie Bild 11 zeigt, in der Nähe der Düse bemerkenswert kleine Werte für 
die Zusatzgeschwindigkeit. Setzen wir beispielsweise d/D=0,5 und diI=1, so ergibt sich 
die Zusatzgeschwindigkeit in der Mündung der Düse für den Quellsenkenkörper zu u = 0,0025 T, 
während diese für einen Halbkörper derselben Dicke 0,01 U beträgt. Im unendlich langen 
zylindrischen Freistrahl dagegen ergeben sich aus Symmetriegründen die entsprechenden 
Verhältnisse, wie sie in Bild 11 (für Bild 9 gilt dasselbe) hinter dem Körper vorliegen (der 
Einfluß der Düse reicht nicht bis dahin), so daß an der entsprechenden Stelle der Achse bei 
einem Quellsenkenkörper # = 0,0056 U ist. Damit haben wir das Geschwindigkeitsfeld der 
@Quell- und Wirbelbelegung bestimmt. 

Zu der Geschwindigkeit «’ liefert nun auch die Quellströmung einen Anteil, wie oben 
bemerkt wurde. Dieser übersteigt den der Zusatzströmung im allgemeinen beträchtlich. 
Um einen Überblick über die Größenordnung zu erhalten, haben wir in Bild 12 die Zusatz- 
geschwindigkeiten in der Kanalachse für die Zusatzströmung, die Quellströmung und deren 
Summe für ein Beispiel aufgetragen. Als Wesentlichstes erkennen wir, daß sich trotz Ab- 
minderung durch die Zusatzströmung vor dem Körper immer noch ein Stau ausbildet. 

Vergleiche mit Messungen konnten leider nicht vorgenommen werden, da Messungen in 
dieser Hinsicht kaum vorhanden sind. Die einzigen Messungen, bei denen die Anström- 
geschwindigkeit mittels Staurohr an verschiedenen Stellen des Kanals gemessen wurde’), 
sind mit kreisförmigen, senkrecht angeblasenen Platten vorgenommen, bei denen die Wider- 
standsverhältnisse durch Bildung des Totwassers gänzlich andere als die hier angenommen 
sind. Es ergibt sich nämlich ein mit größeren Durchmessern der Scheiben abnehmendes c,„ 
im Gegensatz zu den Beobachtungen an Luftschiffkörpern '*), bei denen c„ mit wachsender 
Dicke des Körpers zunimmt. Es ist jedoch nicht sicher, wieweit diese Erscheinungen durch 
Herstellungsungenauigkeiten hervorgerufen wurden, so daß ein Vergleich mit unseren Ergeb- 
nissen nicht vorgenommen wurde. 

Herrn Dr.-Ing. habil. W. Flügge sind wir für die Überlassung seiner Untersuchungen, 
sowie Frau Dr.-Ing. I. Flügge-Lotz für wertvolle Hinweise bei der Weiterführung der 
Arbeit zu besonderem Dank verpflichtet. Bei der Leitung der umfangreichen Rechenarbeiten 
hat uns Fräulein stud. math. I. Multhopp in dankenswerter Weise unterstützt. 


IV. Zusammenfassung. 

In der vorliegenden Arbeit ist der Einfluß der Düse oder des Auffangtrichters auf 
Widerstandsmessungen in einem Windkanal mit- freier Meßstrecke untersucht. In einem 
kurzen Bericht werden die versuchten Lösungsmöglichkeiten erörtert, und es wird auf 
die noch bestehenden mathematischen Probleme hingewiesen. Um mit den bekannten 
mathematischen Hilfsmitteln zu einer Lösung zu kommen, werden als Düse die Wände 
einer Druckkammer mit Stutzen gewählt (Bild 2). Zur Erfüllung der Randbedingungen, 
welche verlangen, daß die durch den Widerstandskörper auf der Kanalwand hervorgerufenen 
Normalgeschwindigkeiten und die auf der freien Strahlgrenze hervorgerufenen Tangential- 
geschwindigkeiten verschwinden müssen, wird die feste Wand mit Quellen und die Strahl- 
grenze mit Wirbeln belegt. Die Randbedingungen liefern für die Belegung eine Integral- 
gleichung zweiter Art, die durch Iteration numerisch gelöst werden kann. Aus der Belegung 
wurde dann das Strömungsfeld um Halbkörper und Quellsenkenkörper beliebiger Dicke in ver- 
schiedenen Lagen zur Düse berechnet, sowie die xraft, welche ein solcher Körper in der 
Zusatzströmung erfährt. Es zeigte sich, daß dieser zusätzliche Widerstand stets sehr klein 
ist und im allgemeinen unterhalb der Meßgenauigkeit liegt, während das Geschwindigkeitsfeld 
in der Umgebung des Körpers durch die Anwesenheit der Kanalwände und freien Strahl- 
grenzen oft merklich gestört wird. 205 


Abgeschlossen März 1940. 
9) siehe C. Wieselsberger: Widerstandsmessungen im freien Luftstrahl und Kanal. Z. Flugtechn. Bd. 6 (1915), 


S. 127, oder Mitt. d. AVA Göttingen 1. Folge, S. 54; vgl. auch Handb. d. Exp.-Phys. Bd. IV, S. 249, 


10), C. Wieselsberger: Ähnliehkeitsuntersuchungen an Ballonmodellen und Versuche über den Einfluß der 
Oberflächenbeschaffenheit. Mitt. d. AVA, 1. Folge, S. 53. 
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Über die Minimalsätze der Elastizitätstheorie. 


Von Constantin Weber in Dresden. 


Im folgenden werden die Minimalsätze der Elastizitätstheorie auf möglichst 
breiter Grundlage abgeleitet und dann nach Einführung gewisser Ein- 
schränkungen ihre Anwendung zur Bestimmung elastischer Konstanten 
dargelegt. 


Il. Der Extremwert der potentiellen Energie. 

Ein belasteter Körper sei im Gleichgewichte. Die Koordinaten eines Punktes des un- 
belasteten Körpers seien «x, y, 2, desselben Punktes des belasteten Körpers =“ + u, „=y-+r, 
=z+mw. Im Teilchen d£dndZ£ sei infolge der Formänderung die Zerrungsenergie ad£dydS 
aufgespeichert, die von der Endform des Teilchens, nicht aber von der Art der allmählichen 
Entstehung dieser Endform abhängt. Das Gesetz zwischen der spezifischen Zerrungsenergie « 
und den Zerrungen (den Größen, die die Formänderung kennzeichnen) sei noch ganz beliebig; 
es kann auch von Ort zu Ort verschieden sein. 

Nun wird dem Körper eine kleine zusätzliche Formänderung erteilt, so daß der Punkt £&, n, Ö 
in den Punkt E+ödu, 7 +0v, © +6ömw übergeht. Die Verschiebungen dw, öv, dw und ihre Ab- 
leitungen können beliebig klein gemacht werden. Im Teilchen d&d»d£ ist nunmehr die 
Zerrungsenergie (a+da)d£dnd£ enthalten. da entsteht dadurch, daß die Seitenflächen des 
Teilchens d£dnd£ sich in verschieden starkem Maße infolge du, dr, ömw verschieben. da 
hängt also von den Ableitungen der öu, dv, dw ab. Verschieben wir die Fläche &= const 


RER A B a i . h 
um dw, die Fläche E+dE2=const um du —+ 48, während alle weiteren Verschiebungen 


“ 


Pr ' j En ; „döu,.,., f 
des Teilchens null sind, so wird die Formänderungsarbeit o:dnd/- JE dS in Zerrungsenergie 
2 , ö du ’ ä a 2 z 
umgesetzt. Folglich wird da = o: 2. Verschiebt man die Fläche n=const um du, die 
1 08 
z, ou } i n N - ; -3r 
Fläche „+dny=const um du+ = d], so wird die zugeführte Formänderungsarbeit r:, d&d{ 
7 
VR)) VR>) 
gu an i e vou an m 
‚, 4»; in Zerrungsenergie umgesetzt. Man erhält dann da=tr;, En; Treten alle mög- 
ur 7] 
lichen du, dv, öm auf, so wird 
s du + r oodu oodr 1 
da = 0- — +4 T- = e)+: ; 
ae 7- "(09 08 u 


(Punkte für Glieder, die durch zyklische Vertauschung entstehen.) 
Der Zuwachs der Zerrungsenergie A des gesamten Körpers wird 


ws rt a8 15 ER 
A=\\\dadzdyd£ \\le ge + tr + E)+|dsanae. 


p- r 

Durch partielle Integration wird dieses Integral umgeformt. Beachtet man die Gleich- 
00: OT, OT; Ä 

-% 3: +X=t0( und zyklisch, und für das Öber- 


gewichtsbedingungen im Inneren: + \ = 
= Os ( 


flächenteilchen do: 

=do=0;docosa+r;,docosß+Tr;.docosy und zyklisch, 
worin X, Y, Z die Massenkräfte je Raumeinheit und Z, A, Z die Oberflächenkräfte je Flächen- 
einheit und a, , y die Winkel der Flächennormalen mit den Koordinatenachsen sind, so erhält man: 


SA=\|\I(Xöu+ Yör+Zöm)dedndi+ \(Zödu+ Hör +Zöm)do  ° 
4 0 


Bringt man alle Glieder auf die linke Seite, so stellt der gesamte Ausdruck links die 
Variation der potentiellen Energie dar. Die potentielle Energie wird folglich zum 
Extremum. Um die Frage zu entscheiden, ob ein stabiles, indifferentes oder labiles Gleich- 
gewicht vorliegt, muß die zweite Variation gebildet werden. Ist diese positiv, die potentielle 
Energie selbst also ein Minimum, so könnte der variierte Zustand nur durch Hinzuführung 
von Energie erreicht werden. Das Gleichgewicht ist also ein stabiles. 

In dieser weitesten Form dient der Satz vom Extremum der potentiellen Energie als 
Ausgangspunkt zur Behandlung von Belastungsfällen mit großen Verschiebungen und großen 
Zerrungen und zur Untersuchung von Stabilitätsfragen. 
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II. Das Minimum der potentiellen Energie für kleine Zerrungen und Drehungen. 

Die Zerrungsenergie im Teilchen d£dnyd£ ist adödnd£. Bezogen auf das Teilchen 
dxdydz, erhält man eine andere Zerrungsenergie infolge der Änderung des Rauminhaltes 
des Teilchens. Der Rauminhalt 1 des unbelasteten Körpers geht über in den Rauminhalt 
(1-+.p) des belasteten Körpers. Die Zerrungsenergie, bezogen auf das unbelastete Teilchen, 
wird darum a(l+ey)da@dyda. 

Die Zerrungsenergie a(l1+ 7) und die Raumdehnung &;- hängen von den Zerrungen des 
Teilchens ab, hierdurch also von den Ableitungen nach «, y, z der Verschiebungen «, ®, w. 
Von diesen neun Ableitungen bestimmen drei (für unsere Untersuchung nicht wichtige) Aus- 
drücke die Drehung des Teilchens, während die sechs Ausdrücke e,, und zyklisch, &, „ und 
zyklisch die Zerrungen geben. Diese Zerrungen sind '): 


P) »] 


N 1 liäe (a0 3 
eu ou ot en | . 

as=3.+t3 IR 3 +5 A E 3.) | und zyklisch, 
ou dv du du  dvdr., dmodm] 


f 
| . 
Ey =tyrn no | und zyklisch. 


= en“ . 4 E; | 

0y 0x dx dy drdy Owody) 

Die Zerrungsenergie a und die Raumdehnung e7- sind Funktionen dieser sechs Zerrungen ; 

die erstere hängt von den Werkstoffeigenschaften des Körpers ab. Bei homogenen Körpern 

müßte diese Funktion für alle Stellen «, y, z dieselbe sein; für isotrope Körper müßte a eine 
Funktion der drei Invarianten J,, J,. ./, der Zerrungen sein: 


J, =exxteyyt &: 


. J,=4eyxtyyt: — Ey" 
ey Exyı Ex: 
| 
J,= Eyx, 2Eyy. Eyz 
Ezxy Eryı Erz 


Wir wollen hier solche Einschränkungen nicht vornehmen. 

Wir führen hingegen die Einschränkung ein, daß alle Ableitungen klein gegenüber I 
sind. Dann kann erstens ej- im Ausdrucke (1-+ ep) vernachlässigt werden. Weiter gehen die 
Zerrungen über: 


N 0 5 1 kli l 
Erx ın Ex und ZVKIISCH, 
x x or . 
r du ov i 
Ex ın Y = + und zyvklisch. 
xu v‚x4 Ri) Yy Ö x A 


Dies sind die üblichen Dehnungen und Schiebungen der Elastizitätslehre für kleine 
Zerrungen und Drehungen. 
Es kommt natürlich darauf an, daß die vernachlässigten Glieder für das Integral 
A =| \ \ ad£dnd£ ohne Bedeutung sind. Das ist wichtig, weil es vorkommen kann, daß 
rs \2 
0 ist, folglich in #,. die Glieder 35) und 3 (*) nur für das Integral bedeutungs- 
los werden. 
Die Zerrungsenergie a (ex, ..&xy, .JdSdnd£ geht über in alex... Yay...)dedyda. 
Wir wollen nunmehr da bestimmen, das infolge der Variationen du, dr, dm entsteht. 
d(u+du) 
02 


es Er OU „. “ i j 
Hierfür setzen wir für e.=s-, die variierte Ableitung usw. ein. 


da 
Dann wird 


ou du, or 
a(e,y)=a (5 +3.) 


und. 
o(u+du) do (u+du) dle+drv) 
y — — - 
a(e,y)+da a Ir ER dy 2; Ei} 
Hieraus: 
, 2a don da sa A 
Eye, da" Ayın\öy dr 
Durch Vergleich mit Gl. (1) erhält man 
da(e,y da(e,y) : ; 
en ) =0,, N on zu weh. . . sr 


1) Love: Lehrbuch der Elastizität, deutsch von A. Timpe. Leipzig und Berlin 1907, S. 59, S. 71. 
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Der Satz vom Extremum der potentiellen Energie nach Gl. (2) gibt: 
\ilda (e,yJ)dedydz \ \(Xdu+Yor+Zdnm)dedydz | 
j 1) 


(4). 
\(Zödöu+Hör+Zon)do=0 | 


Die Massenkräfte X, Y, Z und die Oberflächenkräfte £, H, Z mögen unabhängig von 
den Verschiebungen, Zerrungen und Drehungen sein. Dieses läßt sich physikalisch begründen: 
Da die Raumdehnung so klein ist, daß sie gegen der eins vernachlässigt werden kann, werden 
sich die Massenkräfte nicht ändern; die Oberflächenkräfte, die von der Größe und Lage der 
Oberfläche abhängen, werden ebenfalls praktisch konstant bleiben. 

Für die Variation der Verschiebungen «, v, w machen wir nunmehr folgende Vorschrift: 
An jeder Stelle der Oberfläche ist entweder die Verschiebung oder die Oberflächenkraft ge- 
geben. Die Oberflächenteilchen do teilen wir ein in die Teilchen do,, für die die drei Ver- 
schiebungskomponenten «, v, w und in die Teilchen dog, für die die drei Komponenten der 
Öberflächenkraft Z, H und Z gegeben sind. Allgemein könnte man in jedem Öberflächen- 
teilchen do ein rechtwinkliges Koordinatensystem anbringen und für einen Teil der Richtungen 
die Verschiebungskomponenten, für den restlichen Teil die Kraftkomponenten angeben. Der 
Einfachheit der Darstellung wegen wollen wir uns auf den obigen Fall beschränken, bei dem 
die Einteilung in die Teilchen do, und dog genügt. 

Bei der Variation der Verschiebungen lassen wir nur solche zusätzliche Verschiebungen 
du, öv, dw zu, die an den Oberflächenstellen do, gleich null sind, da dort die Verschiebungen 
gegeben sind. 

Die potentielle Energie //, die zum Extremum wird, läßt sich dann wie folgt schreiben: 


HI = \a(s,y)dedydz - \(Xu+Yr+Zm)dedydz | 
v v 


-\ (Zu + Hv + Zm) dog = Extremum | 
ON 
Die zweite Variation von // besteht nur aus der zweiten Variation der Zerrungsenergie 
A=\a(ey)dedydz:; diese Variation hängt von der spez. Zerrungsenergie a (e, y), genauer 
r 
von deren zweiten Variation ab. 
du odu , ddr 
: di; : 


Die zweite Variation von a (e,y) wird, wenn wir de, 


Öy,, und 
Ö x x, oy Ri) x ixUY 


zyklisch schreiben: 


"a 


se d’a +2 
a= 30Ex ä 
Ri) ex r ß) Ex Ö Ey 


[9 9x 


 \der 


RN. #r 2 (6). 
Öfx : x Ö8y+-- x Syaut|örct:- ) 
OEy ! oy x 3 


—=00,08:4+''0T,ydyxy+t 


Wir machen nunmehr für a, einer Funktion von &x,..yxy; --, die Annahme, daß diese zweite 
Variation ö°a für alle möglichen e und y stets größer als null ist (stabiler Werkstoff). 

Um die Bedeutung dieser Annahme durch ein Beispiel zu zeigen, nehmen wir an, daß 
bei einem gezogenen Stabe von der Länge 1 und vom Querschnitte 1 in a nur die Zerrung ex 


\ Leo da ü ; 
vorkommt; ebenfalls folgt hieraus nur die Spannung %=7, Der Zusammenhang zwischen 
x 
0% 6x und e, sei durch Bild 1 gegeben. Die Zerrungsenergie a 


ist die Fläche unterhalb der Spannungszerrungslinie. Geht 
man von der Stelle e, zur Stelle ex+d:,, so ändert sich 





























‚ o, um Öo,. Das kleine Rechteck Ö*,:50, ist die zweite 
Variation. Ist diese positiv, siehe Punkt 1 des Bildes 1, so 
wird einer Vergrößerung der Dehnung eine Vergrößerung 

| der Kraft entsprechen. Ist sie hingegen negativ, Punkt 2 
|| des Bildes, so entspricht einer Vergrößerung der Dehnung 

Am TT T €, eine Abnahme der Kraft. Ist die äußere Kraft hingegen 
Ötr der konstant, so kann kein stabiler Zustand eintreten. 


Bild ı. Ist ö°a größer als null, so gilt dieses auch für 6°? A und 
für ö° II; die potentielle Energie wird zum Minimum: 


pa 
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A \(X\u+Yv+ZmwdV—-\(Zu+fAv+Zn)dog=Min . . ....M. 

v OR 

" 


Das heißt, vergleichen wir Verschiebungszustände, die mit den gegebenen Randverschiebungen 

(4). im Einklange sind, so wird der Ausdruck, der auf der linken Seite der Gl. (7) steht (die 
potentielle Energie), für die „wahren“ Verschiebungen zum Minimum. 
Das ist der Satz vom Minimum der potentiellen Energie. 


on ’ 
wi III. Eindeutigkeitsbeweise. 
en 
ine Wir machen für die Verschiebungen den Ansatz 
f w=au, +Pfu, +yw+': 
ft: 
pe- v=av, +Pv, +yvV, +: 
eT- n=aw, +ßn,+yn,;,-+':. 
ler i Fr bu “ 
. Hierbei soll jeder Teilansatz «;, v;, mw; den gegebenen Randbedingungen genügen. Nimmt 
ei man a+ß+y+::-=1, so erfüllt auch der Ansatz «, v, w die Randbedingungen. 
er Nun berechnen wir A, \ (Xu + Yr+Zm)dV, \(Zu+Her+Zw)dog und hieraus die 
m v og 
potentielle Energie. Aus der Bedingung, daß diese Energie ein Minimum sein muß, kann 
u man nach dem Ritzschen Verfahren «a, ö, y... bestimmen und findet für «, v, w die opti- 
malen Ausdrücke. 
en - a z en i 4 . : Q 
Hierbei erhält man für «a, ß,y... nur einen Satz von Werten. Nehmen wir an, es gäbe 
- wenigstens zwei. Die hieraus sich ergebenden Verschiebungen wollen wir u, ), 2, tw,» 


und u,PV, 0,0), w,®P) nennen. Wir bilden 
u= au," + (1- a) u,eP®, 
>). v=au, 9 +Hl— au, ®®, 
w= amw,’! +1 a) w,PY, 


Für a=1 und a=0 wird die hieraus gebildete potentielle Energie zum Minimum. 


= Ändert man a kontinuierlich von O0 bis 1, so muß hierbei die potentielle Energie // einmal 

zum Maximum werden. Dieses ist aber nicht möglich, da ö°?// überall größer als null ist. 

Auch die allgemeine Lösung für //=Min, ö//l=0 gibt nur eine Lösung. Gäbe es 

zwei Lösungen, so könnten wir aus beiden eine Linearkombination bilden, wie oben aus «,'®P' 
1d und «,PV; für jede Einzellösung wird // zum Minimum. Der weitere Beweis, daß dies nicht 

möglich ist, wird wie oben geführt. 

IV. Satz vom Minimum der Ergänzungsenergie. 
) Nun wollen wir die Spannungsenergie definieren. Die Annahmen der vorherigen 
® Abschnitte II und III bleiben bestehen. 
Die Zerrungsenergie war eine Funktion der Zerrungen: a=4(&x...yxy; ..); die Ab- 

leitungen gaben die Spannungen als Funktion der Zerrungen: 

te da ( z 
3 K=i, Ty- und zyklisch. 
8x zu ’ 
12} , BEN? i i r . . 
i Aus diesen sechs Gleichungen sind umgekehrt die sechs Zerrungen als Funktionen der 
x k ü 5 R 
Ä Spannungen bestimmt: &,=£x (9, T), Yry=Yx„(0,r) und zyklisch. 
n Ä Wir bilden nunmehr folgende Funktion der Spannungen 
a g b(0x,..Ta- )=0r8x 4° + TeyYyt'’— U lex, Yryı-.)- 
ıt 3 ’ U: ’ R : i a u 
h s Alle Zerrungen werden hierin durch die Spannungen ausgedrückt. Bilden wir die Ableitung 
a H nach einer Spannung, so erhalten wir 
o ob Ö 
= Ox&xt''TeyYx a a) 
g 06, 5, xcx xuYxyTt 
9 e 
Öfr OYxv da dEr 0a Oyry 
g { —=]-e, +0, = t + Tyy Fr Aus3 . = te. - ZIEN 
n "00, "00% O8, 00x Oyxy 00x 
—£,(o,r) und zyklisch, 

d ob 


—=y.,„(o,r) und zyklisch. 
Ri) Tr y IE ) . 
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Haben wir, als Beispiel, nur eine Zerrung &, und eine 
Spannung o,, so sind a und b im ox, — e, Diagramme nach 
Bild 2 Flächen, die sich zum Rechtecke &, 0, ergänzen. 
Im belasteten Körper stellt sich ein bestimmter 
Spannungszustand ein mit den Spannungen 0o%,..Txys +», 
a die Funktionen von «, y, 2 sind. Wir betrachten einen be- 
nachbarten Spannungszustand, der mit den Massenkräften 
X, Y, Z und den variierten Oberflächenkräften Z+6%, 
H+öH, Z+ödZ im Gleichgewichte ist. Es ist aber nicht 
€&r notwendig, daß dieser Spannungszustand für diese Kräfte 
Bild 2. die „wahre“ Lösung darstellt, d. h. daß er das gegebene 
Spannungs-Zerrungsgesetz befriedigt. 





Für die Variation der Spannungen gelten dann die Gleichgewichtsbedingungen im Innern: 


0 00x 
0x 07] 


Oörzy | ddrz, 


x —=0( und zyklisch, 
02 . 


und die Gleichgewichtsbedingungeu für das Oberflächenteilchen do: 
dZ2do=d0,docosa+ ÖrT,„docosß-+ Ört,,docosy und zyklisch. 


Wir bilden nun die Variation der Spannungsenergie für den ganzen Körper 


"t/0b ob _ 


—=||\\(er dog ++ yaydTayt: )dedydz 
1: 


„Allee Br, dv 


Fr na löny +: Jdzdgas, 


R Oy urY” | 

Durch partielle Integration erhält man bei Berücksichtigung der Gleichgewichts- 
bedingungen: : 
öB=\(Zu+öHr+dZm)do. 

[77 
Nun nehmen wir wie vorher an, daß für einen Teil der Oberflächenpunkte die Kräfte, für den 
anderen Teil die Verschiebungen gegeben sind. Für die inneren Punkte sind stets die Kräfte 
gegeben, so daß wir X, Y, Z nicht variiert haben. Die Variation der Spannungen nehmen 
wir so vor, daß dort, wo die Kräfte gegeben sind, die Variation von öZ gleich null wird, 
d.h. daß die abgeänderten Spannungen mit den vorgegebenen Kräften im Gleichgewichte sind. 

Wir erhalten 

öB \ Zu+öHer+ödsZn)do,=d. 
®p 

Die zweite Variation ö?B der Spannungsenergie wird gleich der zweiten Variation der 

Zerrungsenergie, da ö?b gleich ö’a, vgl. Gl. (6), wird. 


eb=do,de, + --zyöyzt Sa... . 2... WM. 
Falls diese, wie angenommen ist, stets größer als null ist, wird 


B-\(Zu+Hv+Zn)do=Min . ......2.2... 
2) 

Das heißt, vergleichen wir Spannungszustände, die mit den gegebenen Massenkräften 
und Oberflächenkräften im Gleichgewicht sind, so wird der Ausdruck auf der linken Seite 
der Gl. (9) (die Ergänzungsenergie) für die „wahren“ Spannungen zum Minimum. 

Dies ist der Satz vom Minimum der Ergänzungsarbeit (Satz von Castigliano). 

Beide Minimalsätze lassen sich nur für kleine Drehungen und Zerrungen ableiten, da 
nur dann die zweite Variation der potentiellen Energie bei stabilem Werkstoffe größer als 
null ist. 


V. Eingrenzung der elastischen Konstanten für besondere Elastizitätsgesetze. 
Mit Hilfe beider Sätze lassen sich Beziehungen aufstellen zwischen einer vorgegebenen 
Kraftgröße und der durch sie hervorgerufenen Verschiebung in Richtung der Kraftgröße, und 
zwar kann man die elastische Konstante, die in der betreffenden -Beziehung auftritt, durch 
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Näherungsansätze nach oben und unten eingrenzen. Hierbei kann das Gesetz für die Form 
änderungsarbeit etwas weiter als das Hookesche Proportionalitätsgesetz gefaßt werden. 
Wir nehmen an, daß a (ex... Yxy,..) eine Funktion der Zerrungen ist, der Art, daß bei 
Multiplikation aller Zerrungen mit » >0, die Zerrungsenergie sich um das n' *“-fache ver- 
größert. Dann wird b(ox,..Tx,...) eine Funktion der Spannungen der Art, daß bei Multi- 
plikation aller Spannungen mit m >0, b sich um das m’ *"“.fache vergrößert (a und b — 
homogene Funktionen der Zerrungen bzw. Spannungen). 
Beispiel: a=(, ex! +°+ 0, le," +++: 0, |yaul' +“ +, wobei a sowohl als auch 
Os Ogas-. Cis,.. positive Konstanten sind, wobei die € für positive und negative Zerrungen 
verschieden sein können. 
Ox 


Hieraus 0%|=(1+a) C,,lexl|®; >—v, 


Ex 
a [ox :(l1+a) A ud 
h= Cr Ox|' ruale.. Ce Teyl' rilan.. 


. u 6 a Suxr 
mit PR (+ aji + Ta 6, Ta USW. 

Die zweite Variation von a und b ist bei diesem Beispiele positiv. 

Für einen belasteten Körper, z. B. für einen 
Balken mit der Einzellast P nach Bild 3 werde der 
Zusammenhang zwischen dieser Kraft P und der Ver- 
schiebung € des Angriffsquerschnittes in Richtung der 
Kraft gesucht. Die „wahren“ Werte wollen wir mit E ) 
P* und C* bezeichnen, da wir es auch mit Näherungs- __ Eu Y- 
ansätzen zu tun haben werden, die Näherungswerte e 
für diese Größen enthalten. 

In den Stützquerschnitten und im Angriffsquer- 
schnitte der Kraft P seien Wände eingebaut, die so weit starr sind, daß sie sich nur in 
Längsrichtung des Balkens verwölben können, ohne daß ein Widerstand auftritt. Auf diese 
Weise erhalten wir eine einheitliche Verschiebung des Angriffsquerschnittes in Richtung der 
Kraft P. 

Wir machen nun bei gegebener Kraft P* einen beliebigen Ansatz für alle Verschiebungen. 
Für die Stützquerschnitte schreiben wir u=r=0 vor, für den Angriffsquerschnitt «—=0, 
Alle Verschiebungen drücken wir durch ihr Verhältnis zu Ü aus, so daß wir außer € 








Bild 3. 





y 


bet, 

| u 
nur noch die Größen r 
Verschiebungsansatze finden wir die Zerrungen und daraus die Zerrungsenergie 


v\ (w N ' z ’ 
); (7) (7) haben, die Funktionen von x, y, 2 sind. Aus diesem 


u 
1-+ı ’ Ü 1-4 ’ 
A=C!tele Jar dedydz=C'+t« 4". 


r 
A’ ist durch Integration aus dem Ansatze zu finden. 


Hieraus erhält man die potentielle Energie 


A \(Eu + Hv+Zn)dog=(!te4’— P*C. 
PR 


i e E u ® w : zus . \ 
Wir wollen annehmen, daß der Ansatz für (77 s (7): (77) nicht den wirklichen Verschie- 


bungen entspricht. Dann ist die potentielle Energie des Näherungsansatzes stets größer als 
die wirkliche potentielle Energie. Die Größe C ist noch frei wählbar. Wir wählen sie so, 
daß die potentielle Energie möglichst klein wird, um der wirklichen potentiellen Energie 
möglichst nahe zu kommen. Wir machen €" *« A’ — P*C zum Minimum und erhalten für den 
optimalen Wert von €: 

Pr" 


Cop * = (1-+a) 4’ 


Für diesen Wert wird die potentielle Energie: 


a e ap*'+tıle 
Br  . ° 


(era — P* Omi = ya Con Sn 237 4 
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Beim Einsetzen von (,,. ist im mittleren Ausdrucke 4’, im letzten Ausdrucke ( beseitigt. 
Wählt man für den Verschiebungsansatz die richtigen Funktionen, so erhält man die 


wirkliche potentielle Energie; diese wird gleich C* P*, entsprechend dem mittleren 


a 
l+a 
Ausdrucke der Gl. (10), wobei C,, in C* übergeht. Die mit beliebigem Verschiebungsgesetz 
gefundene potentielle Energie ist stets größer als diese wirkliche. Wir vergleichen beide 
und nehmen hierbei den letzten Ausdruck der Gl. (10): 


a p*®: + 1/a a (" . 
(1-+ a)! + !!e J’tla= l+a PT. 
Hieraus 

P* 


<< 
(1a) 4’ * 


Das Gleichniszeichen gilt nur dann, wenn man an Stelle des Näherungsansatzes den wirk- 
lichen gewählt hat. 

Bei Verwendung des Satzes vom Minimum der Ergänzungsarbeit nehmen wir (* als ge- 
geben an. Dann wählen wir für die Kraft P einen Spannungsansatz, der mit P und den 
Lagerkräften im Gleichgewichte ist. Hieraus berechnen wir 


BR ei en EA 


0x Try 


pP’ p ..)dedydz= P!+t\ueP', 


B=-\\\bdrdydz P!+!a\\\d 
In a 
B’ ist ein aus dem gewählten Spannungsansatze durch Integration berechneter Wert. Die 
Ergänzungsenergie wird: 
B (Zu+Hr+Zw)do,=P' rag —- PR”, 
pn 

Wir. wählen nun P so, daß die Ergänzungsenergie zum Minimum wird. Dieser Wert von P 
sei Pop. Man erhält: 


a C* 
Pop = i. 
dt (1-+a) B 
Für diesen Wert wird die Ergänzungsenergie 
(P! +1la B' P(*) a 1 p (* a: (*i ta 19, 
mın l+a opt (di + a)* ra B« . . . . ° a). 


Im mittleren Ausdrucke ist B’, im letzten Ausdrucke P,,: beseitigt. 
Wählt man für den Spannungsansatz die richtigen Funktionen, so erhält man die wirk- 


liche Ergänzungsenergie; diese wird gleich P* (*, entsprechend dem mittleren Aus- 


1 
A+e 
drucke der Gl. (12), wobei P,,» in P* übergelit. Die mit beliebigem Spannungsansatze ge- 
fundene Ergänzungsenergie ist stets größer als die wirkliche. Wir vergleichen beide und 
nehmen hierbei den letzten Ausdruck der Gl. (12): 


a«lc*!tra j 1 


Urt iger °©- 
Hieraus 
<pe(1+) PR Be ar a 
\ [71 


Durch die Ungleichungen (11) und (13) ist das Verhältnis An eingegrenzt: 
= ip 
darge pslitz) Be. 


Für a=1, also für das Hookesche Gesetz entsprechen die Ergebnisse meinen früheren 
Untersuchungen ?). 


VI. Annahme eines allgemeinen Elastizitätsgesetzes. 


Um die Verhältnisse für ein beliebiges Elastizitätsgesetz für a (jedoch mit d’a >0) zu 
klären, untersuchen wir das bisherige Gesetz, wobei a eine homogene Funktion der Zerrungen 
war, etwas eingehender. 


2, Z. angew. Math. Meeh. Bd. 11 (1931), S. 244 bis 245 und Bd. 18 (1938), S. 375 bis 379. 
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Die Ungleichung für einen beliebigen Ansatz für die Verschiebungen lautet mit be- 
liebigem €: 


r a gu . 
er 2 -PcC2 ., 
.+« 
R P* J1j« 
Setzt man C=C,.=|,— ‚, ein, so erhält man 
(1+ a) 4') 
P* 
(1+a) 4’ 


Wir können aber auch ('=(* einsetzen und erhalten 


2 Be a N. 
Fi + a 4' P* (* (* P*. 
l+a 


Zu beiden Seiten fügen wir P* (* hinzu und erhalten wie vorher: 
P* 
(l-+a) 4' 
Bei einem beliebigen Elastizitätsgesetze mit ö°a >0 erhält man verschiedene Ergebnisse 
für C=C. und C=(* 

Wir lassen die Kraft allmählich von 0 bis P* wachsen und bezeichnen einen dazwischen- 
liegenden Wert mit p*. Hierbei ändert sich die Verschiebung in Richtung der Kraft von 0 
bis C*; ein dazwischenliegender Wert sei c*. Dabei ist p* eine Funktion von c* oder c* eine 
Funktion von p*. Für die wahren Verschiebungen wird das Integral A=\\\ad.rdydz gleich 

c F 
der hinzugeführten Formänderungsarbeit | p*de*. 
0 

Wenn wir die Kraft P* wirken lassen und einen beliebigen Verschiebungsansatz mit 

der Verschiebung €‘ machen, so erhalten wir die Ungleichung: 


(* P* 
A(C) P*C "\p*de* P=0r \e*dp* ae ar > 


Hierin ist A eine Funktion von €. Das € wird zum Optimum, wenn die linke Seite zum 


Minimum wird, also für 
d Ale) 


a ee le rt ED ER 
d( opt 

Hieraus ist (',,.. als Funktion von P* bestimmt. Man erhält: 

P* 

a. EEE 4 


0 


AAO) 
A (Cop) + Cop ER 
(Vopt) T Vopt IC 
Setzt man hingegen in Ungleichung (14) C=(* ein, so erhält man: 
Ay 
ZEIEIETE 2% Er er AT 
o 
Ungleichungen (16) und (17) besagen etwas Verschiedenes. In Ungleichung (16) kann man 
P* 
die linke Seite als das Integral | e,,.dp* auffassen, wobei c,,. eine Funktion von p* ist, ent- 
oO 


da (Copt) 


dA (Copr) 


sprechend der Gleichung \\\ den dedydz dem —»p*. Denn es ist dann 

P- P* ' 

Ra | dp* E d’ A (Cop) 

Var = an en denen "erde 

C pt 
dA(CH dA (Cop dA 
- \ opt) : \ ( ( } a Copt en CO - A (Co) . 
d Copt A Cop d Copt 


0 


Dieser Ausdruck entspricht der linken Seite der Ungleichung (16). 
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Wir können über p* sowohl cp: als auch c* auftragen, Bild 4. Die Größe p ist in 
P* 


diesem Bilde nach oben, die Größe ce nach rechts gerichtet. Die Fläche | c*d p* 
5 ö 

wird nach Ungleichung (16) stets größer als die Fläche \c,,: dp* für gleiche obere Grenze P*. 
oO 


In Ungleichung (17) kann man die linke Seite wie folgt schreiben: 


c* 
*\ Y% k Ö A (c*) I „* 
\\lac )dedydz=A(C )=\ Jet de*. 
pe ; 
RER : 8 Ale*) Mi BWRR- i 3 d = h 
[rägt man p* und p: Sc über c* auf, Bild 5, so ist nach Ungleichung (17) die 
C* C* 
Fläche \ pdec* stets größer als \p*de*. 
In den Bildern 4 und 5 
sehen die Linienzüge gleich 7 






aus. Es sind aber verschie- 
dene Größen und Flächen, pe 
diebei den Ungleichungen (16) 
und (17) verglichen werden. 

Entsprechend liegen die 
Verhältnisse beim Minimal- 
satz der Ergänzungsenergie. 
Wir schreiben die Verschie- 
bung C* vor und lassen die 

















Kraft P unbestimmt. Dann Be] c 
erhalten wir Bild 4. Bild 5. 
P* (* 

\\\b(P)dedyd2- P*=B(P)-—P*z\*dp - *P=—-\p*de*. 
v° 0 o 

i „. dB(P,..) ; 

Für P=P,p gilt d : "—=(*, Hieraus 

opt 


© 
AB(Pay) _ | 


B (Pop) + Pop P,. <\» 0 BEURÄF FE | 
Für P= P* erhält man 
P* 
et. A Pa rn | 
© 
ge r i j Fu i 2 N BB 
Die linke Seite der Ungleichung (18) schreiben wir | p,,:d c*, wobei =e 


ö dp.p 
c* 
Trägt man p* und psp über c* auf, so ist nach Ungleichung (18) das Integral | p,,: d.* 
c* c* z 
stets größer als \ p*dce*. Im Gegensatz hierzu war nach Ungleichung (17) \pde* stets kleiner 
c* : 
als \p* de*. 


Mus r ’ dB(p*) 
Trägt man nach Ungleichung (19) e= d ” 
p* p* . p* 
\cedp* stets kleiner als \c*dp*. Im Gegensatz hierzu war nach Ungleichung (16) | co dp* 


0 o o 


p* 
stets größer als \c*d p*. 


0 


und c* über p* ab, so ist das Integral 


Hiermit ist zur Bestimmung der elastischen Konstanten noch nichts gewonnen, denn 
aus der Ungleichung der Integrale folgt nicht eine entprechende Ungleichung der Integranden. 
So folgt nieht aus Ungleichung (19), daß bei diesem Ansatze € < C* ist. 


*% 
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ın Man könnte nunmehr annehmen, daß bei einer geschickteren Verwendung der grund- 


legenden Annahmen es vielleicht doch möglich sein könnte, eine unmittelbare Beziehung 
zwischen €, C* und C,,.: bzw. P, P* und P,,: zu finden. 
Daß dieses nicht der Fall ist, soll an einem Bei- 


























e spiele gezeigt werden. Die Kraft P* ziehe an zwei j Fe 
Stäben, Bild 6, wobei infolge der Kreuzkopfführung die ae ih 
Verlängerungen beider Stäbe gleich groß sind. Beide l 
Stäbe haben die Länge 1 und den Querschnitt '/.. Für er 
die linke Hälfte des oberen Stabes laute das elastische / m 
1 A as? auge 
Gesetz für diese Belastung aq,=75 &?, für die rechte 9 — 
die Hälfte a, = —&,', für die linke Hälfte des unteren Stabes Bild 6. 


sei a,= z &,, für die rechte Hälfte a, = J Raı 5 
. PR 
Die Spannungen werden 0, =e,, ,=&', 1,=&,, 

o, = &,'!*; Spannungen und Zerrungen in Richtung V 
der Kraft P. 

Erst werde der wirkliche Zusammenhang 
zwischen P* und C* ermittelt. Wirkt im oberen DA, als AruonC® 
Stabe die Kraft P,, so wird die Spannung 2 P, > 
und hieraus die Verlängerung 












2Py, als Fkt en” 


N 
*=>;[P,+@P,)"']; 


2A0ls fkt wnC*\ 

















wirkt im unteren Stabe die Kraft P,„, so wird Copr 2ls Pkt von Tue \c*ols Arvonn“ 
i _ >P i ’erlä . oder 
- die Spannung er und die Verlängerung I. 
0 == > [2 P, + (2 P.)1- Teoor als Fkt vonc* 

Aus beiden Gl. für C* folgt P,=P,.*' und C*=P, 

+P„=P*. Im Bilde 7 sind 2P,.2P,„ und ihr (Azar] Le 

Mittelwert P* als Funktionen von (* dargestellt. | 

Bild 7. 
Nunmehr machen wir den Näherungsansatz 
s—=czt: € IT; e,=c; 4Atld> —c’+ c® — (e? +8 [ad 
Mc Mn re Da LE Berl 5 

5). . u \ i ee 

Der Satz vom Minimum der potentiellen Energie gibt: 

P* 
1 (Cop) — P* com — \c*dp* 
4 ‚opt opt <= y’ ] 

9). ö 

mit 

AA (Coy) 1 
opt > f 4 1/a * 
7 l&Copt + Cop Cop "|>=P". 

st. d Cop: { [2 Cop: pe + Cop] 1 


Für jedes p* erhält man ein c,p:, das in Bild 7 über p* aufgetragen ist. Man sieht, 
* P* Pr 


C % . 
daß hier \c„dp*=\ec*dp*, trotzdem aber c,p nur im unteren Teile kleiner als c* ist. 
0 0 
er Machen wir den Spannungsansatz 0, =0,=0,=0,=p, so wird 
1/1 4 1 1 
[3 FraB 3 5 & 
B(P)=T spPtrzp "+yrtgr 
l ‘ Der Satz vom Minimum der Ergänzungsarbeit gibt: 
y* i B (Popt) = Popt C* m \ p* d c* 
} mit 
dB (Popt) l i 
s . — [2 Popt + Pont! + Pop] = c*. 
R d Pop 4 [2 Port + Pont" + Popr"] 
n % 
1. Ä Für jedes c* erhält man ein pop, das über c* aufgetragen ist. Auch hier wird psp nur im 


unteren Teile kleiner als p*. 
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VII. Eingrenzung der elastischen Konstanten für den allgemeinen Fall. 


Für diesen allgemeinen Fall des elastischen Gesetzes möchte ich nunmehr einen Weg 
angeben zur Eingrenzung der elastischen Konstanten, d.h. zur Aufstellung von Ungleichungen 
zwischen P* und (*. (Die Genauigkeit wird hierbei jedoch geringer als bei dem früheren 
Verfahren für die Fälle, wo dieses durchführbar ist.) 


Wir addieren die Ungleichungen (16) und (18) und erhalten: 
P* C* 
A'(Copı) + Com P* — B(Pop) + Po CS |c*dp* + p*de*= P* C*. 
Hieraus 


Aid - En P’ + Bin) - Pu” +eP20.. ...... (MM. 


Y* 


In der linken Seite sind C,» und P,)» Funktionen von P* bzw. ( 
Gleichungen: 


entsprechend den 


d A (Copt) B > 9 
1 ea ee N Bpeenge, 
d B (Popt) 2% vr 99 
Fu re 


Wir denken uns (,,. und P,, nach Gl. (21) und (22) in Ungleichung (20) eingesetzt, so 
daß wir eine Ungleichung für P* und C* erhalten. Diese Ungleichung ist nur in einem be- 
stimmten Gebiete der C* — P*-Ebene erfüllt. Um die Grenzen des (zebietes festzustellen, 
lösen wir das Gleichnissystem: Ungleichung (20) mit dem Gleichniszeichen, Gl. (21) und 
Gl. (22). 

Ebenso könnte man die Ungleichungen (17) und (19) addieren. Die hieraus entstehende 


Ungleichung ist jedoch immer erfüllt. 


Für das Beispiel von VI werde die Eingrenzung der Beziehung zwischen P* und (* 
durchgeführt. Ungleichung (20) lautet mit den früheren Ansätzen: 


1 1 - 4 £ ” e 
[Pop + = Popt” le = Popt?'* Popt (* + (* P* >—0 j 


( 2 1. .$ C 5/4 Ber x p*rı 
opt T 5 “opt 75 Lopt opt 4 


1 
t 
Gl. (21) und (22) geben: 
| ’ 1 x 
4 [2 Copt + Copt' T Copt''*] =P, 4 1? Popt + Popt* + Popt” *] =UÜ. 


Setzen wir (* und P* in die Gleichung ein, die aus 
der Ungleichung zur Festlegung der Grenzen folgt, so er- 
halten wir pe 


RrUE. RER 1 
[Con #e 5 Cop + 5 Copt” ı 
£ B | 





Pop’ + 5 Popt” f B Pop“) 


1 eg 
2 = | Cop + 2 Copt* ra Copt''* 











I 
a 


| l 1 

| f f 

| Popt + = Popt’ + > Popt' =V, 
| ö & 








Für jeden Copt- Wert erhält man zwei Popt- Werte und um- 
gekehrt. Dasselbe folgt für die Grenzen der C*. und 
P*-Werte. Bild 8 zeigt das Ergebnis; innerhalb des Bild 8. 
schraffierten Gebietes muß die (*- P*. Linie verlaufen. 227 
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Das Knicken schwerer Gestänge. 
Von Fr. 4. Willers in Dresden. 


Während man gewöhnlich bei der Bestimmung der Knicklast eines Stabes 
sein Eigengewicht vernachlässigt, ist die Berücksichtigung dieses (rewichtes 
bei der Untersuchung des Knickens langer Gestänge, etwa von Bohrgestängen. 
wesentlich. Ausgehend von dem bekannten Minimalprinzip erhält man auch 
hier Gleichung und Randbedingungen. Die Integration führt auf Bessel- 
sche und Lommelsche Funktionen, aus deren asymptotischer Darstellung 
die Knickkraft für sehr lange Gestänge gewonnen wird. 


1. Einleitung und Ergebnisse. Bei kurzen Stäben wird das Eigengewicht die Knick- 
kraft kaum beeinflussen. Anders bei langen Gestängen, für die im folgenden ganz allgemein 
die Gleichung für die Knickkraft in Form einer Determinante gegeben wird. Diese Kraft 
trägt einen Teil des Gestänges von der Länge b, der also auf Druck beansprucht ist, während 
der weiter oben liegende Teil Zug erfährt. Rechnerisch ergibt sich aus dieser Determinante 
die Knickkraft, wenn man die Länge des Gestänges in Vielfachen oder Bruchteilen der Länge b 
in sie einsetzt. Es ist also nicht ohne weiteres möglich, für eine gegebene Länge ! die Knick- 
kraft zu berechnen, sondern man muß zunächst für vorgeschriebene Vielfache von b die Knick- 
kraft P bestimmen. Aus dieser ergibt sich mit den Konstanten des Gestänges die zugehörige 


> 


YET(ug) als 


Länge ! und man kann nun etwa wie in Bild 1 die dimensionslose Größe P- 


8 / 
Funktion der ebenfalls dimensionslosen Größe i=ı| = auftragen (E Elastizitätsmodul, 


I Trägheitsmoment des Gestängequerschnittes, u. g Gewicht der Längeneinheit des Gestänges). 
In Bild 1 sind so Punkte für verschiedene Aufhängungen und verschiedene Lagerungen am 
oO e fe) je} 


o 4 2 3 RE 7 8 
I I | 








Bild 1. Die Knickkraft schwerer Ge 
stänge als Funktion der Länge. 








I: Oben eingespanntes, unten freies Ge- 
stänge, 





II: oben und unten drehbar befestigtes 
Gestänge, 

IlI: oben eingespanntes, unten drehbar 
befestigtes Gestänge, e 

IV: oben drehbar befestigtes, unten ein 
gespanntes Gestänge, 

V: oben und unten eingespanntes Ge 
stänge. 




















unteren Ende eingetragen und durch Kurven verbunden. Auch sind dort die zur l-Achse 
parallelen Asymptoten dieser Kurven eingezeichnet. Für Gestänge von einer Länge, die 
größer als 2b ist, ist die kritische Kraft höchstens um 12 bis 15° größer als für sehr lange 
(Gestänge. Ferner sieht man, daß es, je länger das Gestänge wird, um so weniger ausmacht, 
ob das obere Ende eingespannt (Kurve I, III und V) oder aufgehängt ist (Kurve II und IV). 
Für kurze Gestänge ist das Eigengewicht gegenüber der Knickkraft nur von geringer 
Bedeutung. Bildet man 

a Pl 

Fra, 

E-I 
so ist das gerade die Größe, die in der Eulerschen Knickformel im Falle I (oben eingespannt, 
nn“ 
‚ 


unten frei) gleich 1 im Falle II (oben und unten drehbar befestigt) gleich °, im Falle III 


und IV (einerseits eingespannt, andererseits drehbar befestigt), die im Falle der Vernach- 
lässigung des Eigengewichtes, unabhängig davon, ob die Einspannung oben oder unten 
stattfindet, zu der gleichen Knicklast führen, etwa gleich 4,49? und im Falle V (beiderseits ein- 
gespannter Stab) gleich 4° ist. In Bild 1 sind die entsprechenden Kurven eingetragen; man 











Z. angew. Math. Mech. 


44 Bd.21 Nr. 1 Feb. 1941 


Willers, Das Knicken schwerer Gestänge 





f ? A ur Z 
sieht vor allem bei den beiden ersten Kurven, für die noch die Knickkraft für — b berechnet 


ist, wie diese Kurven die Kurven der Knicklast für kleine Längen sehr gut annähern; und 
zwar müssen unsere Kurven, wie man sich leicht überlegt, immer größere Werte für 
die Knickkraft ergeben als die asymptotischen Kurven P P=c, so daß sie in Bild 1 unterhalb 
dieser Kurven verlaufen. 


I. Aufstellung der Gleichungen. 


2. Drückt man einen senkrecht hängenden schweren, geraden Stab durch eine an seinem 
unteren Ende in der Stabrichtung wirkende Kraft P so weit zusammen, daß die neutrale Gleich- 
gewichtslage gerade etwas überschritten wird und hält dann die Enden fest, so leisten beim 
Ausknicken die äußeren Kräfte kein@ Arbeit; es findet dabei nur eine Umlagerung der Form- 
änderungsenergie statt. Die zur Krümmung des Längenelementes dx erforderliche Energie 

zw 
. EI E n 


2 
d =) dx wird dadurch gewonnen, daß sich dieses Element beim Ausknicken wieder 


2 


\d« frei wird, wo 


dw dw 


/ 2 
um den Betrag (y 1+ () 1)de dehnt, so daß die Energie 5.u(s = 
II die an dieser Stelle wirkende Druckkraft ist. Dabei ist angenommen, daß es sich beim 
Ausknicken nur um kleine Formänderungen handelt, so daß sich die Endpunkte des Ele- 
mentes dx senkrecht zu der mit der geraden Stabachse zusammenfallenden x-Achse um das 
kleine Stück w verschieben. Der Unterschied der Formänderungsenergie für die gekrümmte 
Lage gegen die nichtausgeknickte ist also, wenn man / konstant nimmt, 
j I 


Eis, A alfa 


J 2 dx 


u u 
Die »-Richtung soll dabei positiv nach oben gerichtet sein und der Nullpunkt in das untere 
inde des es gelegt werden. Somit ist, wenn man die Stabmasse pro Längeneinhei 
Ende des Stabes gelegt 1 Somit ist die Stabma >» Län inheit 
mit « und die am unteren Ende angreifende Druckkraft mit P bezeichnet 


I=P-—ug«. 


Soll die ausgeknickte Lage eine Gleichgewichtslage sein, muß 














I 1 
El /d?w\? L{ (dw\? 
J= > \(4) dx 7 \ (P-ug 2) dxe=0 (2,1) 
0 0 

ein Minimum sein, d. h. es besteht die notwendige Bedingung, daß 

I 

„(dw döw 
n. \aa' da? ©* 

yon : 

: lm dö 

dw döm 
P\ (P— ug) 2 dx 
I) 
r (2,2) 
‚„„?wdön! ‚„d’w ! „(dw 
.. dx’ dx& o u. das? "L en ı\ dat? le 
2 —0 
of | 1 
dw ( dw 
> ö re > PR ö 
(I 4IE) ir on ey. \a.,t ng %) 25 öwda 
1) 





werden muß. In bekannter Weise ergibt sich daraus zunächst, durch Nullsetzen des Inte- 
granden, die Differentialgleichung des Problems 


d’w 1 4] 


dw 


.£ Ber 2; 
de *E1da | “900 a) 
Setzt man hier zur Abkürzung 
u 4 an 5) 
Br und we —a 2,4), 
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wo b die Länge des unter Druck stehenden Stabteiles ist, während a die Dimension 1:1” hat, 
und führt eine Integration aus, so lautet die Gleichung 


d’w dw 
Br h 95 
dm tr (b— r) iz" NE ar Dh a a A 


Die Randbedingungen liest man ebenfalls aus (2,2) ab. Als geometrische Randbedingungen 
ergeben sich 


1. für den nicht verschiebbaren Endpunkt v—( 
2. bei unveränderlicher Endrichtung = =0 
und als dynamische Bedingungen ! 26). 
3. falls keine Endmomente wirken .. 0 
4. falls am Ende keine Querkräfte wirken ee +a(b— x) eo 





Ist ein Stabende freiverschiebbar, so ist daher in (2,5) die Integrationskonstante k gleich 
Null zu setzen, so daß man in diesem Fall mit der verkürzten Gleichung auskommt. Will 
man sich auf die Bestimmung der Knicklast nicht zu langer Stäbe beschränken, kann man 
für die Lösung von (2,5) Potenzreihen mit unbestimmten Koeffizienten ansetzen und durch 
Einsetzen in die Gleichung Rekursionsformeln für diese ableiten. Da hier auch asymptotische 
Ausdrücke für den Fall sehr langer Gestänge gewonnen werden sollen, empfiehlt es sich, um 
auf bekannte Funktionen zu kommen, die Substitution 


. - r 
S=yyaylb- x)’ und "uud 10 Be ee 2,7) 


zu machen. Es ergibt sich dann nach einfacher Umrechnung 


dp 1d ) RR 
+rarteli-ggeV or: 88 


a z = m. re z = 
ar 2-08 92 9a? fa 


II. Das an einem Ende frei verschiebbare Gestänge. 


3. Gestänge endlicher Länge. Da in diesem Falle e,=0 ist, kann man die Lösung 
der verkürzten Differentialgleichung. einer Besselschen Gleichung, sofort hinschreiben ; 
es wird 

o=c,J({)+c,J ı(Ö) 
6 3 
und damit 


dw ” - Re: 2 
da VII O+aV JH Mb re, 


. . ö r i 1 E53 R 
Tabellen für die Besselschen Funktionen mit dem Index — und — findet man z. B. in 


dem Buche von Watson'). Da man aber diese Funktionen noch mit Y Z multiplizieren muß 
und außerdem Ableitungen und Integrale gebraucht, ist es meist vorteilhaft, auf die Reihen- 
entwicklung der Besselschen Funktionen 

:=\9% 
” t,\“ 
» \ r (3) 


Zur: ---- +: 
k=u 


„@=6) 


zurückzugehen. Daraus ergibt sich, wenn man für { wieder den Wert aus (2,7) einsetzt und 
wenn man zur Abkürzung noch 





2 3 
e=Valb—a, dh.l=32.......... (8) 
einführt, 
dw, ı/2 1 2* 2’ 2’ in ’ 
a Era: 3.111.4138.511.4.7 3.,311.1.7.90 | kW. (3,4, 
3) 


1) G.N. Watson: A Treatise on the Theory of Bessel Funetions. Cambridge 1922. 


2) z.B. R. Weyrich: Die Zylinderfunktionen und ihre Anwendungen. Leipzig und Berlin 1937, S. 38. 
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und ebenso für 


6 „9 N 
. 


Ir.ataata 5 Far... Wh. 


de, v2 ı E 
dr r(Z)' 3- 


Nun braucht man bei der Behandlung des Problemes nicht nur die Ableitungen von w, sondern 
auch w selbst. Durch Integration ergibt sich aus (3,1) zunächst einmal ohne Berücksichtigung 


3 
9} 
der Grenzen, wenn man beachtet, daß j tde=- 3 „a: ist, 
8 .) E n 8 +) „ 
w= —( — J E dE e — \, CE) d£E ©: r P . . e 3,6). 
| AEG V 3a aOdlte (3,6) 


Mit diesen Formeln lassen sich die Knicklasten von an einem Ende freien Gestängen be- 
stimmen. Aus den Randbedingungen ergeben sich drei homogene Gleichungen für die c,, e, 
und e,, aus denen sich diese Größen nur berechnen lassen, wenn die Determinante ihrer 
Faktoren Null ist, und das gibt die Knickdeterminante. 

Weyrich behandelt so z. B. den Fall des unten eingespannten Gestänges, das unter 
seiner Eigenlast ausknickt®). Da sich das Gestänge gerade selbst trägt, ist am oberen Ende 
x=b, also z2,—0, am unteren «=0, also z2,=yab. Die Randbedingungen sind für 2,=0 
d’w ; ’ u 
ia 79 d.h. wie man aus (3,4) und (3,5) erkennt, wenn man beachtet, daß dz- Yadı 
de? 
ist, es muß c,—=0 werden. Am unteren Ende muß 


c,) JAHR) —(), 


8 


= 
©, — \ OD)dC+ec,=0 
| 3a pi . 


sein. Aus der ersten Gleichung findet man den Wert von £„ für den J_1({.)=0. Aus den 
6} 


Tabellen von Watson’) oder aus den Tafeln von Jahnke-Emde®) findet man „= 1,8663 
und damit naclı (2,7) und (2,4) 


! = 1,9863 | El und P=19863 YEI (ug). 
u 4 


Die zweite Gleichung, die hier nicht weiter interessiert, dient nach Einsetzen der Integral- 
grenzen zur Bestimmung von e,:6,. 

Zum Vergleich sei ein Gestänge behandelt, das oben eingespannt und unten frei ist und 
auf das in der Stabrichtung von unten her eine Kraft P drückt. Das Gestänge habe zunächst 


’ u ; E dw ’ Do : 
einmal die Länge b, dann wird für ,=0 w=0 und 2 —=0(,. Die erste Gleichung dient zur 


Bestimmung von e,:c,, die zweite ergibt, wenn man (3,4) und (3,5) beachtet, e,=0. Am 


BR ) 
unteren Ende wird 73 =0. Benutzt man die Reihe (3,4) und läßt den konstanten Faktor 


fort, ergibt sich 2, als erste Nullstelle der Reihe 


9 


12 


[07 


2° z°® 2 i 
3.11-1 '3°.21.1.4 3°.31.1.4.7'3%.411-4-7:0°° 


Durch Eingabeln findet man z,=1,53, und damit nach (2,4) 


a a 
’ 


2 EI A . sm 5 
I= 1,53 Pr und P=153yEldu VE 


Das ist eine wesentlich kleinere Knickkraft als im ersten Beispiel; eine Einspannung unten 
wirkt eben ganz anders auf einen schweren Stab als eine Einspannung oben. Für die gleiche 

} a a b 3b N 
Befestigung sind die Knicklasten noch für I=—, 5- und 2b berechnet. Die Knickkraft er- 
gibt sich dann aus der Determinante 


3, R. Weyrich:|.e. S. 131 ff. 


4, Jahnke-Emde: Funktionstafeln. Leipzig und Berlin (3. Aufl.) 1938, S. 167. 
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f» (2) 13 (20) 
N 5 . 
dz dz i 
8 1 3 1 8 3 . l l 
wo 2„=yab und z, bzw. —yab, 5yab, yab,.d.h. also gleich —2,, — 5 24, —2Zu 


ist. Durch Eingabeln erhält man für 2, die Werte 2,25: 1,245; 1,113. Da es sich hier nur 
darum handelt, eine Übersicht zu gewinnen, ist die Eingabelung nicht sehr weit getrieben. 
Die letzte Stelle der hier und weiterhin angegebenen Werte für 2, kann ungenau sein. In 
Bild 1 sind diese Werte markiert und durch die mit I bezeichnete Kurve verbunden. Aus 
dieser Kurve kann man jetzt zu einem ! des eingezeichneten Bereiches das zugehörige P ent- 
nehmen. 


4. Sehr lange, unten freie Gestänge. Man sieht aus Bild 1, daß die Knickkraft 
mit wachsender Länge des Gestänges sich sehr wenig ändert, so daß man für sehr lange Ge- 
stänge mit ausreichender Näherung den Wert von P für <>» benutzen kann. Die Werte 
von £ für <>b werden imaginär, während, wie aus (3,4) und (3,5) folgt, f,(z) und f, (z) reell 
bleiben. Man muß also in (3,1) bei YZ denjenigen Wert der Y i wählen, für den w reell 
bleibt. Man benutzt am einfachsten die Hankelschen asymptotischen Reihen für die Bessel- 
schen Funktionen°). Setzt man 








M—ı 
a E „(r, 2m) Be 
P,()= Y (-ıy Brent 5 ab: FE 
m—\ 
wi 
r 7 v‚2m+1) a 
Q,(8)= N ( een u u 12% mal cl EEE * 3 
m=(\ 
n ri 1 it Pr 1 b; , I 
so wird mit v= 7 bzw. 3 
SUR 2 Eu}, a Pen FOR ,. a e 
10=] „eos? 2) © sin (£ 3) 0) - Me © 
. _& u, VRR ee er r 
,0=] „eos (2 Pd sin(Q 75) 2,0) en an 


Es würde hier genügen, wenn man sich auf das erste Glied dieser Ausdrücke beschränkte, 


da aber nachher auch das zweite Glied gebraucht wird, seien für {=in hier gleich die beiden 
ersten Glieder hingeschrieben. Es wird 





: RE WO EN 
Jin) >| ain cos (in is) sin (i7 3) 66, | 
Sr * a ei = 
r 2 fh )e* z 
>| an2 67)° 


’ dw ir . ag 
um einen reellen Wert für 7; 5 zu bekommen, haben wir hier also fin=yne °® zu setzen 


und erhalten 


eo 5 2 E Ki. . 1 \ \ 
rl Yn | an 27 7 una yga)|ı 60) De or 
* in 
Ähnlich findet man mit Yin=yne® 
dw, , [a el I\ = 1 „ 
2 | > tt) not) a N N 


Somit ergibt sich als Knickdeterminante aus (3,7), wenn man nur das erste Glied von (4,6) 
und (4,7) benutzt und den gemeinsamen Faktor Yn9(n) fortläßt 


5), R. Weyrich:l.e. S. 47. 
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N + 


d : d e 
fz (Zu) dafs (2u) 


dz 


(4,8), 


2 1 2° 2 F 
3 
2.y2 32° 62° 92° 


sr) 3.2 732.91.2.5 3.312.5.8 


| 
._ 


Die Werte der beiden /-Funktionen kann man aus den Tafeln von Jahnke-Emde ent- 
nehmen®). Man findet dann den Wert z,=1,0188. Für sehr lange Gestänge ist also die 
Knickkraft 

P=1,0188 YE I(ug)°. 


Damit ist die zur l-Achse parallele Asymptote der Kurve I bestimmt. 


Ill. Gestänge, deren Enden nicht seitwärts ausweichen können. 


5. Gestänge endlicher Länge. Ist das eine Ende des Gestänges nicht frei ver- 
schiebbar, so wird die erste Integrationskonstante c, nicht gleich Null und man hat das 
Integral der vollständigen Gl. (2,8) zu benutzen. Nun pflegt man die Integrale der Gleichung 


[73 1 ’ er 1 - 
n ++ [1 Zln=ua Er \ . 


als Lommelsche Funktionen zu bezeichnen ’). Für sie hat man die Potenzreihenentwicklungen 


g Zen: Cu +1 = i E [u +3 [6 a 
am + (6): (u+1)’- » [ua+1)? — r?] [iu + 3)? — ) 18,2), 
Sud: Be ö ı (u — T y? sr [(u 1)? )- I - 3)” »?] ® | (5,3), 


- - 


die man nach der Methode der unbestimmten Koeffizienten finden kann. Zwischen ihnen be- 
steht der Zusammenhang) 


Su,v() = Su,» (E) + 24 ır(5 u Eee s)I(Su+2 + 5) | 
nr (5,4). 


un. 1 Er Ge 
|sinZz (u—r)aJ,() eos (ur) vol) 





’ ’ : 1 
Aus (5,1) ergibt sich die Gl. (2,8), wenn man «=0 und r—= 3 setzt. Man erhält dann 


unter Benutzung von (5,2), wenn man wieder die Abkürzung (3,3) einführt, 


dw, 3,/3/. 1, | A Eu, 
7 3 (2 0° 30.56? ha. 2.2... GB) 


Durch Integration findet man daraus ähnlich wie (3,6) 


> 


w,= | 30 \9.:0aL ver ee Se 


Hat man nun z. B. ein an beiden Enden eingespanntes Gestänge, so ergeben sich die vier 
Randbedingungen 


6, Jahnke-Emde:|.e.S. 14. 
?), Watson:].e. 8.345. Vgl. a. Nielsen: Handbuch der Theorie der Zylinderfunktionen. Leipzig 1904, S. 48. 


s) Watson:|]. e. S. 347; Nielsen: ]l. e. S. 228. 
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,w,(2) +92) +6,m() +0,=0 
5), S d d d 
i c, ds m, (z,) +0 2 te) 76 3; w,(2,) u 
; TE RR een, ars RE 
? d VRFCh 
© dz w,(zu)+ € 2 w,(Zu) TC, da”: (2.) — 
j em, (Zu) + 0,0, (2u) + 0, m, (Zu) ra, =0 
x Die Determinante der Faktoren der c, gibt die Knickdeterminante. Ist das Gestänge an 
i einem Ende nicht eingespannt, sondern drehbar gelagert, so treten hier an die Stelle der 
h ersten Ableitungen die zweiten. Das kann natürlich auch an beiden Enden der Fall sein. 
i Die Determinante läßt sich sofort in eine dreireihige umformen. Aus ihr sind unter Benutzung 
- i der Reihen (3,4), (3,5) und (5,5), ihrer Integrale, am einfachsten von z2=0 aus genommen, 
ie f 1 j 
> und Ableitungen die Knickkräfte für = b, b, 5b und 2b, also für ,=5 2.0, -— 5 2. und 
; -2, durch Eingabeln berechnet. Für die z, ergeben sich die folgenden Werte: 
i 3 
Länge des Gestänges. . .. 2». .2... a 6 h 5b 2b 
a beiderseits eingespannt . . . u E= 4,15 3,77 3.54 
R unten eingespannt, oben drehbar gelagert 3,18 3,54 3.39 
T- 4 unten drehbar gelagert, oben eingespannt 3,10 2,60 2,35 
as K beiderseits drehbar gelagert. ...... 3.59 2.64 228 215 
ng N 
N H Die Potenzreihen sind bei der Berechnung zum Teil bis z°' benutzt, doch ist die Rechnung 
: nirgends so weit getrieben, daß die. letzte Stelle garantiert werden könnte. Die Knickkraft ist 
& 8 I 
£ n " N 2/0 a a ag 1 kl El 
en ’ in Jedem Falle P=z,yEI(ug)’, die Länge des Gestänges , = 2, RR Me 
= . ‚ r zb 2 ug ug 
x ‘ e ‘ 
8 s/ 3 
3 El El 
-) ! Is = 57 2u bzw. I,,=22, . Die Punkte sind ebenfalls in Bild 1 eingetragen 
s 2 2 ug ug " u. 
N und durch die Kurven V, IV, III und II verbunden. 
3), . 
A 6. Sehr lange, an beiden Enden gelagerte Gestänge.. Für sehr lange Gestänge 
)e- i x £ e . > .. de, dm, 
; braucht man wieder die asymptotischen Formeln, deren erste Glieder für | ° und " 
” uUx 
bereits in (4,6) und (4,7) angegeben wurden. Um für s,.ı(£) aus (5,4) eine asvmptotische 
u 0» . I 
3 
Formel zu gewinnen, braucht man noch eine solche Formel für die Neumannsche Funktion. 
Es ist® 
N). 
) ; Dn / In 
NO=)Y-;, lin -5)AQ+esic-5)40]|: - - - -.- &D. 
:O=V ml B)AQres(li-75) 40] 
ın Setzt man diese Formel und (4,3) in (5,4) ein, faßt die trigonometrischen Funktionen mittels 
re nn iR 2 
der Additionstheoreme zusammen, und beachtet ferner, daß / (3) 3)”, 3 ist, so er- 
(& 3) ! 
hält man 
>). { \ 
7 9 a 7a 
» if #4 Ysnılr > (F nal? + r 2) 
S,.! = — R = l cos -— EN —B. 1) - . 
8,1(8) vs] 7, sin gr) Pit )+ (: 2) 0 0:6) (6,2) 
3 Dieser Wert muß multipliziert mit Y£ für imaginäres [=in reell werden, wie aus (5,5) zu 
>). ersehen ist. Das ist der Fall, wenn man für Yin den Wert —iyn nimmt. Beschränkt man 





sich auch hier auf die beiden ersten Glieder der Hankelschen Reihen, so wird 


er 
: Bin in in, 
dw, a2, DE 1 
de ya! n] Se | ten) 63) 
(6,0). 
48. BETT NG üb ame. 1-5 | 
=) 2 i ya on) 
4 
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® Tor € 

Für die Integrale hat man, da das Integral der asymptotischen Näherung unter sehr all- 


gemeinen Bedingungen, die hier erfüllt sind, eine asymptotische Näherung des Integrals ist, 


3 
/ 3a 
wenn man etwa von „=0 ausgeht und dp =] > Vnd setzt, 


N 
. 


aı/?2 1 i 
w,(y)= ' s| 3a \s Di 6 \dn ED en A 
0) 


] 


und entsprechende Formeln für w, und w,. 


Ausführlich geschrieben lautet jetzt die Knickdeterminante für sehr große n 


ar/3t{ N TarHN N BE 7,05 
N ;| TAEAULL 6) | TArZOLL 6) +] alt) 
0 u 1) 
a, a aaylı & BR > alı- 1 2/7 sodlır 0 =V (65% 
‚ ' 374 67 v9 67 VISNn 67 
dw, (2) dw, (2) dw, (2) 
| dz , = zu | dz z=2u | d; h- zu ’ 
w, (Zu) w,(2,) w; (Zu) | 


Entwickelt man nach den Gliedern der letzten Spalte, erhält man zwei dreireihige Deter- 
minanten, deren erste man wieder nach den Gliedern der ersten und zweiten Reihe in vier 
Determinanten aufspalten kann. Setzt man dabei der kürzeren Schreibweise wegen 


I. “1 
G)=\g@)dn, ODE IDER A re 


und nimmt aus den Determinanten die gemeinsamen Faktoren heraus, so erhält man, wenn 
8/7 
5) 


man noch durch V 32 dividiert: 
« 





+ N +7 rt 
v3 W v3 x 
Gy yngm) + n | +1 +6, 7 ln) i ad +1 +1 N 
v3 N v3 k 
dw, dw, dm, dw, dw, dm, ; 
dz dz dz dz dz dz 3 
2 a7 + + 
v3 + + v3 ° & 
Y ” 7 ‚ ı ) IT (6,7). ; 
+, vn gm) 778 I +1) +6 aim v3 + +1 
dw, dw, de, dw, dw, dw, 
dz dz dz dz dz dz 
a I \ ) L\ 
AL = ae ne 1+45,) 
/Ba 
a vnaln) dw, dw, dw,  =0 
= dz dz dz 
w, w, w, 


Die erste und vierte Determinante in (6,7) werden Null, weil die erste Reihe darin gleich der 
8 / 





zweiten ist. Dividiert man nun durch 6, (n) 2 g(n), so tritt als Faktor der zweiten und 


der dritten Determinante, die sich ja nur durch das Vorzeichen unterscheiden, 
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-all- j F _? G,(n) 
N 3 ik - 
; Ist, | 1) @,(n) 
3 a 7) 


auf und als Faktor der letzten Der erste Faktor wird für sehr große n von der 


2 G,(n)' 


1 A ı 
Größenordnung es der zweite dagegen von der Ordnung v '. Hätte man in der asympto- 


’ 
6,4) As 
tischen Näherung noch weitere Terme genommen, so würden diese mit wachsendem n eben- 
falls stärker gegen Null gehen, als der Faktor des ersten Gliedes. Dividiert man nun noch 
durch F(n) und läßt „> gehen, so ergibt sich: Es muß für unendlich lange Gestänge 


a + +1 

/3 a 

v3 

1 +1 =2 er . (6,8) 
dw,(z.) dw,(z,) 

dw, dw, dw, dz dz 

-0 (6,58 dz dz dz 


werden. Das ergibt unter Benutzung der Reihen (3,4) und (5,5) 2, =3,09, also P=309 YE I (u q). 
. 8 . o . E r . ” ) ’ J 
Dieser Wert ist unabhängig davon, ob das Gestänge im Unendlichen eingespannt oder dreh- 
x RER 1 u x gesp 
bar aufgehängt ist, wie man leicht unter Benutzung der Formeln 


1 1 
(3’ Z,(2))=e”’ 2,_,(@); - 





ter- 
ier 


(.""Zz,@)=-s"Z,.,(e), 


dz dz 


die für alle Besselschen Funktionen Z, gelten, erkennt®). Es ist also damit die zur !-Achse 
parallele Asymptote der Kurven V und IV gefunden. 


.6) ; Setzt man in (6,8) statt der ersten Ableitungen von w die zweite ein, bekommt man 
die Knickkraft des unendlich langen, unten drehbar gelagerten Gestänges. Es ergibt sich 
z„=1,88, also P=188 EI (ug). Dieser Wert bestimmt die zur l-Achse parallele Asymptote 


u der beiden Kurven III und II des Bildes 1. 


TEEN UETATLET KEIN, 


IV. Einfluß des hydrostatischen Druckes auf die Knickkraft. 


7. Mit Hilfe der Kurven unseres Bildes 1 läßt sich z. B. die Frage nach dem Einfluß des 
hydrostatischen Druckes auf die Knickkraft beantworten. Die Folge des Eintauchens des 
Gestänges in eine Flüssigkeit ist die Verminderung des Gewichtes der Längeneinheit von ug 


auf ug. Infolgedessen erhält man zu einem gegebenen ! einen Wert !,, der um den Faktor 


RE > 


be 





LTE LTE 


E ug . 5 Ei " E ie 
j | Fr kleiner ist als das ! für das nicht eintauchende Gestänge. 
a | Zu dem kleineren I, gehört ein P,, das größer ist als das zu / 
| gehörende P. Berechnet man nun P aus P,, so ist dieses mit einem 

’/ (ug\® 
Faktor zu multiplizieren, der um | (u) kleiner ist, als der mit 


dem P zu multiplizieren ist. Man würde in beiden Fällen zu dem 


gleichen Wert P kommen, wenn die Punkte P, und P auf einer 
Kurve 





a2) 





% 


Bild 2. Zur Änderung der Knick- P. u Re = u El Er ee 2er 0 7 A 
kraft durch hydrostatischen ” 
Druck. 


lägen. Diese Kurven verlaufen aber überall wesentlich steiler-als 
unsere Knickkraftkurven. Infolgedessen schneidet die durch P, gehende Kurve auf der 
Parallelen zur P-Achse durch T ein O<P ab. Da also 


P,-E=0P<P.P 


ist, führt nach Gl. (7,1) P, zu einer kleineren Knickkraft als P. Durch den hydrostatischen 
1d Druck wird also die Knickkraft vermindert. 183 
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Statistische Schätzungen auf kombinatorischer Grundlage. 
Von Hermann von Schelling in Berlin. 


Durch ein einheitliches, rein kombinatorisches Verfahren gelingt es, die 
Sicherheit des arithmetischen Mittels von n Beobachtungen einer stetigen 
Größe zu beurteilen und die Differenz zweier solcher Durchschnitie zu be- 
qutachten. Durch Spezialisierung ist es möglich, die Zuverlässigkeit der 
beobachteten Häufigkeit eines alternativen Merkmals zu schätzen und den 
Unterschied zweier solcher Häufigkeiten zu untersuchen. Die bisher üblichen 
Formeln erfahren teils eine wichtige Verallgemeinerung, teils aber auch 
eine wesentliche Berichtigung. 


I. Die Genauigkeit des Mittelwertes bei Stichproben aus einer endlichen Gesamtheit. 
A" einer Gesamtheit von N Gliedern greifen wir n Elemente heraus und messen an jedem 
eine stetig veränderliche Größe x, so daß wir die Beobachtungen «,, &,, ..., &,„ erhalten. 


Aus ihnen bilden wir das arithmetische Mittel M = „eutrat+ + &%,„) und fragen, wie 


a ' ’ te i 
stark es höchstens von dem arithmetischen Mittel «= wett +tay) abweichen kann, 


das wir fänden, wenn wir die Messungen auf alle Glieder der Gesamtheit erstreckten. Aus 


der N-gliedrigen Gesamtheit #,, &,, ..., &y denken wir uns alle (1) Stichproben von je 


n Elementen entnommen und für jede das arithhmetische Mittel M berechnet. So entstehen 


I\ 


(„)W erte für die statistische Größe M. Es ist möglich, die ersten drei Momente dieser 
M-Verteilung zu bestimmen. Zunächst haben wir: 

















| | f \ | 
EiM)=— - -—-- si 
n BEN. „cn 
1 s—ı 
in N n y | 
ee: | 
E(M°’)= n? N ’ \ u T N 3 \ Ki (1). 
ı) „iger (3) ’ 
N rn I RE: 
E(MI=; N \ a ‚N 3 \ ar 2+ N 6 \ ER 7 
= 2|. j ’ ‚3. k 
(1) 2) 9 3) 9 | 
Die Einführung folgender Potenzsummen liegt nahe: 
N 
S= Il (2). 
= 
Dann gelten nach bekannten Sätzen über symmetrische Funktionen die Beziehungen: 
N N 
2 2 uuj=S’—S,; 3 N ar; =3(8, 85, —S); 
1,7 3 
) ) , 
j (3). 
N 0; =28,—38,8,+5° 
ij. 
1 
Nach kurzen Umformungen finden wir: 
1 z 
E(M= N 3, =2 
E(M’)= (N—-n)S,+n—1)S,?] (4) 
e nn N(N— 1)" ’ ’‚ ; 
1 - 
NY - Fur N N—2n)S,+3(N- Ss, oo — 2) 5,’ 
EM) NIN DA alt n)(N n)S, +3(N—n)(n—1)8,8,+(n — lin )S,®] 
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In üblicher Weise erhalten wir die Momente um E(M)=x: 


N—n 1 (S,\? N—n u,(xy) 
’ 2 a kn er = [Bla 
E(M—-a?=o*(M) n(N—1)|IN > 5 N—1 n | 





1 N—n N—2n 4 


1 
? ne . i 
De 


vS-5S+l) 


(). 
_N-n N—-2n u,(zy) | 
"nt mi n? 








Hierbei bedeuten «,(zy) und u,(&y) die unbekannten zweiten und dritten Momente um 
innerhalb der N-gliedrigen Gesamtheit. Eine ähnlich einfache Formel für das vierte Moment 
gilt leider nicht mehr, weshalb auf diese Größe nicht eingegangen werden soll. Dagegen muß 
der Begriff der Schiefe 0 (ey) = u, (ex): [u, (ax ]: noch herangezogen werden. Der angeschriebene 
Ausdruck gilt für die Gesamtheit. Für die Schiefe o(M) der M-Verteilung finden wir nach (5) 
den Wert 

N—2n /N—-1 o(&y) 


N r 
N 2 N . , R ’ . . . » a E (da). 


ER n yn 


In den bisherigen Formeln für die Momente von M kommen noch die unbekannten 
Momente zweiter und dritter Ordnung 1, (zy) und u,(xy) der XN-gliedrigen Gesamtheit vor. 
Sie müssen eliminiert werden, was allerdings nur durch optimale Schätzungen möglich ist. 


’ en is N\. ’ j \ 
Wenn wir für alle »n-gliedrigen Stichproben, deren Anzahl 3 ist, die zweiten und dritten 


Momente bilden, so erhalten wir aus diesen Stichprobenverteilungen folgende strenge Er- 
wartungswerte: 


n—1 N 
El[u,(#,)] = er (ax) 
Aare 
a .—1 »—2 N N 
Elu,(#,)]= er, er 5 Ms len) 
Daraus entnehmen wir die optimalen Schätzungen: 
n N—1 
u, (ey) = a a A (x) 
(6a). 
n n N—-ı1ıN—2 
u,lay)= TEE (2,) 


Die erste Zeile von (6a) liefert für sehr große N die bekannte Behauptung, daß die „wahre“ 


Streuung gleich dem -fachen der beobachteten Streuung ist. Wir haben sie hier 


n 
n—1 
durch unser kombinatorisches Verfahren bewiesen, ohne Annahmen über die Verteilung in 
der Gesamtheit zu machen. Allerdings gilt die Beziehung nicht streng, sondern sie stellt nur 
eine optimale Schätzung dar, was vielleicht manchmal vergessen wird. Unser Ausdruck ist 


ken. 
NV bewirkt, daß für die un- 
eigentliche Stichprobe »=N sich, wie es sein muß, die Identität u, (ey) = u, (ey) ergibt. 
j N—-1N-2,. i 
Entsprechendes gilt für den Faktor N .- yon der Formel für das dritte Moment, für 


etwas genauer als der übliche, da der hinzutretende Faktor - 


die mir übrigens kein Beleg aus der Literatur zur Hand ist. 
Nach diesen Vorbereitungen ist es nun endlich möglich, optimale Schätzungen für Er- 
1 ar i ’ 
wartungswert, Streuung und Schiefe von M = = (2, +2,+:::+%,) bei einer Entnahme dieser 
n-gliedrigen Stichprobe aus einer Gesamtheit von N Elementen anzugeben. Man findet: 
E(M)= x 
0 (M) — Non ‚s’(an) 
3 N n—l (7 
ee a an a ‘), 
N—2n N—-2n n—-1 ol) 
N N—n n-2 yn—2 


o(M) = 
wobei s? (x,) = n > (x; MM)” die beobachtete Streuung und o (x„) die beobachtete Schiefe 


i=1 r 
der Stichprobe sind. Zunächst betrachten wir die Schiefe. Sie verschwindet, wenn n=7; ist. 
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Daraus wird ersichtlich, 


. . o(X 
Für große N und nicht allzu kleine » ist genähert o (M)= % 2) 


ı—2' 
daß die Verteilung von M auf jeden Fall wesentlich symmetrischer ist als die beobachtete 
Stichprobe. Diese Bemerkung ist wichtig für die Beantwortung der naheliegenden Frage, ob 
man berechtigt ist, an M einfach positiv und negativ das k-fache des mittleren Fehlers von M 
anzubringen. Im allgemeinen wird die geringe Schiefe ganz zu vernachlässigen sein, so daß 
die Methode gestattet ist. Da das vierte Moment nicht berechnet ist und je nach der Gesamt- 
heit verschieden ausfallen würde, führt das beiderseitige Anbringen des dreifachen mittleren 
Fehlers nicht streng zu der Restwahrscheinlichkeit von 0.27°/o, die man von der normalen 
Verteilung her gewohnt ist. Jedoch pflegt erfahrungsgemäß die Spanne von 60 bei gering 
asymmetrischen Verteilungen stets einen sehr hohen Prozentsatz zu umfassen, so daß man 
mit außerordentlicher, wenn auch numerisch nicht genau definierbarer Wahrscheinlichkeit 
sagen kann, es müsse A 
M-B30(MSr <M+30o(M) 
oder ausführlicher 


n 


N 1 = al/ N 
0-3] . sn) ZSe=2M+3 ar 


-s(x2,.) el Fe DIR - A 
n—1 1 Pr 


gelten. 

Damit sind die Muturigsgrenzen für das arithmetische Mittel x der N-gliedrigen Gesamt- 
heit gefunden, aus der n Elemente vermessen worden sind. Für die uneigentliche Stichprobe, 
also für a= N ergibt sich streng richtig M=x. Wenn N über alle Grenzen wächst, so er- 
hält man die im Schrifttum übliche Beziehung. 

Für den Praktiker ist unsere Formel aber viel wert. Bei einer anthropometrischen Aufnahme 
bestehe z. B. die Absicht, an allen Männern zwischen 20 und 59 Jahren, die in einem Dorfe ansässig 
sind, die Körpergröße zu messen und dann ihren Durchschnitt 7 zu bestimmen. Einige junge Leute 
sind aber zum Wehrdienst eingezogen und nicht erreichbar. Es ist somit im Augenblick nur möglich, 
den Durchschnitt M für die anwesenden Männer zu ermitteln. Die Gl. (8) setzt uns in den Stand, 
Mutungsgrenzen für den unbekannten 7-Wert aufzustellen, die wesentlich enger ausfallen, als 


wenn wir uns der üblichen Formel bedienten. Nach dieser hätte die Stichprobe W=n- v_, statt 


n Glieder enthalten müssen, wenn die gleiche geringe Unsicherheit herauskommen sollte. 

Wird nicht eine stetige Größe « gemessen, sondern geprüft, ob ein bestimmtes alter- 
natives Merkmal vorhanden ist oder nicht, so können wir der Größe & im ersten Fall den 
Wert +1, im zweiten Fall den Wert 0 zuschreiben. Dann ergeben die Ausdrücke M und x 
die Anzalıl der Merkmalsträger in der Stichprobe und in der Gesamtheit. Die Besonderheit 
des Falles besteht nun darin, daß sich die Streuung u, (xy) in der Gesamtheit durch den zu- 
gehörigen Erwartungswert x streng ausdrücken läßt. Es ist nämlich u, (@y)=#(1— x). So- 
mit folgt nach (5) aus 

M+ko(M)=x, 


N—n z(l—%) 
N—1 n 


|| 


@—NY=l*: 


1 ’ 
- ‚so haben wir: 
Non 
z(1—x) 


G-M=k. 


Benutzen wir wieder die Abkürzung n’—n-- 


Die Auflösung der quadratischen Gleichung liefert: 


.2 „2 
8: vB Mn+y&k] n’M( +) TER 
Wächst N über alle Grenzen, so geht n’ nach n. Dann ergibt die Gl. (9) die von R. Prigge [3]') 
aufgestellten Mutungsgrenzen für den wahren Anteil der Merkmalsträger. Ich habe diese 
Grenzen für k=1,2,3 in ausreichendem Maße tabuliert [4b]. Es ist beachtlich und wird hier 
zum ersten Male erwähnt, daß diese Tafeln auch dann angewendet werden können, wenn die 
n-gliedrige Stichprobe nicht einem unendlichen Kollektiv, sondern einer Gesamtheit von 


N Elementen entstammt. Man läßt den beobachteten Prozentsatz M =":=p, der Tafeln 


N—1 


der Merkmalsträger unverändert, benutzt aber statt n» die Größe n"’—=n: u. 


1) Die Zahlen in eckigen Klammern verweisen auf das Schrifttumverzeichnis am Schluß der Arbeit. 
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Als Beispiel betrachten wir folgende Aufgabe. Von einem Stammelternpaar sind bisher 31 Kinder 
und Kindeskinder entsprossen. Von ihnen konnten 25 daraufhin untersucht werden, ob sie ein be- 
stimmtes Merkmal aufweisen oder nieht. Man fand 10 Merkmalsträger. Für die restlichen 6 bereits ver- 
storbenen Familienmitglieder war eine Entscheidung nicht mehr möglich. Die Mutungsgrenzen der 
wahren Häufigkeit 7 des Merkmals innerhalb dieser Sippe sind gesucht. 


10 
Es ist N=3, n=8, M=2=40%. 
N—1 Ss 
Wir finden N"—=n-— 125. 
N—n 
Für =: = Very dar Tee), 


lesen wir aus der Tabelle ab: 

30.0 u, <E< 533 %,. 
Wenn man den Umstand nicht berücksichtigt, daß die Sippe 31 Köpfe zählt, so entnimmt man aus der 
Tafel für k=3, n=25, n, = 10 die Grenzen 

23.14 9, < T< 59.63 0/0 
die wesentlich weiter sind. Die auf Grund unserer kombinatorischen Methode erzielte Verfeinerung 
kann sich also sehr nützlich auswirken, wenn es sich um Stichproben aus klar definierten endlichen 
Gesamtheiten handelt. 


II. Verfahren zur Beurteilung der Differenz zweier solcher Mittelwerte. 


1. Methode für den Fall, daß die Lage der Differenz unbekannt ist. Soll bei Messungen 
einer stetig veränderlichen Größe das arithmetische Mittel einer Stichprobe von n, Gliedern 
aus einer Gesamtheit von N, Elementen mit dem arithmetischen Mittel einer n,-gliedrigen 
Stichprobe aus einem Kollektiv des Umfangs N, verglichen werden, so liefert auch das 
kombinatorische Verfahren in aller Strenge die: bekannte Formel 0° (M, — M,) = 0° (M,) + 0° (M,). 
Allerdings müssen «ie auf der rechten Seite der Gleichung stehenden Streuungen nach 
Formel (7) und nicht nach den nur für sehr große N,, N, gültigen üblichen Ausdrücken 
berechnet werden. 

Ergiebiger wird diese Betrachtung, wenn die Messungsergebnisse in den beiden Reihen 
nur aus den Werten +1 und O bestehen, wenn wir es also mit einem alternativen Merkmal 
zu tun haben. Statt auf die Gl. (7) können wir dann auf die exakte Formel (5) zurückgehen. 
Wir schreiben also: 

N, m, ; My (En u N, N, Ma(YN2) 
N —1 n X,—1l u: 


1 


0° (d) = 

ze 
wobei wir die Messungen in der zweiten Reihe mit y bezeichnen. Nun gilt, wie schon er- 
wähnt worden ist, für Alternativen streng u, (ay,)=#(1— #), u, (yx.) = y(1— y). Die Gleichung 
dy.on =(#- M + ko(d) lautet also ee 





e n, z(l-r) N,—n, yi—p) 
[dneon. (A W}: N, 1 n, a l a 1 ; N, | : 
. . . .. ’ N, r 1 ’ N, . 2 
Mit unseren üblichen Abkürzungen n, =n;- Vonrmanıy haben wir einfacher 
1 er. ? 
u. 
Idyeon.  (eF=h* | u. r N, | 


Ausmultipliziert ergäbe sich eine allgemeine quadratische Gleichung in x und 4, welche bei 
geometrischer Deutung in einem rechtwinkligen Koordinatensystem eine Ellipse in nicht aus- 
gezeichneter Lage darstellt. Jeder Punkt (x, y) innerhalb dieser Mutungsellipse liefert ein 
Wertepaar &, y, das mit einer gewissen noch zulässigen Wahrscheinlichkeit den wirklichen 
Werten entsprechen könnte. Uns interessieren bei dieser Betrachtung aber nicht die Werte 
= und 9 einzeln, sondern nur die Differenz d=r—y. Die Gleichung x y=D mit un- 
bestimmtem D stellt geometrisch eine Schar paralleler Geraden dar, die mit der x-Achse den 
Winkel 45° einschließen. Zwei Gerade dieser Schar, 


z-y=D wd x—-y=D" 
berühren die Mutungsellipse. Es gilt dann für die Differenz d=x— y der Koordinaten eines 
jeden inneren Punktes der Ellipse die Berater 
D'<dsD". 
Somit sind D’ und D’ als die Mutungsgrenzen der echten Differenz d=xr — y erkamnt. 
Die tatsächliche Berechnung dieser Grenzen ist etwas langwierig aber elementar. Ich 
will hier darum nur das Resultat angeben. Man hat die Geradenschar mit der Ellipse zum 
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Schnitt zu bringen und darauf zu achten, für welche Geraden die beiden Schnittpunkte zu- 
sammenfallen. Das sind die gesuchten Tangenten. Für sie gilt 


D' | Eu 1 | In “ 
D’I (n,+n,))+k? | 7 Rz ) Apeon. 


(10). 
n '£n, ‚ ‚ D ‚np P} D D l 
+k % ia >, In, +n,) —tn,'n, d’yeon. kn’ +m, I; 
ı "a ; 
Eine wesentliche Vereinfachung tritt für n,’=n,'=n'’ ein. Dann hat man 
a PER. / kyan (a) + 10: 
pi In pen Ayeon. FkYy2n( 8 rn 5 
Wenn über die Größe der Gesamtheiten nichts bekannt ist, wird man N, und N, als unend- 
lieh voraussetzen. In diesem Falle ist n,’=n,, n,=n, 


Als Daten sind außer dem Umfang der beiden Reihen nur die Differenz d,.op. = Xyeon. 
Yeon. Selbst, nicht auch con. Und Yyeon. einzeln benutzt worden. Diese Größen pflegt man 
meistens zu kennen. Ein Verfahren, das nicht alle Informationen ausnutzt, hat stets einen 
verminderten Wirkungsgrad. Die durch die Gl. (10) bestimmten Grenzen sind daher in der 
Regel zu weit. 

Hin und wieder ergibt sich aber eine Gelegenheit, bei der man auf die Formel (10) wirklich an- 
gewiesen ist. Eine Firma berichte etwa, sie habe im letzten Jahre nur 150 Stück einer bestimmten 
Maschine abgesetzt, während es im Vorjahre 180 Stück gewesen seien. Da aber der Anteil der höher- 
wertigen Ausführung um 23°/, gestiegen sei, habe der Umsatzwert keinen Rückgang erfahren. Der Zug 
zur Qualität werde vermutlich anhalten. Mit n’=n,=150, n’=n,—=180, dyen.—=+ 0.23 und k=3 
erhält man nach (10) DD’ =+6.45%,, D’"=-+38.33°/.. Auch die untere Grenze ist positiv. Der Zug 
zur Qualität beruht also nicht auf Zufall, das kann man trotz der absichtlich verschleiernden Dar- 
stellung des Berichtes auf diese Weise nachprüfen. 

2. Durchführung eines Ansatzes vonR. A. Fisher. Zum Vergleich zweier arithmetischer Mittel 
gibt es noch ein anderes Verfahren, bei dem alle Daten voll ausgenutzt werden. In der Praxis wird sehr 
oft die Frage gestellt, obzwei Reihen #,,%,, .... 2m und Y,, Ya... .,Y%„ als Stichproben aus der gleichen 
Gesamtheit angesehen werden können oder nicht. Lautet die Antwort bejahend, so sind die Mittel- 


m 


1 ” 1 DE t , } i 
werte M.= — \ x; und M,„= — \ y; Zufallsvarianten des arithmetischen Durchschnitts 
x i y ij 
; Do ; N " 


im! j=1 

dieser Gesamtheit. Diese Überlegung hat R. A. Fisher [1,2] zudem Gedanken angeregt, die beiden 
Reihen zu einer einzigen Reihe 2,,2,,..., 2y zusammenzuwerfen, wobei N=m+n ist. Dann 
denkt er sich die Reihe auf alle N!:(m!n!) Arten in zwei Abschnitte der Längen m und n 
zerlegt und für diese jeweils die Mittel M, und M, berechnet. Er stellt sich vor, man zähle 
ab, wie oft die Differenz d= M, — M, extremer auffällt als bei dem empirischen Reihenpaar. 
Auf diese Art läßt sich beurteilen, ob die gegebenen Reihen als zufällige Stichproben aus 
einer einzigen Gesamtheit aufgefaßt werden können oder nicht. 

Wirklich durchführbar ist dieses Verfahren natürlich nur für ganz kurze Reihen. 
Fisher scheint aber übersehen zu haben, daß sich die ersten drei Momente von d einfach 
bilden lassen. So ist also ein Vergleich zwischen d,..,. und dem mittleren Fehler o (d) mög- 
lich. Ich habe diese bisher offenbar nicht bekannte Tatsache schon kurz angezeigt [4d] und 
will hier den Beweis folgen lassen. Man findet zunächst: 





m /n N 
E(d) = 2 (3) 1) >; 
5 mn a 
le) le N 
sr, lu > HN me tn mul. 3 -. 
et - 
(m) i | . (3) (3) ! (11). 
Be); u 1 : Na: 4 > Si. mn mn), 3 PRF 
e (F) > — . "sn m’ n? m’n mn I Ko“ 
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Wir führen wieder die Potenzsummen 
N 
v y r 
ee VE 
kl 
ein und erinnern an 
N N 
Y y ‘ Y it v ‘ 
> zk 2: Ss: — 8, s> &k zı=3(S, 8, S,) 
kl k,l 
1 1 
(3). 
\ 
” v oc 5 y Y ' v3 
6N22=28,—38,8,+S8 
nur 
1 





Damit erhalten wir nach einigen Umformungen: 





E(d) =0 
nJ 2 N N (S; S,? 
EAN ma|ln N) (12) 
.. N N (w—m)-N(S, „S, 8, „S, 
EA)=y_i na wm N Sn Nt2N) 
’ . ’ . a 1 
Bezeichnen wir das zweite und dritte Moment der z-Reihe um M. = ywatzt tz) 
mit u,(2) und u, (2), so ist 
E(d\ =V 
i N N 
BI=H—Z Te 
ua N N (n-m):-N (l2a) 
E (d’) = N-iNS mm Mid 
EN Jori... yv\—1 n 
e (d) "yVmn N-—2 e 12) 


Bemerkenswert ist vor allem, daß die Schiefe o(d), die für n=m und für o (2)=0 
verschwindet, selbst bei großem o(z) praktisch vernachlässigt werden kann, da beide Fak- 
toren von o(z) stark verkleinern. Die Verteilung von d ist also nahezu symmetrisch, so 
daß der Ansatz 

N N 
d+ko(=d+ | V 


; 11, (2 
lmn'?’ 


erlaubt ist. 
Um den Anschluß an die übliche Bezeichnung zu erhalten, schreibe ich das zweite 
Moment in der Form: 


hiBal „ N ms. —+nsy d? ® 
E (d?) = 0°? (d) = en ie =; + Nr: (13), 
m n 
R 1 ” a a { 
wobei 8,’ = a M,” und s/)’= = \ (yj - M,„)” ist. Dieses Ergebnis weicht von 
— 
i=1 je! 


der bisher üblichen Formel 


N msf.+tns,’ 


"(d=Y 9 (14) 


mn 
entschieden ab. Unwichtig ist zwar der Unterschied um eine Einheit im ersten Nenner, be- 
deutsam aber mein Zusatzglied + d’:(N — 1). Vor allem beim Vergleich kurzer Reihen macht 
sich sein Einfluß stark bemerkbar. 

Im Sonderfall der Doppelalternative sei in der ersten Reihe von m Gliedern das Merk- 
mal @„, mal vorhanden, a@„, mal nicht vorhanden, in der zweiten n - gliedrigen Reihe werde 
es a„, mal beobachtet, a„, mal nicht beobachtet. Wir haben also die 2X 2 Tafel 
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Amı Ama| M 
Aynı One n 
A, 4A|N. 
Für x und y setzen wir bei positiver Beobachtung genau wie früher den Wert +1, 
bei negativer Beobachtung den Wert 0 an. Dann ist: 





a, a Am An 
M, a2 ER M, Br d- - M; M,= 1 Y a 
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ae Imı | Amı ), FRI... 1 =), 
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Die Formel (13) liefert dann: 
0: (d) - (Nms’+N ns,’ +tmnd) 
mn(N— 1) “ Y 
Ei, Bl a te ' 
n Na Pas, yna = ma N Ay mo. zmna a 
mniN 1) mi nı mn m ı nı mi ni mı Anı 
l ! 
- (NA 1,?) 
mn(N—1) z ’ 
1 4,4, 
N-1l mn 


Bis auf die Bezeichnung (ich habe sonst m =a,, n=a, geschrieben) ist das die von mir 
zum ersten Male [4c] klar herausgearbeitete strenge Formel für die Streuung der Differenz 
zweier beobachteter Häufigkeiten, die in etwas anderer Weise von mir bereits früher [4a] 
abgeleitet worden ist. Wenn der Ansatz von R. A. Fisher überhaupt einer Rechtfertigung 
bedarf — dieser Ansicht ist E. S. Pearson [2] —, so wird diese dadurch geleistet, daß die 
allgemeine Formel für den Spezialfall der Doppelalternative zu dem richtigen Ausdruck führt, 
was bei der im Schrifttum üblichen Gl. (14) nicht der Fall ist. Die Angriffe E. S. Pearson’s 
gegen die kombinatorische Methode halte ich für unberechtigt. Das Verfahren zeigt sich bei 
vielen Aufgaben sehr wirksam. Ein schönes Beispiel liefert die Berechnung der Streuung 
der Lexisschen Zahl durch Fr. A. Willers [5]. 


Nachtrag bei der Korrektur: 


Inzwischen habe ich erfahren, daß E. I. G. Pitman in einer sehr lesenswerten Arbeit 
«Significance tests which may be applied to samples from any populations> (Suppl. I. R. statist. 
Soe., London, Bd. 4 (1937), S. 119 bis 130) sich ebenfalls in positivem Sinne mit dem Ansatz 
von R. A. Fisher beschäftigt hat, welcher die Beurteilung der Differenz der arithmetischen 
Mittel zweier stetigen statistischen Größen zum Ziel hat. Zu meiner einfachen Formel (13) 
gelangt er aber nicht. Dagegen geht aus seiner Arbeit hervor, daß OÖ. N. Anderson in 
seinem Buch «Einführung in die mathematische Statistik» (Berlin 1935, insbesondere S. 219) 
bereits Ausdrücke entwickelt hat, die den Formeln (7) dieser Arbeit gleichwertig sind. 211 


Schrifttum. 


[II R. A. Fisher: J. R. anthrop. Inst. Bd. 64 (1936), S. 57; zitiert nach 

E. S. Pearson: Some aspects of the problem of randomization. Biometrika, Cambridge, Bd. 29 

(1937), S. 53 bis 64. 

3| R. Prigge: Fehlerrechnung bei biologischen Messungen. Naturwiss. Bd. 25 (1937), S. 169 bis 170. 

4a) H. von Schelling: Die Beurteilung der Treffsicherheit von Prognosen. Deutsch. Statist. 
Zentralbl. Bd. 27 (1935), S. 141 bis 148. 

(4b! H. von Schelling: Die mathematisch -statistische Bewertung von Stichproben und deren Be- 
deutung für die Beurteilung von Tierversuchen. Zweite Mitteilung. Die Abgrenzung von Mutungs- 
bereichen. Arb. a. d. Staatsinst. f. experim. Therapie z. Frankfurt a. M. Heft 37 (1939), S. 28 
bis 54. 

(4e] H. von Schelling: Über die exakte Behandlung des Zusammenhanges zwischen biologischen 


Merkmalsreihen. Arb. a. d. Staatsinst. f. experim. Therapie z. Frankfurt a. M. Heft 39 (1940). 


S. 35 bis 71. 


[4d] H. von Sehelling: Fehlerrechnung bei biologischen Messungen. Naturwiss. Bd. 28 (1940), 


S. 430 bis 431. 
5] Fr. A. Willers: Prüfung der Einheitlichkeit eines Massenfabrikates. Z. angew. Math. Mech. 


Bd. 14 (1934), S. 77 bis 84. 
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BUCHBESPRECHUNGEN 


FirstReportonViscosityandPlasti- 
city {prepared by the Committee for the Studs 
of Viseosity of the Academy of Sciences at Amster- 
dam). 1. Abschnitt. Teil XV, No. 3. 2. Aufl. VIIL + 
273 5. Amsterdam 1939, N. V. Noord-Hollandsche 
Uitgevers Maatschappij. Preis geb. 12,— f. 

Dieser „erste Bericht“ enthält eine Reihe von 
Darstellungen der wesentlichen Probleme auf dem 
(rebiet der Viskosität und Plastizität und ihrer Be- 
deutung. Es beginnt mit einem Kapitel von J. M. 
Burgers (Delft). in welehem nach den üblichen 
Detinitionen und einführenden Erläuterungen all- 
gemeine Betrachtungen über den Mechanismus 
der inneren Reibung zusammengestellt werden. 
Neben dem Maxwellschen Relaxationsansatz wird 
die Vorstellung der Aufteilung eines Stoffes in 
elastische und flüssige Glieder mit Zähigkeit und 
Fließgrenze erörtert. wobei die Behandlung von 
Modellen auch mit komplizierten Parallel- und 
Reihenschaltungen einen großen Raum einnimmt. 
Ausführlich wird auch das Pran.dtlIsche Modell des 
plastischen Verhaltens erörtert. Es folgen Betrach- 
tungen über die Strukturelemente und die Viskosi- 
tät der Flüssigkeiten. Das Kapitel schließt mit 
der Theorie von Andrade über die Berechnung der 
Zähigkeit am Schmelzpunkt. 

Das zweite Kapitel widmet derselbe Verfasser 
der experimentellen Untersuchung der Fließerschei- 
nungen. Bemerkenswert ist der Oberbegriff „fluid“ 
für einen Stoff, der bei kleinsten Spannungen (also 
ohne Fließgrenze) fließt. Zu ihnen gehören einer- 
seits die „liquids“, für welche ein festes Fließ- 
Er; 2 i 
gi eine Funktion 
der Schubspannung r, also D=f (r)] gilt. insbe- 
sondere also die Newtonschen Flüssigkeiten mit 
dem Gesetz r 7-D: andererseits aber Stoffe. 
deren Verhalten nur durch eine zeitlich veränder- 
liche Funktion dieser Größen darstellbar ist. die 
also z. B. Thixotropie oder Selbstwiederherstellung 
self-restoring) oder sonstwie ein Gedächtnis zei- 
een. Neben dem „fluid“ gibt es das elastische 
Verhalten |—F(y)| und das plastische, bei dem 
zum Eintritt des Fließens erst eine Fließzrenze 
überschritten werden muß. Mit der Erörterung der 
Verwendbarkeit der verschiedenen Meßverfahren 
für die Zähigkeit (Kapillaren, Couette, Ausfluß, 
Plastometer) werden als Beispiele Messungen an 
Chatterton-Compound und an Gummi besprochen. 
wobei auch die allgemeineren Tensorbeziehungen 
kurz angedeutet werden, 

H. G. Bungenberg de Jong (Leiden) berichtet 
über Viskositätsmessung besonders für die Kolloid- 
ehemie. Mit Bezug auf die Notwendigkeit. zur Prü- 
fung der Einsteinschen und Smoluchowski- 
schen Formeln über die Abhängigkeit der Zähigkeit 
von der Konzentration die „relative Viskosität“ recht 
genau zu messen, behandelt er zunächst die syste- 
matischen Fehler der Kapillarviskosimeter beson- 
ders in Abhängigkeit von den Kapillarenabmes- 
sungen und dem Meßbereich, sowie die Bestim- 
mungsmethoden. An hydrophilen Solen werden so- 
dann zwei Typen unterschieden: Iyophobe Solen. 
deren Viskosität sich nur wenig von der des Dis- 
persionsmittels unterscheidet, und lyophile Solen. 
die auch bei kleinen Konzentrationen beträchtlich 
höhere Viskosität aufweisen. Die Gültigkeit des 
Poiseuilleschen Gesetzes als Ausdruck desN ew 
tonschen Charakters der Flüssigkeit kann mittels 
Ostwaldschen Viskosimeters geprüft werden. Bei 
Nichtgültigkeit besitzt das weniger genaue Couette- 
Viskosimeter Vorzüge. "Für die Stoffe ergibt sich 
die Dreiteilung: a) Poiseuille gilt. b) Poise- 
uille gilt nicht, Hysteresis besteht nicht: e) Po- 
iseuille gilt nieht und Hysteresis besteht. 


gesetz |[Fließgeschwindigkeit D 


Für die hydrophilen Solen der ersten Gruppe. 


die dem Poiseuilleschen Gesetz folgen. 
ergibt die viskosimetrische Analyse Abwei 


Einsteinschen Formel. 
für die drei Ursachen als möglich erwähnt 
werden: nichtkugelige Form der dispersen Teile. 
Schwellung der Teile (Solvation) und Porosität 
(Vollsaugen) der Teile. Die Bedeutung der Vis 
kositätsmessungen für die Solvationsforschung 
und für die Erforschung der elektrischen Ladungen 
der dispersen Teile wird an typischen Beispielen 
erörtert. 


chungen von (der 


Es folgt ein Kapitel über Viskosität und Plastizi- 
tät vom technischen Gesichtspunkt von U. J. van 
Nieuwenbureg (Delft). der auf die Unzulängliel 
keit der einen Kenngröße n (Zähiekeitswert im 
Newtonschen Gesetz) bei nicht newtonischen Flüs 
siekeiten hinweist und (die Wiedergabe des ganzen 
Zähigkeitsgesetzes D ir) befürwortet. wie dies 
analoge z. B. auf dem Gebiet der Werkstoffkunde 
längst geschieht. Entsprechend findet eine Grup- 
pierung der Stoffe nach dem Verlauf des Gesetzes 
Anwendung. d. h. danach, ob » e = konstant ist 
oder mit wachsendem r wächst oder fällt und ob 
eine Fließgrenze vorhanden ist (plastisch) oder 
nicht (flüssig). Die vielfältigen technischen Meß 
verfahren finden insbesondere bei Besprechungen 
der Petroleumderivate Erwähnung, getrennt nach 
solchen für niedrige konstante Zähigkeit und sol- 
chen für hohe und veränderliche Zähigkeit bis zur 
Penetrationsmessung und Messung durch Scher- 
biegebeanspruchung von Stäben. Mit den Eigen- 
schaften von Gummi finden nochmals Plastometer 
messungen Erörterung. Weiter werden für Kunst 


harz, Farbtünche. Ton. Glas und andere Stoffe 
typische Besonderheiten hervorgehoben. Einige 
allgemeine Betrachtungen über eine mögliche 


Superposition von Zähigkeit und Plastizität sowie 
Bemerkungen zur Nomenklatur schließen das 
Kapitel. 

Die Plastizität kristalliner Substanzen, insbeson- 
dere der Metalle. behandeln W.G. Burgers(Eind- 
hoven) und J.M.Burgers (Delft). Im ersten Teil 
dieses Kapitels findet man die übliche Behandlung 
der Zerreißkurven von kristallinen Gefügen und 
Bemerkungen über Reibungswirkungen als Ursache 
des Fließwiderstandes, worunter einerseits die 


Zähigkeit (— Viskosität!) im elastischen und im 
plastischen Gebiet, andererseits die Fließgrenze 
verstanden werden kann. Die grundlegenden 


Fließgesetze von Andrade und Prandtl-Deut 
ler werden wiedergegeben, der Hystereseverlust bei 
schwingender Beanspruchung erörtert. Der zweite 
Teil ist der plastischen Verformung des Einzel 
kristalis gewidmet. Ausgehend von den bekannten 
ersten Ergebnissen über Gleitflächen und Gleit- 
richtung werden die Verfestigung der kritischen 
Schubspannung und ihre Abhängigkeit von der 
Temperatur bei Metallen beschrieben und die merk- 
würdigen Erscheinungen an Steinsalz in Zusammen- 
hang mit der Smekalschen Lockerstellenvorstel- 
lung dargestellt. Von besonderem Interesse ist die 
dann folgende Zusammenstellung der Vorstellungen 
über die Möglichkeit plastischer Vorgänge im Kri- 
stallgitter. wobei sich als Grundgedanken heraus- 
schälen die Theorie der ortsfesten Lockerstellen 
(Smekal) die Theorie der wandernden Ver- 
setzuneen (Taylor), die eine Verfestigung lie- 
fert. die Theorie der thermodynamischen Atom- 
sprünge (R. Becker) eine Vereinigung der 
Becker-Smekalschen Gedanken (Orowan). 
die Zeiteffekte (Zähigkeit) liefert, und schließlich 
Anregungen zur Vereinigung der Taylorschen 
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und Orowanschen Theorie, die Verfestigung init einer wertvollen Übersicht über den Zusammen- 


und Zähigkeit zugleich liefert (Burgers). 

Im letzten Kapitel werden Viskositätseffekte in 
lebendem Protoplasma und in Muskeln von H. J. 
Jordan (Utrecht) erörtert. Nach einführenden Be- 
merkungen über die diesbezügliche Fragestellung 
für lebende Stoffe und über die Viskosität und 
Elastizität von Protoplasma bilden vor allem (lie 
elatten Muskeln niederer Organe und gewisser 
wirbelloser Tiere, die solchen Organen gleich zu 
werten sind. den Gegenstand der Ausführungen. 
Die Wiedergabe experimenteller Untersuchungen 
über das Verhalten unter Last zeigt den typischen 
Unterschiel zwischen elastischen und plastischen 
Muskeln. Als allgemeine Merkmale der Plastizität 
beim Streekvorgange werden der Einfluß der Tem- 
peratur und der Hydration auf die Dehnungs-Zeit- 
Kurve, die Streckgrenze und die visko-elastische 
Erholung dargestellt und an einem mechanischen 
Modell erläutert und mathematisch durchgeführt. 
Von den Änderungen des viskosen Widerstandes 
während der Dehnung wird besonders eine Zunahme 
(„Sehneepflugeffekt“) besprochen, die als abhängig 


von der Geschwindigkeit nicht unmittelbar der 
Verfestirung bei Metallen gleichzusetzen ist. Der 


Bedeutung des tonischen Spannungszustandes wird 
durch weitere Ausführungen Rechnung getragen. 
Die Anwendung auf die transversalgestreiften 
Muskeln der Wirbeltiere wird an Hand eines 
elasto-viskosen Modells angedeutet. 

Im ganzen stellt dieser „erste Bericht“ einen 
wertvollen kritischen Überblick über den heutigen 
Stand unserer Kenntnisse und der Forschung auf 
allen mit Fragen der Viskosität und Plastizität zu- 
sammenhängenden Gebieten dar, der überall an die 


wiehtiren Grundfragen heranführt und auf die 
einschlägige Literatur hinweist. Wenn dabei 
naturgemäß auch einige offensichtlich schwache 


Methoden ohne Beanstandung wiedergegeben wur- 
den. so ist doch beherrschend eine gute kritische 
Behandlung festzustellen. die durch ihre weg- 
weisende Ergänzung einen erfreulichen Fortschritt 
darstellt. 


Danzig-Langfuhr. H. Fromm. 247 
Dr. phil. WILHELM SPÄTH, Physik und 
Teehnik der Härte und Weiche. VIl+ 


250 S. m. 214 Textabb. Berlin 1940, Verlag Julius 
Springer. Preis geb. 19,50 M. 

Das Ziel des Buches ist es, als neuen Begriff 
in die Werkstofflehre das Verhältnis einer Span- 
nung zu einer bleibenden Formänderung einzu- 
führen. Beim Zugversuch erscheint so neben dem 
Rlastizitätsmodul ein „Plastizitätsmodul“ als Ver- 
hältnis von Zugspannung zu bleibender Dehnung. 
bei der Härteprüfung eine neue „Härte“ als das 
Verhältnis einer mittleren Eindruckspannung zu 
einer mittleren Eindrucktiefe. Der Kehrwert der 
„Härte“ heißt „Weiche“. 

Diese Neuerung erscheint wenig <lücklich, da 
sie nicht nur die Prüfune umständlicher macht. 
sondern auch begrifflich keinen Vorteil bietet und 


sachlich keine typische Kenngröße liefert. Be- 
sonders aber würde die Bezeichnung .„Härte“ für 
diesen neuen Begriff abzulehnen sein, da der 


Grundbegriff der Härte doch sehr viel sauberer 
festliegt. als dies der Verfasser darstellt. wenn 


auch ihre meßtechnische Erfassung noch nieht 
vollkommen geglückt ist. Leider verwechselt der 
Verfasser wiederholt Grenzwerte mit laufenden 


Werten oder den Funktionswert mit der Ableitung, 
wie z. B. Verfestigung mit Abnahme der Verfesti- 
geungsfähigkeit und ist auch sonst in seinen voran- 
gestellten „Grundbegriffen“ so verwirrend, daß es 
schon einer Sicherheit des Lesers bedarf, wenn er 
keinen Schaden nehmen soll. 

Im übrigen findet man eine Zusammenstellung 
und Beschreibung der üblichen Härteprüfverfahren 


hang der verschiedenen Härtewerte (bisher üblicher 
Art) untereinander und mit anderen Werkstoffeigen- 
schaften, ohne dabei auf erundsätzlich Neues zu 
stoßen. Die Bedeutung der Kalthärtung wird kurz 
gestreift. Sonderprüfungen an unvulkanisiertem 
Kautschuk {Plastometer), vulkanisiertem Kautschuk 


(Rückprallversuch), Kunststoffen und Holz sind 
kurz. behandelt. 
Danzig-Langfuhr. H. Fromm. 246 


Prof. Dr. F. RAUSENBERGER, Theorie der 
Rohrrücklaufgeschütze. 2. Aufl. bearb. v. 
Prof. Dr. ©. von Eberhard. IX + 175 S. m. 64 Abb. 
Berlin 1939. Verlag Julius Springer. Preis geb. 
16.20 M. 

Das vielen Artilleristen und fast allen Geschütz- 
konstrukteuren bekannte Buch ist aus den Vor- 
lesungen des 1926 verstorbenen Lehrers der Waf- 
fenkonstruktionskunde an der ehemaligen militär- 
technischen Akademie zu Charlottenburg. und spä- 
teren waffentechnischen und artilleristischen Vor- 
standsmitgliedes der Krupp A.G. hervorgegangen. 
Die vorliegende 2. Auflage rührt von dem jüngst 
verstorbenen Chefballistiker bei Krupp. dem älte 
sten Schüler von ©. Cranz, Prof. ©, v. Eberhard. 
her. Er hat das Werk entsprechend den Fortschrit 
ten der Geschütztechnik erweitert und vergrößert. 
wobei er auch einige Fachgenossen zur Mitarbeit 
in besonderen Fachgebieten hinzugezogen hat. Das 
gilt insbesondere für die Bestimmung der Brems- 
kräfte und Bremswiderstände und der Erwärmung 
der Rohrbremsen, Untersuchungen, an denen Dip!.- 
Ine. Geore Hayn mitgearbeitet hat. Ein ausführ- 
licher und völlig neu bearbeiteter Abschnitt bezieht 
sich auf die Vorholmittel. insbesondere die Berech- 
nune der Vorholfedern, von ©. v. Eberhard selbst 
herrührend. Um die zur Erzielung eines langen 
Rücklaufes bei großer Rohrerhöhung erforderliche 
Rückverlegung der Schildzapfen des Rohres wie 
der auszugleichen, sind besondere Feder- oder 
Luftausgleicher anzubringen, mit denen sich Dipl.- 
Ine. F. Böminghaus und F. Kemnitz beschäftigen. 
Der eigentlichen Geschütztheorie sind drei Anhänge 
angefügt. von denen der erste innerballistischen 
Inhalts ist und sich mit der Nachwirkung des Gas- 
druckes der Pulvergase beschäftigt. Die dadurch 


bedineten Änderungen der Geschoßgeschwindig- 
keit und der Gasgeschwindiekeit an der Rohr 


mündung werden nach Formeln bestimmt. die von 
Dipl.-Ing. H. Kratz zusammengestellt worden sind. 
Im zweiten Anhang werden kurz Geschütze mit 
Doppelbremse betrachtet. Zuletzt kommen noch ein 


paar Bemerkungen über den Hochsprung eines 
Rohrrücklaufzeschützes hinzu. insbesondere über 
den Einfluß der Bremskraft auf die Größe des 


Hochsprunges. 

So ist das alte klassische Werk wieder auf die 
seeenwärtice Höhe der Geschütztechnik zebracht 
worden und wird seinen Zweck erfüllen, dem jun- 
gen Ingenieur und dem für technische Fragen inter- 
essierten Artillerieoffizier eine Hilfe dabei zu bie- 
ten, wenn er sich über das Gebiet des Geschütz- 
baues unterrichten will. 

Berlin. R. Rothe. 


©. CRANZ und O. von EBERHARD, Die neu 
zeitliche Entwicklung der Schuß- 
waffen. 2. erweiterte Aufl. 72 S. m. 23 Bildern. 


243 


Berlin 1939. VDI-Verlag G.m.b.H. Preis broseh. 
2M. 
Wer die theoretische Ballistik nicht nur vom 


mathematischen Standpunkte aus, sondern auch im 
Hinblick auf die Technik der Gewehre und Ge 
schütze studieren will, ohne sich dabei in fern- 
lierende konstruktive Einzelheiten verlieren zu 
müssen, dem sei dieses Heft aus der Feder zweier 
Fachleute ersten Ranges auf das angelegentlichste 
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empfohlen. Es ist sehr zu beklagen, daß Herr von 
Eberhard, einer der ältesten und besten Schüler des 
Altmeisters Cranz, der Chefballistiker von Krupp. 
vor einigen Wochen seine Kameraden und Fach- 
genossen im Tode hat verlassen müssen. Das aus- 
gezeichnete kleine Werk wird jedem, der sich über 
die Technik der Feuerwaffen in geschichtlicher und 
physikalischer Hinsicht kurz unterrichten will, ein 
ebenso interessant geschriebener wie anregender 
und zuverlässiger Berater sein. Die vorliegende 
„weite Auflage führt die Waffentechnik bis zu den 
neuesten Errungenschaften, soweit diese überhaupt 
einem öffentlichen, nicht militärischen Leserkreise 
ohne Gefährdung der notwendigen Geheimhaltung 
mitgeteilt werden dürfen. Zu den zwei schon vor- 
handen gewesenen Teilen über die Entwicklung der 
Handfeuerwaffen und über die Entwieklung der 
Artillerie ist ein dritter hinzugekommen (Schluß- 
bemerkungen allgemeiner Art), der zunächst über 
die Entwicklung der Theorie des Schießens in den 
letzten 25 Jahren in treffender exakter Weise die 
wesentlichen Fortschritte der äußeren Ballistik be- 
richtet und sodann den theoretischen Wirkungsgrad 
einer Schußwafle ins Auge faßt. Während früher 
das Treffen bei der Artillerie im wesentlichen auf 
ein Einschießen abgestellt war, die Schußweiten 
also a priori nur ungefähr bekannt zu sein brauch- 
ten, ist heute die Möglichkeit des Einschießens nur 
noch in seltenen Fällen zegeben. Vielmehr soll 
jetzt gleich die erste Salve genau in einem be- 
stimmten Planquadrat liegen. Infolgedessen sind 
heute ganz andere Anforderungen an die Ballistik 
zu stellen, als es vor dem Weltkriege nötig war. 
Mittlere Schußweiten von 10 km kann man heute 
auf rund 15 m genau, und Schußzeiten von 50 sec 


auf etwa 0.04 sec. genau bestimmen. Wenn wir 
heute in einem Heeresbericht lesen, daß unsere 


Ferngeschütze von der flandrischen Küste aus ein 
feindliches . Schiff aus einem Geleitzug heraus- 
reschossen haben, so ist ein Gutteil solcher Prä- 
zision sowohl den exakten Methoden der heutigen 
Ballistik wie auch der technischen Höhe unserer 
Schußwaffen zuzuschreiben. 

Berlin. R. Rothe. 244 

Dr. KARL GEY, Oberstudiendirektor in Leipzig, 
Dr. habil. HORST TEICHMANN, Dozent für Physik 
a. d. Techn. Hochschule, Dresden. Einführung 


in die Lehre vom Schuß (Ballistik). 
Math.-Phys. Bibl.. Reihe II, Abrisse aus dem 
Gebiete der Mathematik und der exakten Natur- 


wissenschaften.) III. verbesserte und erweiterte 


Aufl. 128 S. m. 69 Fig. u. 2 Tafeln. Leipzig und 
Berlin 1939, Verlage B. G. Teubner. Preis geb. 
3.60 M. 


Die neue Auflage des bekannten verdienstvollen 
Büchleins enthält einige Erweiterungen auf folgen- 
den Gebieten: Innere Ballistik, Schießen nach be- 
weglichen Zielen, Einfluß der Erddrehung auf die 
3estimmung der Flugbahn, Erweiterungen, die 
manchem Leser willkommen sein werden. Es darf 
jedoch nicht verhehlt werden, daß eine Beschei 
dung des Zieles an manchen Stellen des Buches, 
das doch immer einen schulmäßig elementaren Cha- 
rakter behalten wird, am Platze gewesen wäre. Dies 
trifft u.a. zu für die Bemerkungen. die sich über die 
Integration des gekoppelten Systems der beiden 
außerballistischen Hauptgleichungen auf S. 54 an 
diese anschließen. Hier hätte es genügt zu sagen, 
daß man heute sowohl rechnerische wie zeichne- 
rische Verfahren kennt. die Differentialzleichungen 
des außerballistischen Hauptproblems mit belie- 
bierer Genauigkeit zu integrieren, wobei das Luft- 
widerstandsgesetz beliebig durch Beobachtung ze- 
wonnen sein kann und sogar die Abhängigkeit der 
Luftdichte und der Schwerkraft von der Höhe be- 
rücksichtiget werden kann. Die Beweise auf Seite 58 
über die Existenz einer lotrechten Symptote und 


einer Grenzgeschwindigkeit bei jeder Schußbahn 
sind nicht fehlerfrei und sollten aus diesem Buche 


lieber ganz wegbleiben. Für ernsthafte Unter- 
suchungen über theoretische Ballistik reichen die 


Grundlagen dieses Büchleins nicht aus. Dagegen 

ist es für den Standpunkt der höheren Schule, ins- 

besondere auch für wehrtechnische Arbeitsgemein- 

schaften durchaus zu empfehlen. 
3erlin. 
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Dr. FR. REGER, Tachymetertafeln für 


neue zentesimale) Teilung. VI-+ 
490 S. Stuttgart 1940, I. B. Metzlersche Verlazs- 


buchhandlung. Preis geb. 19,40 M. 


Bei der tachymetrischen Punktbestimmung mit 
Hilfe des Fadenentfernungsmessers für topogra- 
phische Zwecke erhält man die horizontale Entfer- 
nune E und den Höhenunterschied h — die Be- 
zeichnungen werden in Anlehnung an den Ver- 
fasser gewählt auf Grund der beiden Näherungs- 
gleichungen 
: sin 2 
Dabei ist «a der Vertikalwinkel und D eine Rech- 
nungshilfsgröße, die man erhält au D=c-+kl, 
wobei € die Additionskonstante und k die Multi- 
plikationskonstante «des Instruments und ! der Lat- 
tenabschnitt ist. 

Für die Berechnung insbesondere von h auf 
Grund der obigen Gleichung gibt es verschiedene 
Hilfsmittel in Form von numerischen Tafeln. gra- 
phischen Tafeln und mechanischen Vorrichtungen. 

Bei der Herstellung einer numerischen Tafel für 
den Höhenunterschied " als Funktion von D und 
a bestehen zwei Möglichkeiten. In dem einen Fall 
wählt man — wie dies der Verfasser im Anschluß 
an die Tachymetertafeln für alte Teilung von W. 
Jordan getan hat — die Rechnungshilfsgröße D 
als Haupteingang in die Tafel. Im zweiten Fall 
wählt man den Vertikalwinkel «a als Haupteingang. 
Diese zweite Art der Anordnung verdient in bezug 
auf die Bequemlichkeit beim Einschalten zwischen 
Tafelwerten den Vorzug. 

Die Herstellung dieser neuen Tafel ist dadurch 
erforderlich geworden. daß nach einem Runderlaß 
des Reichs- und Preußischen Ministers des Innern 
bis zum 1. 4. 1945 die neue Teilung oder Zentesimal- 
teilung in Deutschland eingeführt wird. 

Stuttgart. P. Werkmeister. 223 


E=Deos?a und h=D 


[44 


Ferner sind bei der Schriftleitung folgende 
jücher eingegangen (ausführliche Besprechung 


bleibt vorbehalten): 


Dr. G. JOOS, Prof. an d. Univ. Göttingen, und 
Dr. TH. KALUZA, Prof. an d. Univ. Göttingen. 
Höhere Mathematik für den Prakti- 
ker, an Stelle einer 6. Aufl. des Lehrbuchs der 
Differential- und Integralreehnung von H. A. 
Lorentz. 2. verb. Aufl. XII+368.S. m. 
85 Abb. Leipzig 1940. Verlag Johann Ambrosius 
Barth. Preis geb. 24.50 M. 

Vorspannung im Eisenbetonbau. 
Grundlagen, Ziel, Zweck und Anwen- 
dung. Beiträge von Prof. Dr.-Ing. L. PISTOR, 
Dipl.-Ing. R. OPPERMANN, Dr.-Ing. W. PASSER, 
Oberhaurat Dr.-Ing. e. h. Dr. F. v. EMPERGER, 
mit einem Geleitwort von Prof. Dr.-Ing. A. KLEIN- 
LOGEL. 80 S. m. 68 Textabb. Berlin 1940, Verlag 
Wilhelm Ernst & Sohn. Preis brosch. 3,80 M. 


Prof. Dr. HUGO DINGLER, Max Planck und 


dieBegründungdersogenanntenmo- 


dernen theoretischen Physik. 3 S. 
Berlin 1939, Ahnenerbe-Stiftung- Verlag. Preis 
0,80 M. 
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Handbuch der Werkstofiprüfung. Herausgegeben 
unter besonderer Mitwirkung der Staatl. Material- 
prüfungsanstalten Deutschlands, der zuständigen 
Forschungsanstalten der Hochschulen, der Kaiser- 
Wilhelm-Gesellschaft und der Industrie sowie der 
Eidgenössischen Materialprüfungsanstalt Zürich. 
l. Bd.: Prüf- und Meßeinrichtungen, bearb. von 
R. Berthold, Berlin, A. Eichinger, Zürich- 
Düsseldorf, W. Ermlich, Berlin. G. Fiek, 
Berlin, L. Föppl. München, R. Glocker, 
Stuttgart, E. Lehr, Augsburg, E. Siebel, 


Stuttgart. O. Vaupel, Berlin, herausgeg. von 
Prof. Dr.-Ing. E. Siebel. Techn. Hochschule 


Stuttgart. XIV + 658 S. m. 763 Abb. Berlin 1940, 
Verlag Julius Springer. Preis geb. 69 M. 


Kalender der Technik 1941. 
Verlag, G.m.b.H. Preis 250 M. 


Berlin 1940. VDI- 


Hydromechanische Probleme des Schiffsantriebs, 
Teil II, Veröffentlichung der Vorträge, die anläß- 
lich des 25jährigen Bestehens der Hamburgischen 
Schiffbau-Versuchsanstalt am 14. Juni 1939 ge- 


halten wurden, herausgegeben von Dr.-Ing. G. 
Kempf. 236 S. m. 102 Abb. München 1940, 


Verlag R. Oldenbourg. Preis geb. 12,50 M. 


NACHRICHTEN 


Max Lagally zum 60. Geburtstag. 


Am 7. Januar 1941 hat Professor Max Lagally 
an der Technischen Hochschule in Dresden sein 
60. Lebensjahr vollendet. Das legt uns nahe, des 
wissenschaftlichen Werkes dieses vielseitigen Mathe- 
matikers zu zeedenken, der dem Fachgebiet dieser 
Zeitschrift so nahe steht wie wenig andere seiner 
Zunft. Der Familie eines Studienprofessors .der 
Mathematik mit weitem Gesichtskreis entsprossen, 
am humanistischen Gymnasium ausgebildet und 
über dessen Bannkreis hinaus früh angeregt, wid- 
mete er sich an der Münchener Universität dem 
Studium für das Lehramt der Mathematik und 
Physik. Es war die von G. Darboux und S. Lie 
ausgehende Richtung in der Mathematik, die ihn 
unter dem Einfluß von F. Lindemann und A. Voß 
besonders anzog. Seine Doktorarbeit über Flächen 
mit sphärischen Krümmuneslinien und eine Reihe 
diesen 


von weiteren Veröffentlichungen aus Ge- 
bieten zeugen davon. Die Annahme einer AÄssi- 
stentenstelle an der Technischen Hochschule in 


München erweiterte sichtlich sein geistiges Arbeits- 
gebiet. Da kam Lagally mit der Mechanik, mit der 
Vektorreehnung mit Ilydrodynamik, Flächenkrüm- 
mung und Photogrammetrie in greifbare Berüh- 


rung. Seine vorangegangene gründliche mathe- 
matische Ausbildung befähigte ihn bald, diesen 


Gebieten neue Gesichtspunkte abzugewinnen, zu 
schwierigen Aufgaben Stellung zu nehmen und ihre 
l.ösung anzubahnen. Es folgte die Habilitation an 
der Technischen Hochschule, deren Wirkung als- 
bald durch den Weltkrieg unterbrochen wurde, an 
dem Lagally als Photogrammeter teilnahm. Seine 
Habilitationsschrift über unendlich kleine isome- 
trische Verbiegungen einer Fläche mit höherer als 
erster Näherung birgt noch mancherlei fruchtbare 
Gedanken über das eigentliche Thema hinaus, die 
bei einer praktischen Verwertung der behandelten 
\ufgabe zur Geltung kommen könnten. Inzwischen 
war Lagally in das Lehramt an höheren Schulen 
eingetreten, was jedoch seine Fruchtbarkeit auf 
mathematischem Gebiete kaum beeinträchtigte. Be- 
sondere Hervorhebung verdient eine neue Seite der 
3etätigung Lagallys, die äußerlich mit seiner Assi- 


stentenzeit zusammenhing, innerlich aber durch 
persönliche Liebhaberei wohlvorbereitet war. 
lagally hatte sich aus dieser heraus von Haus 


aus begünstigt. beträchtliche Kenntnisse im Reiche 
der beschreibenden Naturwissenschaften verschafft, 
und da war es nur natürlich, daß er mit Freuden 
die Gelegenheit ergriff. sie bei den Gletscherunter- 
nehmungen des Alpenvereins zu verwerten, wobei 
ihm außerdem seine erhebliche Bergsteigerfähig- 
keit nützte. Wie immer beenügte er sich auch 
hier nicht mit einer beiläufigen Anteilnahme; er 
führte z. B. ganz selbständig die photogramm- 
metrische Konstruktion der Karte des Sulden- 
ferners nach dem damals noch üblichen mühseligen 
Verfahren durch. Als ihn dann 1920 ein Ruf als 
ordentlicher Professor an die Technische Hochschule 
in Dresden in ein weites Arbeitsfeld versetzte, fand 


er trotz der vielen neu aufgetretenen Verpflich- 
tungen die Zeit zur Entfaltung eines reichen 


wissenschaftlichen Schaffens. in dem die in seiner 
Lehrzeit in München gelegten Keime zu selbstän- 
diger Blüte und Reife gediehen. Ihre Frucht zeigt 
sich vornehmlich in drei Hauptwerken, die Lagallys 
unermüdliche Arbeit krönen. Es sind dies der zu- 
sammenfassende Bericht über die „Idealen Flüssig- 
keiten“ im Handbuch der Physik von H. Geiger 
und K. Scheel, sein Lehrbuch der Vektoranalysis 
und die „Mechanik und Thermodynamik des statio- 
nären Gletschers“. Allen ist neben tiefer Beherr- 
schung des Stoffes, Strenge des Ausdruckes und 
bohrendes Eindringen in die letzten Folgerungen 
der Lehre eigentümlich. Sie sind nicht volkstüm- 
lich, aber dankbar für ernstes Studium. Die Arbei- 
ten Lagallys über Gletscher haben seinen Namen 
in weite Kreise der Geographie und Geophysik ge- 
tragen und die diesbezügliche Forschung maß- 
gebend beeinflußt. Nichts wäre jedoch falscher. 
als wenn man angesichts dieses erfreulichen Er- 


folges Lagallys wissenschaftliche und lehrhafte 
Leistungen auf mathematischem Gebiete zurück- 
treten ließe. Ihr Reichtum und ihre Eigenart 


schützen sie vor solcher Schattenstellung. Überali 
hat Lagally Höchstes erstrebt und Bleibendes er- 
reicht. 

München. Seb. Finsterwalder. 248 

Zum 60. Geburtstag von E. Salkowski. 

Am 3. Februar 1941 vollendete Dr. Erich Sal- 
kowski, der Ordinarius für Geometrie an der 
Technischen Hochschule Berlin, sein 60. Lebens- 
jahr. Sein Name ist nicht nur mit der geometri- 
schen Forschung der letzten ‚Jahrzehnte, sondern 
darüber hinaus mit allen Bestrebungen verbunden, 
die auf eine Betonung der anschaulich-geometri- 
schen Denkweise, auf eine organische und gegen- 
seitig befruchtende Verbindung zwischen Geometrie 
und Anaiysis abzielen. Der geometrische Unter- 
richt in Schule und Hochschule hat die stärksten 
Anregungen durch die Arbeiten Salkowskis er- 
fahren, und wenn sich die alte „Darstellende Geo- 
inetrie*“ an unseren Hochschulen heute zu einer 
lebensnahen .„Konstruktiven Geometrie“ wandelt. 
bei der nicht mehr das .Darstellen“ vorhandener 
mathematischer oder technischer Raumformen, son- 
dern vielmehr ihre konstruktive Erfassung „in 
statu nascendi* im Vordergrund steht, so liegen 
die Wurzeln für diese Entwicklung bei den Ge- 
dankengängen Salkowskis. Der 60. Geburtstag 
wird und soll für den Jubilar kein Abschluß sein: 
mögen noch viele weitere Jahre rüstigen Schaffens 
folgen! 

Danzig. R U. Graf. 250 

Ende Januar verstarb der ord. Prof. für Ver- 
messungswesen und Photogramn:etrie an der Techn. 
Hochschule Dresden Dr.-Ing. R. Hugershoff. 

Am 3. Februar 1941 verstarb der ord. Prof. der 
Mathematik an der Techn. Hochschule Hannover 
Dr. G. Prange. 
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ZUSCHRIFTEN AN DEN HERAUSGEBER 


R. Sauer, Über Interpolation von Kurven- 
scharen mit Anwendung auf die Berechnung von 
Geschoßflugbahnen. Z. angew. Math. Mech. Bd. 20 
1940), S. 280 bis 284. 


Unter Anwendung von Affintransformationen ge- 
langt Herr Sauer zu einem einfachen Verfahren, 
das eine bequeme Linearinterpolation zwischen den 
Kurven einer Schar von Parameterkurven ge- 
stattet. Die Anfangswerte der affin transformierten 
Kurven enthalten bis zu den Ableitungen 2. Ord- 
nung den Parameter A nicht, sie haben also im 
Anfangspunkt für jeden Parameter drei unendlich 
benachbarte Punkte gemeinsam. Im besonderen 
führen die Sauerschen Betrachtungen bei der An- 
wendung auf das außenballistische Problem auf 
Formeln, die ich näher behandelt habe (Z. angew. 
Math. Mech. Bd. 19 (1939). S. 361 bis 371). 

Man kann nun die Genauigkeit bei der Interpola- 
tion in vielen Fällen noch erheblich steigern, indem 
man Transformationen einführt, welche die Para- 
meterkurven einander so annähert, daß sie im An- 


fangspunkt eine beliebige Anzahl von unendlich 
benachbarten Punkten gemeinsam haben. 

Gegeben sei eine Schar von Kurven 

SEEN: TEEN: 5.50’. DD 
mit den Ableitungen 

u=ult,A);v=dl(l,));d=Ult.i); 

2 {n) (mn) x (in) (m) . 
ser di... wer, di er; 


Die Anfangswerte seien Ug, Vg Ur. U U Ugerıı. 
U, ®„, die noch von A abhängen. Wir verlangen 
nun, daß allen Parameterkurven im Anfangspunkt 
rn +1) unendlich benachbarte Punkte gemeinsam 
seien. Dazu braucht dann die von Herrn Sauer 
abgeleitete Affintransformation nur auf das System 


(n-2) (n-2) _ +2 ei. . (n) (m) i 
u—_ ul); = u KA), Veul,f); , 
(2) 

(n—2) (n—2) r (n—1) (n—]) ? (n) (nm) “ 

= ED > N, ru 


angewandt zu werden. Man gelangt so schließlich 
zu folgender Transformation: 





‚un 
Pre... N’ Bd | 
” eos SD PR (9 . 
. in—ı 
Nn-COs w-Itang d- Un v' ' 
ef (n -»)! 
=» in—v 3 
ee —i2 5; (3). 
ne (n—»). | 
en in—ı 
n=n-sin u a M2 Um-7" m y)! 
„=? k 
von j} 
R 
= v Un r fin N 
Aw 
v‚—? zu v) g | 


Wendet man diese Formeln auf die außenballisti- 
schen Gleichungen an, so ergibt sich das Formel- 
system (4). Dabei wurde vorausgesetzt, daß im 
Anfangspunkt eine vierpunktige Berührung unter 
den Parameterkurven stattfindet (Parameter 
. 1+ sing Fe 
= = . 9 =Erhöhung). 


- / 


1]-sing-sing 


cosgy 
— Y-[n (v,) - 1]-sin 9; 
„= —[X-tangg — Y]-[n (v,) — 1 


wo 


na "rel 


sing =sing- |1- ! 
; n | g-|n (vo) 1] 


f’(v,) 
N (Vg) == Yg* . 
! ( ®o) 

(0, = Anfangsgeschwindigkeit: (A). 


4 = Temperaturgradient). 


Dabei ist: X=r -—- 1,008 9-8: 


Y =y- v,-sing +t = .f? 
Ferner: SV; H=V: &o c-f( u); 
> e!(v,)-c! (%) 
n=V; = 0: u= ®; 
FE G 
£ Yaeas „. -ct(v,).In(o,) — 1 
®n 





Bei Verwendung anderer Gesetze für die atmo- 
sphärischen Bedingungen ändert sich nur sing in 
diesen Formeln. 

Es sei noch bemerkt, daß es im Falle der Bal- 
listik nicht ratsam ist, noch höhere Berührungen 
zu verlangen, da die Anfangswerte der Ableitungen 
sehr bald unübersichtlich werden. An einigen 
numerischen Beispielen wurde festgestellt, daß der 
Interpolationsfehler durch diese neue Transforma- 
tion gegenüber der von Herrn Sauer und mir zuer:t 
beschriebenen um etwa die Hälfte sinkt. 
Schließlich sei noch eine Bemerkung zum Schluß- 
absatz der Arbeit von Herrn Sauer gestattet. Dort 
wird auf eine Verallgemeinerung der Affintransfor- 
mation hingewiesen, die bei ihrer Anwendung auf 
die Ballistik gestattet, aus einer Erdschußtafel die 
zugehörigen Flugbahnen (Luftschußtafel) abzuleiten. 
Ich habe auf diese Möglichkeit bereits im Schluß- 
absatz der schon zitierten Arbeit (Z. angew. Math. 
Mech. Bd. 19 (1939). S. 361 bis 371) hingewiesen und 
ihre Darstellung in einer gesonderten Abhandlung 
angekündigt. Obwohl diese aus bestimmten Grün- 
den bisher nicht erscheinen konnte, habe ich das 
Prinzip eines solchen Verfahrens in seiner «in- 
fachsten Form bereits für die Ballistik des Bomben- 
wurfes behandelt (Luftf.-Forsche. Bd. 17 (1940) 
Heft 7, S. 216). 


Berlin. Dr. Hermann Athen. 242 
Die Zuschrift des Herrn Dr. Athen hat Herrn 
Prof. Sauer vorgelegen. W. 


H. Athen: |Interpolationsverfahren zur Be- 
rechnung von Flugbahnscharen und ihre Verände- 
rung durch Variation der ballistischen Grundwerte. 
Z. angew. Math. Mech. Bd. 19 (1939), S. 361 bis 371. 


In dem Beitrag von H. Athen über ein Inter- 
polationsverfahren zur Berechnung von Flugbahn- 
scharen wird für den maximalen Fehler der zwei- 
gliedrigen Näherung die von der Flugzeit ?! unab- 
hängige Abschätzung in Prozent gegeben: 


a 2 (A, in? 
1 Wmaxl 4) = 1W- 11 K, +4, 
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Die Formel ergab sich durch Vergleich der Ent- 
wieklung 


v=w+twitwit-.- 
mit der harmonischen Reihe 

v—=w,-\1l+ . + +), 

2 3 

woraus folgt: 

Ws W, 11 

aka > et = A Be 
Tatsächlich ist aber Aw noch stark mit der 

Flugzeit ? veränderlich, und zwar wächst Aw für 


größere ! sehr schnell an. So betrug in einem Fall 
= MW m/s, 9 = 235°, %=45°) das Fehlerprozent 
für die Schußweite X der Erhöhung # =35° bei 
einer Flugzeit von t=65 Sekunden das Fünffache 
des obigen Ausdruckes. Die Schußweite unter- 
schied sich tatsächlich um 150 (m) von ihrem 
exakten Wert, wogegen jene Abschätzung einen 
Fehler von nur etwa 30 (m) erwarten ließ. Für 
größere Werte von $ und ? verliert demnach obige 
Formel ihren Wert auch als rohe Abschätzung. Zur 
Abschätzung des Fehlers von interpolierten Schuß- 
weiten ist sie bei größeren Erhöhungen (d. unt. 
Winkelgruppe) mit Vorsicht zu verwenden. 

Schließlich sei erwähnt, daß der Fehler der drei- 
gliedrigen Näherung 

v—w—+t w, + w,4? 


wegen der dann bestehenden Beziehungen: 


d wg: by dry Ag" Wy 
Sw— —(Ki,t ii, + 4) W 
dw, (++ 4) w, 
sich proportional zu 
8B=— (A — A) (A—4,) (A — 4,5) 


ergibt, wenn A,,/,. /, die den drei benutzten Aus- 
gangsbahnen ww’, w?, zugehörigen Para- 
meterwerte bedeuten. Ist etwa für /, der Fehler 
der dreigliedrigen Näherung bekannt, so ergibt sich 
der Fehler einer Bahn mit dem Parameter / zu 


w‘3 


4) 
44) 


A—a)(A—:))(4 . 
sul)— Den 2 -sw(h,): 
"RR - AH 

} +4 : i “ 

Setzt man 4, — -" I ?,so erreicht der Fehler für 


1 


6 (8 + Y3)- A. +ß@— v3) A.;] 


i= 


die maximalen entgegengesetzt gleichen Werte 


5 Was = + 0,05 (A, — A,)*-W,. 

Für das Verhältnis der maximalen Fehler der 
zwei- und der dreigliedrigen Näherung ergibt sich 
} . . u o fäl 

Gr eefähr- 
jieraus mit w, = 4;3 ungefähr: 
Öl max do—h 
„max —03-2 Er. 
A) u max ha + I 
hy—h . ’ 
wo; 77; in der unteren Winkelgruppe un- 
o— "ij 
günstig mit 0,3 angenommen werden kann. 
Setzt man den bei der zweigliedrigen Näherung 
äußerst möglichen Fehler mit 1% an, so ergibt 


sich der Fehler der dreigliedrigen Näherung hier- 
nach in der Größenordnung 10%. 


Essen. Walter Correll. 234 


Erwiderung. An einem numerischen Beispiel 
stellt Herr Correll fest, daß die von mir in der er- 
wähnten Arbeit durchgeführte Fehlerabschätzung 
zu kleine Werte liefert “und folgert daraus, daß sie 
für große Erhöhungen und Flugzeiten nicht zutriift. 

Die OR 
la Pr 
;,) e=zmw. .() 
r Fi) (a=z,y ( 
nimmt keine Rücksicht pi die Flugzeit, obwohl 
tatsächlich der Fehler stark von ihr. abhängt. Die Ab- 
schätzung (1) soll, wie die Art der Ableitung auch 
zeigt, nicht den Fehler seiner Größe nach liefern, 
sondern nur einen ganz rohen, der ersten Orientie- 


141%) = 10- 


rung dienenden Anhalt für den zu erwartenden Iı- 
terpolationsfehler ergeben. Für genauere Abschät- 


zungen trifft (1) nicht zu, und insofern ist der Ein- 
wand von Herrn Correll berechtigt. — Man weist 


nun leicht nach, daß bei der linearen Interpolation 
n—==2 


Vin, 
nz 
n—ü 


w (1,1) = w,(t) + 4-W, (f) 


die Reihe w (4,1) = w„(t) durch die Summe 


ersetzt wird, wobei stets 
n= x 
w (t > Dani, () und W(Ü)< w,(t) 
n==1 
wird. Andererseits ist offenbar die Abschätzung (1) 
wegen ihrer Herleitung aus der harmonischen Reihe 
gebunden an die Voraussetzung 


(2) 


Ian“) 2) ) 
er ‚diern 3) 
5 Sg Mr... 





Diese Ungleichung ist aber für große Flugzeiten in- 
folge der verschlechterten Konvergenz der Reihen 
nicht mehr erfüllt. Vielmehr gilt dann 
w'r E= ww? 
D) 0 r 
— Pr Er won >1. 
0 
Die Abschätzung (1) führt daher in den meisten 
Fällen zu besseren Ergebnissen, wenn man setzt 
Aa (%) = u- Aa (%) RE I, 
Herr Correll leitet ferner - Grund der Vor- 
aussetzungen für (1) ab, daß das Fehlerverhältnis 
der quadratischen zur linearen Annäherung in der 
Größenordnung 10? liegen müsse. Diese Tatsache 
besagt, daß die dreigliedrige Annäherung in allen 
Fällen Werte liefert, die innerhalb der geforderten 
Genauigkeit liegen. Diese Aussage läßt sich auch 
durch exaktere "Betrachtungen gewinnen und ng 
dann für den größten Fehler im Intervall a 
zu folgender Abse hätzung 





(Na — Da — a) 
da (%) = 10. [F7) u 3 | a) .([e 
a’ (na +3) + a”-(3na+ 1) 
(a=y,2,w) 
wo 
F + hy 2 05 0< N (5). 
Na TI = . (® = 0,5: genauer 10) ) 
r 1— w.4,/ ee ı 
oder 
/ 2 \o N @\r 
(0) (md): nyr (y y ) ‚(Tr 4) 
Dabei wurden näherungsweise die Kurven 
a (A, t)t=const. als Kurven 2. Ordnung in / ange- 
nommen. 


Die Formeln (4) und (5) liefern in allen Fällen ge- 


nügend genaue Werte, wie mehrere durchgerech- 
nete numerische Beispiele zeigten. — Ich werde 


gelegentlich einer wahrscheinlich demnächst abge- 


schlossenen größeren Arbeit ausführlicher auf Kon- 
vergenz- und Fehlerfragen u. a. eingehen. 
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